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ABSTRACT 


This thesis is concerned primarily with an invest- 
igation of the influence of eigenvectors and eigenvalues on 
linear system characteristics such as disturbance rejection 
and eigenvalue insensitivity to system parameter variations. 

Existing design procedures for eigenvector/eigen- 
value assignment are briefly reviewed and the results re- 
interpreted using range space restrictions. Such analysis 


ae order controllable system with 


readily shows that for an 
r inputs, complete state feedback allows arbitrary assign- 
ment of n eigenvalues and up to r elements of each of the n 
eigenvectors. 

The importance of eigenvectors in the design of 
control systems is illustrated by an investigation of the 
'disturbability' characteristics of multivariable systems. 

A linear system is defined as ‘undisturbable' with respect 
to.a particular input variable if the state or output var- 
iables of interest are not disturbed by arbitrary varia- 
tions in that input. Undisturbability is closely related 
to the system properties of uncontrollability and struc+ 
LUtal wunconterollabi lucy, sbutoica not dented BNecossary 
and sufficient conditions for k S$ r state variables to be 


h e e 
disturbance, require 


undisturbable with respect to the 4° 
that the system matrix (or equivalently the matrix of 
eigenvectors) be quasi-triangular with a kxn-k off- 


diagonal partition of zero elements and that the corres- 


ponding k elements of the apse column of the input 
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disturbance matrix be zero. These results together with 
existing eigenvector/eigenvalue assignment techniques provide 
a Simple, straightforward design procedure to produce 
undisturbability in linear systems. The problem of asymp- 
totic setpoint tracking in closed-loop undisturbable systems 
was also considered and it is shown that a solution to this 
problem is almost always possible. 

Feedback and feedforward controllers designed to pro- 
duce undisturbability were evaluated by experimental applic- 
ation to a computer controlled, pilot-plant evaporator. 

The results were superior to conventional multiloop control- 
lers, and comparable to controllers designed using optimal 
quadratic techniques. The design method was also applied 


: and 20th 


to ‘aie order models of two different binary dis- 
tillation columns and evaluated by digital simulation. 

For the class of systems which do not satisfy the 
necessary and sufficient conditions for undisturbability, a 
design procedure is proposed that uses state feedback control 
to: (1) minimise the effect of external disturbances on 
system outputs of interest and (ii) carry out arbitrary 
eigenvalue assignment in the closed-loop system. An exper- 
imental evaluation of this procedure on the pilot-plant 
evaporator demonstrated the practicality of this approach. 
Using duality this procedure is also applied to the design 
of full order observers to minimise the effect of unmeas- 


urable external disturbances on the state estimate of in- 


terest. 
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A set of simple, constructive conditions for achieving 
eigenvalue invariance to arbitrary variations in specified 
system parameters is derived in terms of the closed-loop 
system eigenvectors. An illustrative example demonstrates 
the superiority of this procedure over conventional pole 


placement techniques using unity-rank state feedback. 
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CHAPTER 1 


INTRODUCTION 


The basic objective of this thesis was to invest- 
igate the effect of external disturbances on specified 
system state, or outputvariables, and to develop design 
procedures for controllers that would eliminate or min- 
imise the effect of such disturbances. Most process 
systems operate in an environment where unknown exter- 
nal disturbances are almost always present and hence good 
regulatory control is usually a prime objective. Signif- 
icant theoretical advances have been made in recent years 
in the control of linear multivariables systems but, very 
few of these new techniques have been primarily concerned 
with the effect of external disturbances on the overall 
System response. In fact, most multivariable synthesis 
techniques do not make full use of information avail- 
able in the process model about the effect of disturb- 
ances, and many design techniques ignore disturbances 
altogether, or consider only special classes of dis- 
turbances such as impulses. Therefore, this study set 
out to look at state space models of systems and to 


determine under what conditions, if any, the state 


and/or output variables of such systems could be made 


invariant to arbitrary, unknown disturbances in one or 
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more system inputs. The main emphasis in the study was 


to be on the structural aspects of system models. 


La Relationship to Previous Work. 


Over the past 10 years a number of research pro- 
jects have been carried out in the Department of Chem- 
ical Engineering concerned with the development and 
experimental evaluation of modern control techniques 
such as optimal multivariable feedback control, (Newell 
[1], Fisher and Seborg [2]); adaptive control (Oliver 
[3],[2])+; multivariable frequency domain methods (Kuon 
[4]); and eigenvalue assignment by state feedback (Park 
[5],Topaloglu [6]). These investigations included ex- 
perimental applications on a computer controlled, pilot 
plant evaporator which showed that most of the modern 
control techniques gave significantly better performance 
than conventional single input, single output design 
techniques. However, some difficulty was experienced 
with the experimental implementation of state feedback 
controllers which had been designed to assign eigen- 
values to specified locations. In many cases the ex- 
perimental responses were unstable and it became clear 
that the design procedure should impose other conditions 


in addition to the location of specific eigenvalues. 
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Therefore, one of the concerns of this investigation was 
the identification of additional conditions that should, 
and could, be incorporated into pole-placement design 
procedures. 

Less nOt dufcicule. to show from aymodal 
analysis of the time domain response of the state var- 
iables of a system, that the eigenvectors as well as the 
eigenvalues, play a key role in determining the time 
domain response; system sensitivity; input/output 
decoupling; and system decomposition. In fact the 
eigenvalues of the closed-loop system serve only as 
exponential weighting factors on the terms making up 
the system time domain response and thus affect the form 
of the response. Thus, the initial direction taken in 
this study was to investigate the role of eigenvectors 
in multivariable feedback control systems and to deter- 
mine what specifications could be put on the eigenvectors 
during the design procedure. 

In addition to the theoretical studies and the 
development of a practical design procedure, the basic 
objectives of this thesis also included an experimental 
evaluation of the theoretical developments and design 


techniques. 
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1.2 Thesis Organisation 


The main contribution of this thesis is the 
development of the concept of system undisturbability and 
the determination of necessary and sufficient conditions 
for a system to display this property. In general terms, 
system state and/or output variables are said to be 
undisturbable by a particular input disturbance if they 
are unaffected by arbitrary variations in that input. 
These theoretical concepts are discussed in Chapter 3 
with the main emphasis being on the structural aspects 
of these results. The design procedure that was dev- 
eloped based on these theoretical results requires the 
use of eigenvectors/eigenvalue assignment techniques 
and hence this project also included the evaluation and 
interpretation of such techniques. 

The necessary and sufficient conditions for 
undisturbability are not overly restrictive and are 
usually not difficult to satisfy. However, for those 
cases where they cannot be satisfied an alternative 
design procedure was developed that minimises the 
effect of external disturbances on selected state and/or 
output variables and also permits arbitrary eigenvalue 
assignment. The investigations of the structural aspects 
of the results for disturbance localisation also led to 


conditions under which specific closed-loop eigenvalues 
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could be made invariant to arbitrary and unknown vari- 
ations in system parameters. The theoretical results 
and design procedures developed for disturbance local- 
isation, disturbance minimisation and eigenvalue invar- 
iance were all evaluated by numerical simulation and by 
experimental application to the pilot-plant evaporator. 
The order of the chapters in this thesis has been de- 
Signed to present the above results in a convenient order. 
Chapter 3 contains the main theoretical results and is 

a prerequisite to the full understanding of all the fol- 
lowing chapters. However, the other chapters are 
relatively independent of each other and may be read 

in any order. 

Chapter 2 presents a brief review of the main 
results on existing eigenvector/eigenvalue assignment 
techniques and conditions which specify how many and 
which elements of the eigenvectors can be arbitrarily 
assigned are given in terms of range-space restrictions. 
The implications of controllability on eigenvector 
assignment is also discussed. 

In Chapter 3, the concept of undisturbability 
is formally defined and related to uncontrollability 
andustructural controllability. Here the necessary 
and sufficient conditions for a system to have undis- 


turbable state or output variables are expressed in 
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terms of the structure of the coefficient matrices of 

the state space model and the structure of the eigen- 
vector matrix. These results are then used in conjunc- 
tion with the eigenvector/eigenvalue assignment algorithm 
discussed in Chapter 2 to arrive at a design procedure 
for disturbance localisation. 

The experimental evaluation of multivariable con- 
trollers designed to produce undisturbability is 
presented in Chapter 4. The performance of the dis- 
turbance localisation controllers is compared with the 
performance of conventional multiloop controllers and 
with those designed using optimal control techniques. 

The design procedure to achieve undisturbability is 


ene 2 


also applied to 20 order state space models 

of two different distillation columns. A discussion of 
the combination of proportional state feedback with 
integral feedback of specific state variables of interest 
is also included in this chapter. 

In Chapter 5, a method is presented for design- 
ing state feedback controllers which assign eigenvalues 
of the closed-loop system to specified positions, and 
also minimise the effect of external disturbances on 
state variables of interest. Results from the exper- 


imental application of this method to the double-effect 


evaporator are also included. An application of the 
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method to the design of full order observers for systems 
with unmeasurable disturbances, is also included. 

In Chapter 6, a method is presented for design- 
ing constant multivariable feedback controllers that 
make selected closed-loop eigenvalues invariant to un- 
known perturbations of arbitrary magnitude in system 
parameters. The sufficient conditions for eigenvalue 
invariance are expressed in terms of the structure of 
the closed-loop system matrices. 

Some overall conclusions and recommendations for 
future work are presented in Chapter 7. 

Appendix A is devoted to the discussion of sel- 
ected geometric concepts which are used in the main 
body of the thesis. Appendix B outlines in detail the 
difficulties associated with eigenvector assignment 
using output feedback. Appendix C contains a descrip- 
tion of the pilot-plant double effect evaporator; the 
associated instrumentation; and the computer equipment. 
The relationship between the concepts of invariant 
Subspaces, invariant zeros and parametric sensitivity 
is discussed in Appendix D. A new algorithm to compute 
the invariant zeros of a system is also presented in 


this Appendix. 
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CHAPTER 2 


EIGENVALUE/ELGENVECTOR ASSIGNMENT FOR 


MULTIVARIABLE SYSTEMS 


Abstract 


Existing design procedures for eigenvalue/ 
eigenvector assignment are briefly reviewed. The prob- 
lem of eigenvalue/eigenvector assignment using state or 
Output feedback is considered using range space restric- 
tions. This readily shows how many and which elements 
of the closed-loop system eigenvector or reciprocal 
basis vector can be specified. The implications of 
controllability in eigenvector assignment is also dis- 


cussed. 
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Since 1962 when Rosenbrock [1] first introduced 
modal control as a possible design aid in the control of 
large chemical plants, considerable attention has been 
paid to control system design by assigning the eigen- 
values of the closed-loop system. Eigenvalue assignment 
techniques have been widely used in the control system 
design literature to improve system stability and the 
speed of response of interacting processes. A wide var- 
iety of different algorithms/methods exist for eigenvalue 
placement via state feedback. Most of the well-known 
algorithms and results in this area, prior to and 
around 1973, have been Smears in surveys by Topaloglu 
Plmandarark 13]\. 

It is well known that for a linear, controllable 
system, eigenvalues can be assigned to arbitrary positions 
by uSing state feedback. Conventionally this has been 
achieved either by using unity-rank state feedback 
(Simultaneous design) or by non-unity rank state feedback 
(recursive design). However, for multi-input systems the 
controller is not unique since more than one feedback 
controller can be designed to achieve the desired eigen- 
value placement. Different feedback controllers result 


in different system performance. Closed-loop system 
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eigenvectors are equally important in determining system 
performance. This becomes evident at once from a modal 
analysis of the closed-loop system. 

For an ye order controllable system with r inputs, 
Only n degrees of freedom are required to assign all n 
eigenvalues to arbitrary locations. The remaining 
nx(r-1) degrees of freedom available in the design of a 
state feedback controller can be utilized to assign com- 
ponents of the eigenvectors. However, very few algorithms/ 
methods are available at the present time that fully 
utilise all degrees of freedom available to arbitrarily 
assign some elements of the eigenvectors as well as all 
of the closed-loop system eigenvalues. 

Srinathkumar and Rhoten [4] have proposed an 
algorithm for the computation of a state feedback matrix 
that assigns closed-loop eigenvalues as well as certain 
elements of the closed-loop eigenvectors. The relation- 
ship between their algorithm and the matrix pseudoinverse 
approach was shown as part of this work (cf. Section 2.2 
and onan et al. [5]). fn 15], application, of simple 
geometric concepts also showed how many and which elements 
of the closed-loop eigenvectors can be assigned arbit- 
rarily. More recently Moore [6] has provided a simple, 
but general, mathematical characterization of the class 


o£ all closed-loop eigenvector sets which can be obtained 
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for a given set of distinct closed-loop eigenvalues. 
Using a minimisation procedure, Ramar and Gourishankar 
[7] have also proposed a method whereby all the closed- 
loop system eigenvalues and some or all of the elements 
of the closed-loop eigenvector matrix can be assigned. 
In deriving an alternative proof for conditions for pole 
asSignability, Kimura [8] has also shown how eigenvalues 
and eigenvectors can be assigned using state feedback 
Control. 

Bigenvalue/eigenvector assignment procedures are 
discussed here because they form an important part of the 
design procedures developed in later chapters. The 
present chapter is organized as follows. The problem of 
eigenvalue/eigenvector assignment using state feedback 
is considered in Section 2.2. The implications of con- 
trollability in eigenvector assignment is discussed in 
Section 2.3. The difficulties associated with eigenvalue/ 
eigenvector assignment using incomplete state feedback is 
considered in Section 2.4 followed by some concluding 


remarks in Section 2.5. 


2.2 Assignment of Closed-loop Eigenvalues and 
Eigenvectors 


Consider the en order completely controllable, 


linear, time-invariant system. 
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with r inputs, m outputs and a state feedback control 


law 
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The resulting closed-loop system matrix, H, can be ex- 


pressed as 


ae 

II 
YP 
+ 


BK = WJV 2rs)) 


where J is the Jordan canonical form defined by the desired 
closed-loop eigenvalues; W is the closed-loop system 


eigenvector matrix; and vew i 


is the matrix of reciprocal 
basis vectors. The following relationships apply when 
the closed-loop system matrix, H or WJV, is fully speci- 


fied 


i) A control matrix K which exactly satisfies eqn. 


(2.3) existorur sand only if 


Range (WJV-A) © Range B (2.4) 
Liye Tiecondition W204) is ‘satistied, thena) control 
matrix K, which exactly satisfies eqn (2.3), can 
be computed from 


je B. (WJIV-A) CONS) 
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* 
where B_ denotes the pseudoinverse of B given by 
‘s Ne ca Meg 
B = (B’B) ~B’ and T denotes the transpose of a 


matrix. 


The manner in which condition (2.4) determines 
how many, and which, elements of W can be specified 
arbitrarily is illustrated by writing eqn. (2.3) in the 


partitioned form: 
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where Adie Bir Kier Wii Jy and Vii are rxr matrices, Aoor 
Woor J5 and Vo are (n-r)x(n-r) matrices and B is assumed 


to have full rank. From the above eqn. (2.6) it is not 
difficult to see that the rxr entries in Wii and the 


rx(n-r) entries in W can be chosen to be completely 


2 

arbitrary if and only if Range (B,) = Re a1. LE and 

1 is nonsingular. Implicit in the above condition 
1 


is the requirement, since V = W , that W be nonsingular. 


only 22 8B 


This in turn implies that the closed-loop system should 
have distinct eigenvalues. Such a restriction, however, 
is not limiting. Note that specifying the first r elements 


of each of the n eigenvectors requires only rxn-n degrees 
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of freedom, i.e. rxn-n elements of K, providing that at 
least one of the r specified values in each eigenvector 
1S nonzero. Once Wii and Wi> 


are computed such that condition (2.4) is satisfied. 


are so chosen, Woy and Woo 


One method of calculating Woy and Wo is to use 
the algorithm presented by Srinathkumar and Rhoten [4]. 
They start with an equation identical to eqn. (2.6) and 
by algebraic manipulation determine the degrees of freedom 
available for specifying elements of W. The algorithm 
uses state feedback control to assign the closed-loop 
eigenvalues as well as elements of the eigenvector matrix. 
The derivation of the algorithm requires some tedious 


algebraic manipulations, as shown in Appendix B. The 


two main equations of the algorithm are: 


W144 - PW, = TWi4 ~ SW, 49, (O57) 
WooJ5 ~ PWo5 = TW, > + SW, 55> (2:8) 
where 
a = . ¥, & fi 
= B.B, Pee ele Ao, ae and P A, SA 


When the right hand terms in eqns. (2.7) and (2.8) are 
known, the algorithm represents a set of n@-nr linear 
equations in as many unknowns. It can be shown that, 
dats Wii and Wio 
calculated by the above algorithm, then the resulting W 


are chosen arbitrarily and Woy and Wo5 are 
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will satisfy condition (2.4). It therefore follows that 
K can be computed using the matrix pseudoinverse approach 
as in eqn. (2.5). Note that an identical K can also be 
evaluated from eqn. (2.6). There is a computational 
advantage in uSing the pseudoinverse approach, that is it 
is applicable generally. For example, when condition 
(2.4) is not satisfied exactly due to numerical round-off 
errors, no control matrix can satisfy eqn. (2.3). However, 
the control matrix computed from eqn. (2.5), using the 
pseudoinverse of B, satisfies eqn. (2.3) in the least- 
squares sense. 

Moore [6] has characterised the class of all 
closed-loop eigenvector sets which can be obtained for a 
given set of distinct closed-loop eigenvalues using 
state feedback. This is a more general result that does 


not require B, to be nonsingular. His result requires 
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ay) w, € span (N, ) 
uh 
ii) for i€n vectors Ww, are linearly independent 


: n 
vectors in ¢€”, and 


iii) Ww, = w. whenever A, = Xs where denotes a com- 


plex conjugate. 
In [6] Moore has made some remarks regarding controllabil- 
ity and eigenvector assignment. While equivalence of 
controllability and pole assignability is well known, the 
relationship between controllability and eigenvector 
assignment, however, is not so well known. The following 
section contains a discussion of controllability consid- 
erations in eigenvector assignment based on the algorithm 


in, (4). 


2090 eOntrol lability Considerations in Eigenvector 


Assignment 


Consider an wee order, r-input uncontrollable 
system. Let the Kalman canonical form of such a system 
be characterised by the pair (Anes where matrices A and 


B in their partitioned form are: 
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From eqn. (2.10) it is clear that the system under con- 
sideration is uncontrollable and that the eigenvalues of 
ASS correspond to the uncontrollable modes of the system. 
Application of the eigenvector/eigenvalue assignment 


algorithm, i.e. eqns. (2.7) and (2.8) from [4] to the 


pair (A,B) as shown in (2.10) leads to the following re- 


sults. 
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where matrices S,T and P have been defined in eqn. (2.8). 
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assigning eigenvectors corresponding to uncontrollable 
modes, eqn. (2.8) will now be considered. For an arbit- 


rary choice of Wy and J. given by 


EGE; 


where Wi and J5 have been compatibly partitioned with 


respect to matrices T and S, the appropriate elements of 


Wo5 are given by: 


i 1 ma 2 
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(Qi) 
Eqn. (2.12) clearly implies that diag (J5) = o(A5,) and 


that Wo is the eigenvector matrix of Ato: Since 


distinct eigenvalues have been assumed to have been 


3 
22 


is intuitively expected. This leaves two matrix linear 


specified, eqn. (2.12) also implies that W = 0, which 


equations to be solved for Wo and Ws which can be 


solved for any arbitrary value of Wi2° (Note that if 
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one or more Vectors in W are selected to be zero then 
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the corresponding eigenvalue in J. should be the approp- 


2 


riate eigenvalue of P). The conclusion to be drawn 
from the above analysis is that the r elements of each of 


the eigenvectors i.e. We and Wr corresponding to the 


uncontrollable modes can be arbitrarily assigned, while 


the corresponding £ eigenvalues remain invariant under 


eC 
Zz 


matrix block to be identical to the eigenvalues of ae 


feedback and thus they should be specified in the J 


The main point to note here is that while con- 
trollability is necessary and sufficient for pole 
assignment it is not necessary for eigenvector assign- 


ment (i.e. for assigning r elements in each eigenvector). 


2.4 Eigenvalue/Eigenvector Assignment Using Output 


Feedback 


Consider eqn. (2.1) in conjunction with the 


output feedback control law, 


(tt) Se=yehyGt), -=- ak ext) (2213) 


the resulting closed-loop system matrix is given by 


(A + BFC) = WJV (2.14) 
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For a given closed-loop system matrix, WJV, the approp- 
riate F exists if and only if condition (2.4) and the 


following condition are satisfied: 
ih T 
Range (WJV - A)” © Range C Zien >) 


conditions (2.4) and (2.15) impose significant restric- 
tions on the specification of the elements of the closed- 
loop system eigenvectors and reciprocal basis vectors, 
and hence in the general case a number of difficulties 
arise in the design of F. These are discussed at the 
end of this section. Normally, n of the rxm available 
degrees of freedom (assuming rxm > n) are used to assign 
the closed-loop eigenvalues, leaving only rxm-n degrees 
of freedom available for assigning elements of the 
eigenvectors or of the reciprocal basis vectors. fThus, 
in the general case, it is probably not worthwhile to 
attempt to assign elements of the eigenvectors or recip- 
rocal basis vectors unless rxm>>n. 

The output-feedback design problem for aie 
special case of an n-input, n-state, m-output (m < n) 
system can be readily handled because condition (2.4) 
is always satisfied for such systems, and hence only 
condition (2.15) need be considered for the design pro- 
cedure. The algorithm for this case is the dual (or 


adjoint) of the one for state feedback control and can 
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be obtained by using algebraic manipulations on the 
dual of the system of eqn. (2.1). Alternatively, it 
can also be derived by considering the following part- 


itioned form of eqn. (2.14) 
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All ! Ajo By F Ci | C.| 
eae La ee 4 
1 
2 ieee =e 
1 a ae \ 
Wa ! Neely ! Soe a. | vie 
EP EL tae 11 ae + mee |) i imine pone (2-516) 
i i 1 
| i i 
BN ooh eee ra: dp sl oi tues 
Where Ady Cyr Wii Jy and Vii are now mxm matrices; 
Ayo: Woo Jo and Vo5 are (n-m) x (n-m) matrices; By is 
an mxr matrix; F is an rxm matrix; and C, and B are 


al 


assumed nonsingular, obtained, if necessary, by reorder- 


ing of the state variables. Bearing in mind that 
Range oe = R™ 


and the (n-m)xm entries in V 


, it follows that the mxm entries in Vii 

5, Can be chosen arbitrarily, 
subject only to the requirement, since W = va that 

V be nonsingular. Once Vi and Vo4 are chosen, Vi2 and 
V55 can be computed such that condition (2.15) is 
satisfied. Algebraic manipulation of eqn. (2.16) to- 
gether with the use of the identity VW = I leads to the 


following two equations for calculating V and Vy? 
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JiVi9 ~ VR = J4V4 32 + VM Coe a7) 
JoV55 - VooR =) J5V512 + Vo4M (2506) 

where Gh crite M=A =A..0, RISA ee ee 
yero! 1 sl ae 22 a Goa 


The derivation of these equations is discussed in detail 
in Appendix B.. ~ Eqns.» (2.17) and (2.178) are linear ian 
and V 


Vio and Vo9 and permit Jy J to be speci- 


uae bp 21 
fied. It can also be shown that V and V obtained 


2 22 
from eqns. (2.17)"and? (2.138) satisfy ‘condirtzon “(2 215) 
and hence this is a practical approach for designing 
output feedback controllers for n-input, n-state 
systems. 
The output feedback matrix can also be computed 
using the pseudoinverse approach. If condition (2.15) 


holds the following equation leads to the same output 


feedback matrix as eqn. (2.16): 


EL) =a Boke eC ee Bs) Ce (2.19) 


where Cc = cliec!) + 


denotes the pseudoinverse 
The restrictions that: conditions (2.4) and (2.15) 
impose on the general output-feedback design problem for 


r<n for eigenvalue/eigenvector assignment can now be 
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realized by considering the following steps in the design 
procedure, further details of which are discussed in 


Appendix B. 


i), Arbitrary Choices or V Voy" W and W is 


ial 11 2 
made subject to the constraint that the identity: 


Wi 11 + Wi2Vi2 = I (from WV = I) be preserved. 


Li) After Vii" V 


and V 


21! Wii and Wid have been specified, 


V can be calculated from eqns. (2.17) 
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and. (2.18)) (to, Ssacisty conditvon (2515) and Woy 


and Woo calculated from the dual form of egns. 


(2.17) and (2.18) te. eqns. (2.7) cand (2.8). 
iii) With V and W determined, a check must be made to 

determine if the identity WV = I holds. If not, 

the entire procedure is repeated with a differ- 


Ww and Wen. 
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ent choice of V 


ley” Dale’ 


Implicit in the above mentioned procedure is the 
assumption that the eigenvalues of the closed-loop 
system can be arbitrarily assigned. This usually implies 
that (r+m-1)2n, Topaloglu and Seborg [9]. If arbitrary 
assignment of all eigenvalues is not assumed then the 
algorithm would become prohibitively complex. MThere- 
fore in general, unless rxm>>n (as remarked earlier) 
and (r+m-1)2n, it is probably not worthwhile to attempt 
to carry out eigenvalue/eigenvector assignment using out- 


put feedback. 
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The problem of eigenvalue/eigenvector assignment 


using incomplete state feedback can also be considered 


th 


by noting that the i column, Ks of the state feedback 


matrix, K, can be expressed as a function of h and ass 


the Age columns of H and A respectively, by: 


ke, = Be a 2.20 


It is easily observed from eqn. (2.20) that if measure- 


ment of Xs is not available then one possible way of 


making kK, = 0 is to have hy = a.. However, it is not 
always easy to choose Wii and Wi> such that the resulting 
he = (WIV) 5 = a,. There are some exceptions to this 


case and these will be discussed in the application of 
the eigenvalue/eigenvector assignment algorithm to the 
distillation column example in Chapter 4, where eqn. 

(2.20) is used to advantage in reducing the number of 


state variables required for feedback control. 


2.5 Conclusions 


The problem of eigenvalue/eigenvector assign- 
ment using state or output feedback has been 
considered using range space restrictions. For the 
state feedback case such an analysis readily shows how 
many and which elements of the closed-loop system 


eigenvector can be specified. The implications of 
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controllability on eigenvector assignment are discussed 
and it is shown that in using the algorithm in [4], 
controllability is not necessary for eigenvector assign- 
ment. An eigenvalue/eigenvector assignment method for 


output feedback control of aie 


order, n-input, m-output 
(mkn) systems is also presented. The restricted set of 
range space conditions and the associated difficulties 
‘ : : th 
in eigenvector/eigenvalue assignment for a general (n 
order, r-input, m-output) system using output feedback 


-are discussed. 
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CHAPTER 3 


UNDISTURBABILITY OF LINEAR MULTIVARIABLE SYSTEMS 


Abstract 


The concept of undisturbability, is formally 
defined and related to uncontrollability and structural 
controllability. Necessary and sufficient conditions for 
a system to have undisturbable state or output variables 
are expressed in terms of the structure of the coefficient 
Matrices in the state-space model and the structure of the 
eigenvector matrix. These results and the eigenvector/ 
eigenvalue assignment algorithm,as discussed in Chapter 2, 
provide the basis of a design procedure for synthesizing 
multivariable controllers which achieve disturbance loc- 


alisation. 
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3.1 Introduction 


Practical control systems must operate in envir- 
onments where unknown external disturbances are invariably 
present. Some typical examples of such disturbances in 
control problems are: load fluctuations in power systems; 
ambient temperature changes and unexpected variations in 
feed conditions in distillation columns, chemical reactors, 
etc.; wind gusts and turbulence acting on aircrafts, etc.. 
Despite the inevitability of such external disturbances, 
Synthesis techniques for multivariable control systems 
usually ignore them or consider only a special class of 
disturbances such as impulses. Consequently, the result- 
ing multivariable controllers may not perform satisfact- 
Orily in the presence of other types of disturbances, 
particularly sustained or persistently fluctuating dis- 
turbances. 

In this chapter a linear system is said to be 
"undisturbable" with respect to a particular input var- 
iable if the state or output variable(s) of interest are 
not disturbed by arbitrary variations in that input. The 
property of undisturbability is of considerable peer est 
to control engineers because it implies perfect regulatory 
control in the sense that the input disturbance does not 


affect the system variable(s) of interest in any way. 
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In the present investigation, the problem of con- 
trolling systems subjected to persistent, unknown external 
disturbances is considered from a system structure point 
of view. Necessary and sufficient conditions for an open- 
loop or closed-loop system to have undisturbable state or 
Output variables are expressed in terms of the structure 
of the state space matrices. An alternative set of nec- 
essary and sufficient conditions for undisturbability are 
also derived in terms of the closed-loop system eigenvec- 
tors. These latter conditions and the eigenvector/ 
eigenvalue assignment algorithm considered in Chapter 2 
provide a design basis for synthesizing multivariable 
controllers which produce disturbance localisation. 

These results are equivalent to existing necessary 
and sufficient conditions on disturbance localisation de- 
rived from a geometrical characterization of the problem. 
However, in contrast to previously known geometric re- 
Sults, the necessary and sufficient conditions developed 
here are easier and more straightforward to use, both 
from a computational, as well as, applications point of 
view. In contrast, the geometric results that already 
exist do not lend themselves easily to computer imple- 
mentation; furthermore, the computational operations 
that they require are subject to the hazards of numerical 


ill-conditioning. 
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3.2 Previous Work 


The problem of controlling systems subjected to 
persistent, unknown, external disturbances, i.e. the 
regulatory control problem has been one of the oldest and 
most fundamental problems of control. The regulatory 
control problem has been widely researched. There is a 
variety of terms used in the literature to describe this 
and related problems. The main terms that will appear in 
the discussion below are: dtsturbance loecalisatton, dts- 
turbance rejectton and system invariance. 

Rutman and Epelman [1] were amongst the first 
researchers to carry out a formal treatment of this sub- 
ject in their study of the effect of external disturb- 
ances on the output of a system. They defined a system 
output to be invariant with respect to a disturbance if 
the output was completely independent of the disturbance, 
and they obtained necessary and sufficient conditions for 
a Single-input, single-output system to be invariant. 
Preminger and Rootenberg [2] defined a system to be 
Steady state invariant or completely invariant to a 
Class of disturbances depending on whether the steady 
state or the entire effect of the disturbances is elimin- 
ated. They also discussed, heuristically, the design of 
compensators to attain complete invariance. Using vari- 


ational methods Rozoner [3] also obtained necessary and 
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sufficient conditions for a single-input, single-output 
system to be invariant. By using Rozoner's result, Wang 
[4] related the concept of invariance to uncontrollability 
and unobservability. Later Cruz and Perkins [5,6] gener- 
alized the results on invariance to the multivariable case 
and also obtained results relating invariance to unobserv- 
aod Lien. 

As early as 1966, C.D. Johnson [7] conceived the 
idea of invariant hyperplanes and discussed their role in 
state space theory. Geometrically, a hyperplane passing 
through the origin is a subspace. Johnson's results were 
actually based on a series of earlier Soviet publications. 
It was a geometric concept related to this, namely that 
of (A,B)-invariant subspaces that finally led to general 
results on disturbance localisation. Ina later public- 
ation, Johnson [8] considered regulation of a linear 
System subjected to constant disturbances using an optimal 
control formulation. In a series of papers [9,10], he 
extended these results to a class of unknown disturbances 
that satisfy a specified linear differential equation, or 
that can be characterized by certain waveform types or 
modes. 

Basile and Marro [11] and Wonham and Morse [12] 


independently introduced the concept of (A,B)-invariant 
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subspace and (A,B)-controllability subspace. These con- 
cepts played a basic role in the derivation of a general 
result on disturbance localisation by Wonham and Morse 
[12]. The main result of Wonham and Morse a a sufficient 
and necessary condition that is also a synthesis result 
for disturbance localisation using state feedback control. 
By considering the problem of zeroing the disturbance 
transfer function in a linear system, Bhattacharyya [13] 
has also derived a set of equivalent geometric necessary 
and sufficient conditions for disturbance localisation. 
Using an algebraic approach Shah et al. [14] - (cf. 
Section 3.5) have also derived a set of necessary and 
sufficient conditions for disturbance localisation. In 
contrast to the synthesis result in [12] the results in 
[13,14] are structural results, and ae would be expected 
they are weaker than the conditions given in [12]. In 
an earlier investigation McLane and Davison [15] have 
also reported algebraic conditions which are sufficient 
for disturbance localisation. An extensive treatment 
of the regulator problem based on an algebraic approach 
can also be found in Silverman and Payne [16]. 

In a more recent paper Fabian and Wonham [17] have 
reported necessary and sufficient conditions for simult- 
aneous disturbance localisation and input-output decoup- 


ling. They have also proved conditions for disturbance 
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rejection and closed-loop stability. Using a geometric 
approach, Chang and Rhodes [18] have also derived neces- 
sary and sufficient conditions for simultaneous disturb- 
ance localisation and pole assignment or closed-loop 
stability. 

In [13] Bhattacharyya derived solvability condi- 
tions for disturbance localisation using feedback plus 
feedforward control. Shah [19] (cf. Section 3.6) has 
derived necessary and sufficient conditions for the 
existence of a feedforward controller that would localise 
the effect of external disturbances. More recently 
Bhattacharyya [20] has described a design procedure for 
designing a disturbance localisation controller that 
uses dynamic feedback and feedforward compensators. 

An alternative design strategy is to design a 
controller which reduces the effects of disturbances on 
selected state or output variables. Willems [21] and 
Muller and Luckel [22] have recently proposed synthesis 
procedures for the design of state feedback controllers 
that minimise or reduce the effect of disturbances on 
Outputs of interest. Chapter 5 presents a Similar 
procedure that minimises the effects of disturbances on 
Outputs of interests and also allows arbitrary pole 
assignment. 


By a geometrical characterization of the regulator 
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problem, Bhattacharyya et al. [23] have derived conditions 
under which a zero output of a system is maintained at 
zero, and any non-zero output due to non-zero initial 
conditions or impulse-type disturbances is made to approach 
zero. Davison [24] and Grasselli and Nicolo [25] have 

also considered the design of output feedback controllers 
to accommodate a class of disturbances that can be des- 
cribed by a set of differential equations. Johnson's 

work [8,9,10] whichwas discussed earlier, also falls in 


this category except that he requires state feedback. 


3.3 Statement of the Problem 
Consider the linear time-invariant dynamic system 


x(t) = Ax(t) + Bu(t) + DE&(t) 


(331) 


y (t) Cxz{t) 


where x € R” is the state, u é€ R is the pM aenehean © 1s, RI is 
the disturbance and y € R" is the output; A, B; Cvand D 
are constant matrices of appropriate dimensions. Ifa 


feedback-feedforward control system of the form 
u(t) = Kx(t) + KF e(t) (2) 
is employed, the resulting closed-loop system is 


x(t) = Hx(t) + LE(t) (3.3) 
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where H and L are defined by 


+a 
Il 


Ar aBe (3.4) 


and 


ens yea yaaa (3.5) 


The disturbance localisation problem can now be posed as 
follows: For the closed-loop system in eqn. (3.3), is it 
possible to decouple the effect of §(t) from the system 
Outputs (or some of the state variables), for all t > 0 


FE 


by an appropriate choice of K (or K and K if feedforward 


control is possible)? 


3.4 Definitions 


The time response, x(t), of the system in egn. 


(3.3) can be expressed as the sum of two terms: seas Clee 


E 


due to the initial state x(0); and x’(t), due to the 


disturbances, i.e. 


(ES oe k(t) (3.6) 
where x° (t) = exp(Ht)x(0) and 
E tc 
xe Ct) a | exp{H(t-t) }L&(t) dt (35.7) 
fe) 


In this investigation, the term, AAG is the primary 
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concern. For simplicity, the disturbance rejection 
characteristics of system (3.3) will be considered in 
terms of a single disturbance, oh Veet! Aaa a) = 
However, the results that are obtained can be easily 
geheralized- to all’q disturbances.” It is clear from 

eqn. (3.7) that the disturbance rejection characteristics 
of system (3.3) with respect to disturbance, emt are 
characterised by the pair hier for the state variables, 


=k th 


and by (Hy L505) fOrmthe rl output where ee is the j 


column of L and ce; is the ine LFOWHOLsC. 


Definttiton ts: “State Undisturbabiiity 


The gen 


Stare varvable, Les of a system charactertsed by 
the pair (H,£ 4) ts satd to be undtsturbable wtth respect 
to dtsturpance, b 4) Lp ep Or MaYrpDTerary Bytes and for all 


La>S). ne, the ae element of ee eattsyies: 


we (t) a (3.8) 


Definition 2: Output Undtsturbabiltity 


the 1° output, Y,, of a system charactertsed by iad Aare 


where ce; ts the pine HOW OfeG ts Satd to pe Uundverurpaple 


with respect to disturbance, Eas tf for arbitrary E y(t). 


th 


and For@art + > 0, Hoe Lhe 4 element of Hee SQttsites 


ye(t) = 0 (3.9) 
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When condition (3.6) (or (3.9)) 1s not satisfied 
for some ee) # 0 and some t > 0, then the sls state 
variable (or output Secs is said to be disturbable. 
When all elements of thesstate (Or output) vector are 
undisturbable or disturbable, the modifier complete is 
used; for example, as in complete state (or output) dis- 
turbability. Complete state undisturbability implies that 
the disturbance has no effect on the system, as would be 
the case with a "perfect" control policy. 

The concept of undisturbability is closely related 


to the concepts of uncontrollability and structural con- 


trollability, as will be discussed in Section 3.10, 


3.4.1 Undisturbability in terms of signal flowgraphs 


State undisturbability of linear systems can also 
be considered in terms of a signal flowgraph formulation. 
Consider the following third order open-loop system with 


the pair (A,d.) in the following form: 


aii ai 0 0 
Artes, oat d, =| 0 (3.10) 
rar Mere | oe 435 


where O35 and all elements of A are nonzero except for 


and a The signal flowgraph of the system charact- 


a whe 


rized, by the. pain. (3.10) 1s shown, in Figure) 3.1. 
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In signal flowgraph terminology, the state vari- 
able, Xin is nonaecesstible to aa if and only if there is 
no possibility of reaching Xs starting Eom and proceed- 
ing in the direction of the arrows along any of the paths 
in the signal flowgraph of the pair (A,d.). Similar con- 
Siderations for the closed-loop system, Hae oas lead to the 


following equivalent definition of state undisturbability. 


Definition 6: State Undisturbability in Stgnal Flowgraph 
Termtnotlogy: 


he ac 


Stare, varvaple., Les of a system character- 
tsed by the pair (HL), ts satd to be undtsturbable wtth 


respect to dtsturbance, E 4 Cais x ts nonaccesstble to 2x 


Fig. 3.1 Signal flowgraph of the pair a) as in 
eqn. (3.10). 
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3.5 Necessary and Sufficient Conditions for Output 
Undisturbability Using an Algebraic Approach. 


The necessary and sufficient condition of Wonham 
and Morse [12] for disturbance localisation using state 
feedback control is based on the concept of an (A,B)- 


invariant subspace. This condition is as follows: 


Dc y* Sieruh) 


where 0D denotes the range space of D and V* is the maximal 
(A,B)-invariant subspace contained in Kernel C. 

It is well known that the (A,B)-invariant subspace 
is spanned by the eigenvectors of the closed-loop system, 
and hence it is possible to express the necessary and 
sufficient conditions for disturbance localisation (or 
decoupling) in terms of the closed-loop system eigenvectors. 

Assuming distinct eigenvalues for the closed-loop 


= (t) of the system of eqn. (3.3), can 


system, the output y 
be represented by 


e 


y(t) =C [W exp{A (t-t)}V]DE(t) dt (3-512) 
fo) 
where W exp{, (t-t) }V = exp{H(t-tT) } 


that is, W is the closed-loop system eigenvector matrix, 


weve 


, and since distinct eigenvalues have been assumed 
A is a diagonal matrix whose diagonal elements are the 


closed-loop system eigenvalues. 


38 


_ bsurrrQed:.t saidep ta as pole os LE 


eftor i 


- 
f . SuisguO 20% _ snot Ried’ anet: dipioi tad bak 
7 


- emast to rs Liuees. Frat 23 e 5 Aue viseasasa 0 - 
/i see J Ta : 7 . 


ce tan iP sf OLE aR, S.A]; We 


s Ba 


noi peel fess 


rm Slt 11 i at ord Ameo dvodbde t 
8,4) me a i Perot 4 oO f 2 a _ ine ans ee 
7 neecue aneiees nuk 


7 we 
| yy 

¢ — » Fi 
C12 2k) : 7 . 


9 Bry y esede sere St Seaenen 0 #te iw 
LiL 48) "4 harm Legioo A = Awe ihn keh aa) 
a ] | a ; 
| ial oe _— WOM 7 - a 


4 te gsovseucgls sie xe tomar 


SQ" mcig “Seeraxs: CU 91.05 RES SS tk sone ‘Bas 
cence) oF 2a0i4\ineo siaenon te 

« — os 
Tete sf ef, BO agset a panto ” wef 


sO? eo! rien rsenicda ls aki : 


my 


SZ 
ger CELE) Voge Fo) meteye 3249 tq 19) “ep Jagan a antag had 


” 
we hesawe bg — 7 
— 
. 


a i 


7 


= 
‘<) & (pi SC(Vatr <1) Aas, Wy | = {2 


Oo 7 
Cir—d) Hgke = Vii tpt} ieee 


| {RISIHN, YOI95 vans fuse qoor tefisgole 4 ati 
7 
‘ 


; OM Aas mAs ape nities torte bh 
= ae - ie = on 


Pe 


=>) 
rm 


wi 


. é 
UY 
~e ay 


iz 
, 
eT 


Hig ate = etnbmate 'f te ae enone 


it es iv 
G" er 


Now define the m x q matrix M(t) by 
M(t) = CW exp{A(t-t) }VD (eg ol) 3) 


It is then clear from eqn. (3.12) that the necessary and 


SUEFIcCIent condition Lor output undisturbability is 
M(t) = 0 LOL Tea SO (3204) 


In terms of element m4 (Tt) of M(t), the condition of 


eqn. (3.14) becomes 


T 
W =’ = 
e pV pdjexPta, (t tt 0 


m5 (T) = 


| ms 


135) 


where the superscript T denotes the transpose of a 


eh th 


vector, ey is the k LOW,oOL Cc, Ws and Vv; are the i 


eigenvector and reciprocal eigenvector, respectively; 


h th 


qd. is the 4° column of D and hs is the i eigenvalue. 


Assuming that the eigenvalues are distinct, and that each 


element of M(t) is equal to zero for all t > 0, it is 
easy to show that the necessary and sufficient condition 


expressed by eqn. (3.15) is equivalent to 


cw vid. =O... ke m, 1 & nyajeas (3.16) 
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abbreviated as k em, ié€n and j € q, respectively. 


Using eqn. (3.16), the following result is stated in the 
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form of a theorem: 


Theorem 3.1 


The econdttton a One = 0, Kk € my 2 © Hy 9 eq. can De 
Ko i ta 

satisfted tf and only tf DCVcKer C, where D denotes 

the range space of D; and V ts the space spanned by at 


Least q ltnearly tndependent etgenvectors. 


Proof of Theorem 3.1 


To prove the theorem, define an nxn real 
matrix G consisting of linearly independent columns 


[gy. Joree-rG] such that 
[gy Jore--/D,] = (4) 1%57%3,---,RE(B,),Im(B,),...] 


in any order, where a, are the real closed-loop system 
eigenvectors and B. are the complex symmetric eigenvectors. 
Re (8; ) and Im(8; ) represent the real and imaginary part 
of Bar respectively. By the original definition of W, we 
represents any one of a78s Or 8. (the complex conjugate 
of B )- 

For necessity suppose that wi F Ker ce dee dhe 
then (cw, )# 0, at least for one value of kK for each i. 
This implies that (vid,) = 9, i en, 4-6 q; thateis vb = 0. 


This is not possible. Therefore some of the eigenvectors, 
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say p, must lie in Ker C. After renumbering these 


a = 0, ° keermyyl£) p, Pletithe 


first p columns of G correspond to the real and imaginary 


eigenvectors such that c 


parts of these p eigenvectors, and define V to be 

the space spanned by these p column vectors, {Gp rFogre-erd te 
Then V © Ker C. Now, for each i ¢ {1,2,...,p} there exists 
at least one k such that clws #0. This implies that 
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columns of G span the real n-space, Se can be written as 


aa 
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From this inclusion relationship it is also easy to see 
that p 2 q, that is V is Spanned by at least q, linearly 
independent eigenvectors. 

With {Gy rForeee Gy} defined as above, and V as 
the Space spanned by these p columns vectors, the condition 
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Remark 3.1: Note that expression (3.17) represents a 
Weaksr Condition than (3.11). The result in (3.17) 1s a 
structural result rather than a synthesis result but in 
conjunction with the eigenvalue/eigenvector assignment 
technique it provides useful insight into the synthesis or 


the design problem. 


3.6 Output Undisturbability Using Feedforward Control 


It is easy to extend the necessary and sufficient 
condition of Wonham and Morse [12] given by expression 
(3.11) to include the case when feedback plus feedforward 
control is allowed. This can be done by simply replacing 
D in expression (3.11) by L, where L denotes the range space 
of Lo=) 0--+ eee (eqn. (3.5).) The necessary and suf- 


ficient condition for disturbance localisation using feed- 


back plus feedforward control then becomes 


pe aal2 oe AND) 
FF 


or, Since L = D + BK , expression (3.18) can also be 


written as 
DCU*r 4 8B (370.9) 


Bhattacharyya [12] has derived just this result by start- 


ing with a transfer function, GAR ics) relating y and 
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—€, and investigating conditions under which G.(s,K,K) = 0 
fOr. alsitts:: 

G, (s,K,K’') = C(sI-(A+BK))"(D+BK"") = 0 
POL alles ae 0) 


However, when feedforward control is used alone 
Bhattacharyya in [12] claims that the solvability condi- 


TOT OT Ls 
Dc B C320) 


Unfortunately this condition is incorrect, as shown by 
the following discussion, (see also Shah [19]). 


When feedforward control only is used the transfer 


FF 


Eunction, G,(s/K ), relating y and € is 


FF BF 


me, many = c(st-a)~!(p+Bx**) (3.22) 


The condition for complete disturbance rejection is that 
yGb)= ONforvall,t°> OPand for every possible €({t). This 


is equivalent to the requirement that 


G.(s:K') £0 for all s (a25) 
Now define the subspace 
n-1 i 
@ = f) (Ker CA’) (3.24) 
i=0 


which is the subspace of unobservable states, that is, 
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(C,A) is observable if and only if 6 = 0. The Cayley- 
Hamilton Theorem insures that 9 is A-invariant (e.g.: 
x € 0 implies A x € 6) (see Appendix 1). In fact 6 is 
the largest A-invariant subspace contained in Ker C. 

| The solvability condition for disturbance rejec- 
tion with only feedforward control allowed, is now pre- 


sented in the form of a theorem. 


Theorem 3.2 


FF 


[Neve errecs, a ae SUC iat G(s,K JOO PO rea? sur 


One only, tf 


DCB +0 


Proof of Theorem 3.2 


Necessity: Bhattacharyya [12] has shown that condition 


(3.23) is equivalent to 


Pe 


CA*(D+BK'") = 0 or CA"L = 0 ie BOS LUO een. 


\eoe92 5) 


This equivalence relationship can be shown by writing 


the characteristic equation of A as: 


n-1 n-2 


|G, (s) | =g -a,s =a Oe Bs -a =0 
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Using a Bezoutian expansion for (sI-A) 2 (Gantmacher [26]) 


it follows that 


G, (s/K’) See tion icp ug! > Cay lal ag Re Sere 
|G, (s) | 
+ CB _jIs + CB. jb} (3.26) 
where 
Bo = a gr a = ea ea aate 5 & 2 juste 


Using (3.26) 41t issnot difficult to show: that condition 


(3,23) 918 equivalent to’ (3.25). From (3525) 


ie oh 


Range (D+BK §)c 6 (3.27) 


Expression (3.27) implies that 
0542.6 58 (3.28) 


which proves necessity. 


Sufficiency: It is possible to chose fea so that 


Range (p+BK  ) es, {6 
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Range (D+BK' *) <= §@ ¢ Ker C 
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hence 


EF; 


CA (Deer MS Oe 4-92 001) ane | 


which is equivalent to 


FF 


H(S,K. Joes OF “for all <s 


This completes proof of the theorem. 


Remark 3.2: The solvability condition expressed by 
(3.28)is more general than the solvability condition (3.21) 


of Bhattacharyya [13] as shown by the following discussion. 
FF 


In the conventional feedforward control problem K is 
determined by solving the equation: 
DiABEKa Gu (ie 29) 


Note that for completely observable linear dynamic 

systems, condition (3.28) reduces to condition (3.21) as 
presented by Bhattacharyya, and both are equivalent to 

eqn. (3.29). However, for the case where 6 is non-zero, 
condition..(3.21)-£ailsebut (8.28) still holds’ and in fact 
shows that the solvability condition (3.28) is more gen- 
eral. As is obvious from (3.28), the solvability condition 
for unobservable systems is weaker than that for observ- 


able systems. 
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In summarising the main contents of sections 3.5 
and 3.6o0ne can say that the necessary and sufficient 
conditions for the existence of state feedback controllers 
that will decouple the effect of disturbances, &(t), 
from outputs, y(t), can be defined in geometric terms as 
follows: 

1) the range space of D must lie within the maximal 


(A, B) - invariant subspace contained in Kernel C 


(Wonham and Morse [12], Bhattacharyya [13]) or that 


2) the range space of D must lie within the space 
spanned by at least q of the closed-loop system 
eigenvectors (where q is the dimension of the 
disturbance vector, &), which in turn must be 


contained in Kernel C (Cf. Theorem 3.1). 


However, in the results to be presented below these 
conditions are not tested directly. Instead it is shown 
via Theorems 3.3 and 3.5 that the necessary and suffic- 
ient conditions for state undisturbability can be expres- 
sed in terms of the structure of the system matrix, H, or 
equivalently, the structure of the closed-loop eigen- 
vector matrix, W. Then a design procedure is presented 
that generates the required controller matrices based on 
eigenvalues and elements of the closed-loop eigenvectors 


chosen by the designer. 
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Parnes tunbabi lity sand Closed-loop System Structure 


The following result forms the basis for design 


of undisturbable systems. 


Theorem 3.3 


A system*characterised by the pair (H,2.) has k(1<k<n) 


undtsturbable state vartables with respect to mae i 


and onty 27, 


1) mR ORK (L.) = l.and the patr (H,£.) Loti thes fOL— 


Lowutng form, ov, can be. brought inte thre 7 orm, by 


@ suttable reordering of the state vartables: 


Ee ! 0 0 
y 3 
iis al 3 
H H 
sl 4 a) 
where H. € Rake H, & eats eee ke 
i 4 es) 


(sien 0) 


and Lys Los +++s Hy are the k undtsturbable state 


Vartaples OF Che pac (Hots). 


Ox 


ti) 02 = 90. “There eorresponds oO ad trivial case. 


To prove the theorem the following lemma will be used, 


Herstein [27]. 


* Note: Hye H3, Hy 


only restriction is H. = £(A,,B, Ko) = 0. 


and oe can be time varying since the 
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Lemma 3.1 


Let M be a nonempty set of nxn matrices in Rr”, 
if W(WER") is an invariant subspace under M e M and is of 
dimension (n-k) over R, then there exists a basis of ee 


over R such that every M e M in this basis is of the form, 


(331) 


where M, € pEXK ona Mae p (n-k) x (n-k) | 


this basis M, is the matrix representation of the restric- 


ELON Of M to. Wand M, is the restriction of M to RV /W. 


The proof of the lemma is straightforward and is given 


Furthermore, in 


in Appendix A. 


Proof of Theorem 3.3 


Let the first k state variables of the system 
characterised by the pair (H,£.) be represented by xl, 


Pe. 
Tr 
Sel) ae [x, (t) Xo (t) eee x, (t) J 


Sufficiency is clear from the structure of the pair 


a any feqri.) 3.30) eerscunce 


te) ="Q0"rGreall t > 0 and arbitrary au las 


7 ] 
lan 
“< aT Boe 3 Motebae aes Jnl SETS 
| 7 

Me "vo eiesd & BteiKe sisht neds .A I9Vve 


- 
,arnot uae] at 21 aleoe aide iy M y WV VY Tava ten 


; _ 
ty 
(£.@) fea 
4 
if 
45 
fm i i x 
i 1 Sais j 2 1 hae 
7 : 
: | 2 7 ' ss it 
@ +4 if Sti Sa €04 7 T) —s ti ,™ oiend phy 
2 2 ie oe 
rh hy ao anism arer eat 2 it ore Y oF wi “mort 
. “22 2) hme a Ro Boog amt 
tar typ | Oris ld rtoe er dh > Bal 4 “ pares Lf - ed! a 
. 
A BEE 
ta x % J 
_— —a 
mare yt &;! ft, f ai rié 4 = 


Lk pix yeh beoeastged ed (23,8) tind ond ¢ beet; came 


. a : 
: a 
. a8 


: v) (poe . (e) + 43) 5) = yin 7 

7 a ad 2 . ; 
a>) 1m | . Gio : 
eb Sc adi Yo sxustoucie wid. more oa bd, ae 


as . 


Necessity Consider a system characterised by the pair 
ee) which has k state variables which are undisturb- 
able with respect to Bae Let the system be brought into 
the form of eqn. (3.32), by reordering the state variables 
if necessary, such that x = [xl 52 1 where xl = [x1 Xos 


ee ie contains the k undisturbable state variables and 


kxk (nai) caine) k. 
Hy e R F Hy, e R and Hoe Gu es 
iS & 
ee (Ae) Hy pola x1?) 
no eho = ---f----- ral la ct sae 
x2© (t) Pettit x2© (t) 
ei) i 4 
(SS) 
Since undisturbability implies x1° (t) =— OVPOrearoLeLary 


B, (t) end, tor aii, « > O,con. (3.32)§ reduces to the <fol— 
lowing two equations: 
E e 

0 = H,x2 Gee ne a ee (3,233) 


Mee = 3 
22 (b= H,x2 (t) + £555, (*) (3.34) 


Let W be the H-cyclic subspace generated by Lar as defined 


below. 
ee Ss Se er ey (3.35) 
j j J 


In fact W is the H-invariant subspace of states which = 


Can control. Let dim(W) = 9p. Then it is clear that 
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p < (n-k). Since the effect of 5 on x is given by egn. 
(3.34), it is correct to say that the matrix represent- 

ation of the restriction of H to W in the standard basis 
is H,- Two cases will now be considered: (i) dim (W) = 


pe="n-k and?” (17) "aim (W)"="p-<(n-k) < 


Case l: p= n-k 
Since the matrix representation of the restriction of 


H to W in the standard basis is H then by Lemma 1 the 


jal lé 


p standard basis vectors: { Bets e,} span the 


“k+l? §K42" 
H-invariant subspace W. This implies that ae = 0, and 
Since x2 is disturbable, it follows from eqn. (3.33) 


that H = 0. 


Case 2: p< (n-k) 
Now Hy can no longer be the matrix representation of the 


restrIiceron Of H-to Wy but 2e Stil. follows ™thacep ort 


the (n-k) standard basis vectors, {ep 444 Cpagr cece ents 
span W. This once again implies that Ah = 0, and as in 
the previous case, that H. = 0. 

In either case ue = 0 and H. = 0 which means 


that the undisturbable system (3.32) exists in, or can be 
brought into, the form of eqn. (3.30). This shows neces- 


sity and completes the proof of the theorem. 
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Remark 3.3: The condition rank ) = 1 is equivalent 
to the condition k < n in the sense that each implies 

the other. The significance of the condition, rank ae = 
l(or k < n), is that a system characterised by the pair 
(H,£.,) has at least one disturbable state variable. This 


conclusion also follows directly from geometric concepts 


and the principle of conservation of energy. 


Remark 3.4: A non-trivial system (i.e. oe # 0) charac- 
terised by the pair (H, 2.) has no undisturbable state 
variables (i.e. it is completely state disturbable) if 


either se or H do not contain any zero elements. 


Theorem 3.4 
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Proof of Theorem 3.4 
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sl sceitb) po eter or ea fee 


Using eqn. (3.26) the above condition can be shown to be 


equivalent to 
m 
c.H toe = 0 mn ce Op ipiciene ple (33.7) 


Both sufficiency and necessity then follow easily from 


eqn. (3.37) and eqn. (3.35) which defines W. 


Remark 3.5 Note that Theorem 3.3 can also be derived 
as a direct consiquence of Theorem 3.4, in the following 


Manner. For system (3.32), where x are 


i Or eooo ge — 
the undisturbable state variables, by Theorem 3.4 the 
following necessary and sufficient condition needs to be 


satisfied: 


Aly 
W Cc Ker [e, re e,] Geo8) 


Die ke 


where W is defined by eqn. (3.35) and Cyr Cor seer Sy are 


the standard basis vectors. This clearly means that 


fe. aye Cpygrocees ora span W. Hence it follows that 
sale =O, and from edie ais eo, Chat H, = 0. 
Eqn. (3.38) also provides an interpretation of V* 


for the case when Vie ars 2 poe ta ke ehOrMenats. cose, 


V* is the largest subspace of states which can be made 


unobservable at y = Cx, where m= k. 
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3.8 Closed-loop Eigenvector Structure and Undisturbability 


The following theorem relates the property of 


undisturbability to the closed-loop system eigenvectors. 


mmMeorem 34.5; 


Constder a system charactertsed by the pair (H,k) where 
H has dtstinet etgenvales; this system has k(1<k<n) un- 
disturbable state vartables wtth respect to disturbance, 


5 3 ees Onl COMETS 


t) the elosed-loop system etgenvector matrix, W, 
ts of the followtng form or can be brought 
tnto the following form through rearrange- 


ment OF Tus CoLunne: 


W4 ! 0 
W = |-------- (3.39) 
1 
i 
oe 
where W, € geek W, ein k)atn-k) and 
wt) ee LS Of acne form: 
0 
ete fe. (3.40) 
J 
iC 
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where Le. = pin-k) 
J 


and the first k elements 
are gero and correspond to the k undisturb- 


ablevetate variaplec. 
To prove the theorem the following Lemma will be 
used. 
Lemma 3.2 


Let M be a set of nxn matrices in R”™ with distinct 


eigenvalues. Then M é M is of the following form: 


{ 
M, | 0 
M = }---4---- (3.41) 
{ 
| 
ee 
where M, ¢€ raise isp M, € Rae pope bean Ro le a Salvi eae 


and only if W, the eigenvector matrix of M, is in the form 


Shown in eqn. (3.39). 


Proof of Lemma 3.2 


With W invthe ormof eqn. (3.39), suLficiency is 


1 and A denotes the nxn 


clear since M = WAV, where V = W. 
diagonal matrix whose diagonal elements are the eigenvalues 
OLaM. 


To prove necessity, note that since M in eqn. (3.41) 
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is a lower quasi-triangular matrix, its determinant is 
the product of the determinants of the diagonal blocks, 
This means that, 


Gantmacher [26], i.e. |M| = |M,|] |m 


Pele 


o(M) = o(M,) U o (M4) 
where o(M) denotes the spectrum of M. For this part of 


the proof assume that A and W are partitioned so as to 


be compatible with the partitions of M in eqn. (3.41): 


1 
en i vated 
A = ----+--- W = ses (e354 '3)) 
i i 
i 
0 | Ay Ws | Wy 
where Ay E cs No E g rk) a {n-k) and they are arranged 
Such that: 
0 (M, ) = diagonal elements of Ay (3.44) 
o(M,) = diagonal elements of A, (3.45) 


Since W is the eigenvector matrix of M, it follows that: 


MW = WA (3.46) 


Bromequs. (o.4.), (3.43). and (3.46), it foliows chat: 


M,W, = WA (3)04.7) 
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and 


MWe t= "WA (3.48) 


Clearly, Wy in eqn. (3.47) is the eigenvector matrix for 


M Since Ay # Nor eqn. (3.48) implies that W, = 0, 


i. 2 


which proves necessity. 
Proof of Theorem 3.5 


In view of the necessary and sufficient conditions 
for undisturbability in Theorem 3.3, it suffices to prove 


that the closed-loop system matrix H is of the form 


5 (3.49) 


am 
II 
i 
I 
I 
! 
—----p---- 


pkxk k) x(n-k) 


where H, € ik and OnLy win its 


uM 


eigenvector matrix W is, or by rearrangement of its columns, 


and H, € p (n= 


can be brought into the form shown in eqn. (3.39). Clear- 
ly, the proof for sufficiency and necessity both follow 


£Erom the result of Lemma 3.2. 


Remark 3.6 The assumption that H has distinct eigen- 
values entails little loss of generality, since almost 
any feedback matrix, K, will yield a closed-loop system 


with distinct eigenvalues (Davison and Wang [28]. 
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Remark 93.7: The extension of the results of Theorem 
3.5 to discrete-time systems is obvious. 

Theorem 3.5 helps to relate the two sets of 
necessary and sufficient conditions considered in 
Theorems 3.1 and 3.3 in the following way. The necessary 
and sufficient conditions for disturbance localisation 
in Theorem 3.1 are expressed in terms of the subspace 
Spanned by some of the closed-loop eigenvectors, i.e. 
the (A,B)-invariant subspace V. The necessary and suf- 
ficient conditions for undisturbability in Theorem 3.3 
are given in terms of the structure of coefficient mat- 
rices of the state-space model. Theorem 3.5 is basically 
an extension of Theorem 3.3 and shows that the two sets 
of necessary and sufficient conditions in Theorems 3.1 
and 3.3 are equivalent. In contrast to previous results 
On disturbance localisation (Wonham and Morse [12], 
Bhattacharyya [13]),the results in Theorems 3.3 and 3.5 
are structural results. However, in conjunction with 
the eigenvalue/eigenvector assignment technique these 
results lead to a very simple design procedure.” In 
contrast to that the results of Wonham and Morse [12] do 
not lend themselves easily to computer implementation. 
One of the main difficulties is the evaluation of V/* 
which is subject to the hazards of ill-conditioning, 


Bhattacharyya et al. [23]. 


* The design procedure based on Theorem 3.5 allows the 
user to specify the elements of Wie which is not possible 
in previous methods [12]. 
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3.9 Output Undisturbability and Closed-loop System 


Structure 


When the system outputs are a subset of the 
state variables then the condtions for state undisturb- 
ability (Theorems 3.3 and 3.5) can be applied directly 
to determine output undisturbability. However, when 
outputs are a linear combination of the state variables 
then conditions for state undisturbability do not apply 
directly. In this case output undisturbability can be 
handled by transforming the system prior to applying 
the conditions of Theorem 3.3. The first step is to 


define a new vector, z, of dimension n, such that (cf.[36]) 


where C is the mxn output matrix of the original system 

(assumed to be of full rank) and E is an (n-m)xn matrix 

chosen arbitrarily such that its (n-m) rows are linearly 
independent of the rows of C. Under these conditions Q 

exists, and the original system defined by eqn. (3.1) 


becomes: 
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The closed-loop systems can now be represented by egn. 


(333) ewhere®: 


L 


1, O40 °-B kK and (3.52) 
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The conditions for state undisturbability can therefore 


be applied directly to the transformed system and since 


{Zz = y,, i=l,...,m} they provide a basis for the design 


of systems with undisturbable outputs. This is illustrated 


in one of the following examples (Section 3.12). 


3.10 Discussion of Theoretical Results 


3.10.1 The relationship between undisturbability and 
uncontrollability 


From a geometric viewpoint the essential differ- 
ence between undisturbability and uncontrollability is 
the following. With W defined by (3.35), then from 
Lemma 3.1 it is clear that an open-loop system character- 
ised by the pair Wetha is “uncontrollable" if and only 
if a linear transformation can bring it into the form of 
(3.30). On the other hand, the pair i has k undis- 
turbable state variables if and only if it can be brought 
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of linear transformations, namely, permutations of coor- 
dinates to reorder the state variables. This class of 
linear transformations also allows the undisturbable 
State variables to be identified. These considerations 
lead to the following conclusions: 

i) complete state controllability of the pair (Ad) 


implies complete state disturbability of the 


Pair eavastey) However, the converse is not always 
true. 
ii) if some of the state variables of the system 


(A,d.) are undisturbable, then it follows that 

the pair (OnE) is uncontrollable. However, the 

converse is not necessarily true. Even if the 

pair ae is uncontrollable it is possible that 
ge) is completely state disturbable. Note 

that analogous results also apply for pinout 

disturbability. 

As is evident from the above discussion, undis- 
turbability of the pair Te) cannot always be determined 
from controllability concepts. However, uncontrollability 
of the pair (Ads) can often be inferred by application 
of the undisturbability criterion; that is, the pair 
(eve is uncontrollable if the pair (A,d.) has Some 
undisturbable state variables (cf. Section 3.10.2). These 
relationships are graphically summarised in the venn dia- 


Gram Of Figures 372 and’3.3. 
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Venn diagram to graphically illustrate 
the relationship between disturbable, 


controllable and structurally control- 
lable systems. 
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Venn diagram to graphically illustrate 
the relationship between undisturbable, 


uncontrollable and structurally uncontrol- 
lable systems. 
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Jwl0L2) Thesmelationshipebetween undisturbabi la tyvand 
structural» controlilabil ety 


The class of undisturbable systems are a special 
class of uncontrollable systems, namely those which are 
not "structurally controllable", a concept recently in- 
troducedubys lin 120ml Structural controllability: isa 
property that is useful in determining the controllability 
characteristics of a system from a knowledge of only the 
zero entries in matrices of the pair (A,B). A system 
which is not structurally controllable can be brought 
into two possible forms through permutation of coordinates. 
One of these forms is identical to the form of (3.30). 
Thus there is a close relationship between undisturbabil- 
ity and structural uncontrollability. Undisturbability 
can also be investigated without knowledge of the exact 
numerical values of the non-zero elements of matrices A 
and D since only knowledge of the zero elements of these 
matrices is required. 

Theorem 3.3, which is based on algebraic concepts, 
can be used to show the equivalence of the following two 
properties for the system in eqn. (3.1). 

i) The signal flowgraph of the nontrivial pair (A,d.) 
does not contain any non-accessible state variables, 


(i.e. the system is completely disturbable). 
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ii) there is no permutation of coordinates which 


brings the pair (Ard) into thertorm of (3.30). 


By contrast, Lin [29] has used a graph theoretic approach 
to study structural controllability and to show the 
equivalence of results similar to (i) and (ii) above 

for structural controllability. The relationship be- 
tween disturbability, controllability and structural 
controllability is summarised in the venn diagram of 


Riguresecy2.and 73. se 


oe 107.29 Theerelationshipebetween undisturbability and 
zeros of Rosenbrock's system matrix. 


In Section 3.10.1 the relationship between 
undisturbability and uncontrollability was considered. 
Earlier we also remarked about the relationship between 
undisturbability and unobservability. In the ensuing 
discussion we will investigate a more definitive mathe- 
matical relationship between these three concepts by 
considering Rosenbrock's systems matrix. In Rosenbrock 
[35] the ideas of controllability and observability have 


been embodied into the following definitions. 


Controrlaprtirty: “The pair (A,B) 1S controllable 1f and 


only if (sI-A) and B are relatively left prime. 
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Observability: The pair (C,A) is observable if and only 


if (sI-A) and C are relatively right prime. 


According to Rosenbrock [35], the eigenvalue 
corresponding to a mode which is completely decoupled 
from an input is said to be an input decoupling (i.d.) 
zero. In the same manner, the eigenvalue corresponding 
to a mode which is completely decoupled from an output 
is Said to be an output decoupling (o.d.) zero. Finally, 
the eigenvalue corresponding to a mode which is completely 
decoupled both from the input as well as the output is 
Said to.sbeyam input @eutput decoupling (1.0.d.) zero. 

With this brief review of basic system matrix terminology 
we now present the following result in the form of a 


theorem: 


Theorem 3.6 


An ee order system characterised by the trtple (H,L,C) 


te Completely output, unateturbable tf and. only. 27, 


where 25 denotes the number of t.d. zeros, Bo denotes 
the number of o.d. zeros, and 255 4s the numoer Of 2.0.0. 
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Remark 3.8: The statement of Theorem 3.6 in effect 
says that (sI-H) and L are relatively (left) prime and 
(sI-H) and C are relatively (right) prime. Note also 
that the triple (A,D,C) can be regarded as a special 
case of (H,L,C) and hence the result of Theorem 3.6 


applies to both open as well as closed-loop systems. 


Proof of Theorem 3.6: 


SULTICIeNCVE, ls Cleareornce Ag Z ae 2 an implies 


that the minimal system order = 0, that is to say that 
Loe 


vis) = e(Si—H) “L=)0. 
Necessity: For this assume that the outputs are a sub- 


set of the state variables. If this is not the case, 
the system can always be transformed into a set of new 
coordinates where the outputs are a subset of the new 
State variables (cf. Section 3.9). Once this is done 
then using the arguments of Theorems 3.3 and 3.4 it can 
be shown that the system matrix P(s) for the undisturb- 
able system characterised by the triple (H,L,C) is of 


the form: 
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With the above form of P(s), the transfer function be- 


tween y(s) and €(s) can be evaluated as: 


aS) ie C(st-H) +n, or 
= -1 
i 0 sI-H, 0 I O;; 1 0 
y(s) = [I 0] x 
0 SI-H, 0 I -H, I/} 0 sI-A, 


Wl 
i=) 


SI-H, 0 0 een er 
x = [I 0] 
0 AL, Lo “Hy AL Lo 


From the above matrix equality it is clear that the i.d. 
zeros are the poles of |sI-H, |, i.e. the eigenvalues of 


Hyyeand=the=o,d = zerocearerthe poles of |sI-Hj)| Or the 


4 


eigenvalues of Hy. Now H is a quasi-triangular matrix 


therefore eigenvalues of H are o(H, )Uo(H,). Since H 


is nt order, Zs + Zy = n- Note that in some cases the 


iseteor ind. zeros itseteot o.d. zeros) “+ 0 tergwat Lo 
has one or more zero rows), in this case the system has 
(i1.0.d.) zeros as well. Because of this Z.+Z, > n and 


hence we need to subtract Ze to satisfy the equality: 


O 
Zot 2a 2k = Noe SioecOnpLe ces. CNS /proot sot ithe 


theorem. 
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Remark 3.9: It is interesting to note that conceptually 
the proof of this theorem becomes readily apparent when 
the definition of undisturbability is considered from a 
Signal flowgraph terminology - just as the concept of 
decoupling zeros becomes clear from a signal flowgraph 

of a system in the Kalman canonical form. For example 

the eigenvalues associated with the non-accessible modes 
or the undisturbable modes can be identified with i.d. 
zeros and the eigenvalues associated with the accessible 
modes (the disturbable modes) can be identified with the 


O.d. zeros since there is no path from € to y. 


3.11 Design Considerations 


It was mentioned in Section 3.5 that the necessary 
and sufficient conditions for disturbance localisation 
reported by Wonham and Morse [12] and Bhattacharyya [13] 
were expressed in terms of V*, the maximal (A,B)-invariant 
subspace contained in Kernel C. Athough V* plays a key 
role in a number of control problems (Morse and Wonham, 
[30]) it is not easily calculated (Bhattacharyya [31], 
Bhattacharyya et al. [23]). By contrast, the result of 
Theorem 3.5 in conjunction with eigenvalue and eigen- 
vector assignment techniques allows insight into the 


design problem and is also useful from a synthesis point 
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In this section a design procedure to achieve 
disturbance localisation will be developed based on 
Theorem 3.5 and the eigenvalue/eigenvector assignment 
technique (Srinathkumar and Rhoten [32] (cf. Chapter 2)). 
From, Theorem 3.5 it is clear that in order to achieve 
disturbance localisation, K must be specified so that 
the closed-loop eigenvectors are in the form of eqn. (3.39) 
and if feedforward control is necessary, then poe must 
be specified so that L is in the form of eqn. (3.40). 

The manner in which the design objectives of arbitrary 
eigenvalue assignment and disturbance localisation in- 
fluence the design procedure will now be considered. 

The closed-loop system matrix can be expressed as: 

H = A + BK = WAV (37.93) 
where V is defined by, V = Wine Bon 2.( 35.53) mcanabe 


partitioned as: 
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In equation (3.54) the partitions of H and W are denoted 
by double subscripts (e.g. Wi4) to avoid confusion with 
the partitions in eqns. (3.30) and (3.29) which have 


different dimensions. Matrices H A 


Tale nee ek ae Se 
and Vil are rxr matrices where r is the number of inputs; 


Hoos A A and V are (n-r)x(n-r) matrices and 


os Minot fiot 22 


By is assumed to be nonsingular. (If necessary the 


state variables can be reordered to make B, nonsingular. 


il 
Note that W and consequently A and V can be complex 
matrices.) From eqn.(3.54) it follows that Range (B,) = 


ee and thus a K exists such that the rxr elements of W 


tae 
and the rx(n-r) elements of Wi can be chosen arbitrarily, 
subject only to the requirement that W be nonsingular 
(since VY = wot). Once Ay and No (i.e. all the closed- 
loop system eigenvalues) and Wii and Wi> are chosen, Woy 
and Woo can be calculated (Srinathkumar and Rhoten, [32]) 


from eqns. (3.55) and (3.56) which follow from eqn. (3.54): 


Wo Ay ~ PWo4 = TWiy + SW,44) (3255) 


(3756) 
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Then K can be computed from eqn. (3.53) by using the 
pseudoinverse of B (Chapter 2, see also Shah et al. [33]). 
Note that the pseudoinverse solution for K provides an 
exacc solutwon: to eqn. (3.53) since Woy and Wo have 
been chosen so as to satisfy eqns. (3.55) and (3.56). 
Specifying the first r elements of each of the 
n eigenvectors requires only (rxn-n) degrees of freedom, 
i.e. (rxn-n) elements of K, providing that at least one 
of the r specified values in each eigenvector is nonzero. 
The remaining n degrees of freedom can be used to arbit- 
rarily assign the closed-loop eigenvalues. Thus in 
general, simultaneous eigenvalue and eigenvector assign- 
ment are possible provided that at least one of the 
first r entries in each eigenvector is chosen to be 
nonzero, and that the final Mcnveatan matrix W is non- 


Singular. 


S-il.1 Steps inethesdesign yorocedure 


The main steps in the design procedure can be 


summarized as follows: 
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After formulating the control objectives, matrix 
D is inspected to determine if the selected state 
or Output variables can be made undisturbable by 
feedback control or whether feedback-feedforward 
control is required. If D has the form in eqn. 
(3.40), then only feedback control is required; 
if not, then feedforward-feedback control must 


be used. 


Matrices Wii and Wi> in eqn. (3.54) are assigned 


according to the result in Theorem 3.5. 


The closed-loop eigenvalues in Ay and Ao are 


assigned. 


Matrices Woy and Wo5 are calculated from eqns. 
(3.55) and (3.56) and a check is made to determine 
whether the resulting eigenvector matrix, W is 


nonsingular. If W is singular, step ii) is 


repeated with a different choice of Wii and Wio- 


Feedback matrix K is calculated from eqn. (3.53) 
uSing the pseudoinverse of B. 


If feedforward control is used, ee is then de- 


Signed so that the appropriate elements in L are 


zero. (cf. Section 3.12 for illustrative examples). 
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Note that the present design procedure suffers from the 
disadvantage that if at is singular, then not all of the 
control objectives can be met and therefore new control 
objectives have to be formulated such that the final By 
obtained after reordering state variables is nonsingular. 
However, there are exceptions to this case; that is in 
some instances the original design objectives can be 

met even if By is singular, as for example, the distilla- 
tion column example considered in Chapter 4. Moore [34] 
has recently proposed a method for eigenvalue/eigenvector 


assignment which does not require B, to be nonsingular 


i 
and thus provides an important alternative to the approach 


of Srinathkumar and Rhoten [32] as considered in Chapter 2. 


3.11.2 Simultaneous disturbance localisation and egal. 
value assignment 


From a design point of view, it is desirable to 
be able to arbitrarily assign all n closed-loop eigenvalues 
as well as achieve disturbance localisation. However, 
it is not always possible to satisfy both objectives if 
the eigenvector/eigenvalue design technique presented in 
Chapter 2 is used (see also Srinathkumar and Rhoten [32] 
and Shah et al. [33]). In particular, when the desired 


number of undisturbable state variables is less than the 


rs 
b osia ‘tq iin Fon roth ,twlveace oi, 


i 7 : 
fotsnoD wu eictetedd! bap’ tem ac 285 4 
© font? sis J3edat toe bateiodne? aA ce asvitoatde 
f 


a 


be 
tefupniecon 2! nefdairey state po! Tebiess seat st 18 sich ] 


res 
om, we) : 
‘ete (eas el of edots ge=re ate ais revewort 
os) Vo 
ac ns thee Boe Ap Leet , oie snbuases amas 
of) .afenees sot Be ,telogn ta ee ga M4 mere tem 


| «.f qptmce denlor ioid 


ihe} oom t sfasr md 
nerbo\eu layne TT tan shoe » beara qideene ape 
Je hyortt axe 4 mpe oy do eelw, soomere keae 
erage sn ie ode des 2eqsi. Ge Rabivaay wae bas 


4) Jsnoc 2b (Sf) asso bne Sgemumity naka ‘to 


hon ih 2 

7 A: 
~1i3j bos “oabteatiegds sat 

. - 
ot elde+thaei-at) Pie weey Se & mort. 


roulévaonte coaal-Beeoloa m Téa AvLERS 4 : inven bebrs _ 
: <a »~ a at 
terevion ~,doryserhsool sorecss 'h Sveu 


1} esvitoshdo geot yRatgen of ol dreaoy ByaweA 


, est a 4 idenakall 


id HSTipestq 3 pi gitses Ao leeb ori liv ie 
foe, nevons sah Sone: pane sofa ‘ac wee at 
eri Al 4 Tr 


7 : sy 


Bevel Zeb SAg nit! 4 
ran - 
aah 


PL: < 
oo. te 
> 


“iat mee "neny a2et [daiae 


number of control variables (i.e. k<r), then simultaneous 
eigenvalue assignment and disturbance localisation is 
possible, but when k=r not all eigenvalues can be arbi- 
trarily assigned. The reason for this limitation will 


be considered in the following discussion. 


In order to make k state variables Xproceen Xe 
undisturbable, the kx(n-k) matrix, Woe must satisfy the 
condition in Theorem 3.5 that W, = 0. However, in the 


eigenvalue/eigenvector design algorithm, arbitrary 
eigenvalue assignment is possible only if there is at 
least one nonzero element among the first r elements of 
each closed-loop eigenvector. These two conditions imply 
that if k<r, simultaneous disturbance localisation and 
arbitrary eigenvalue assignment can be achieved by 
setting Wo = 0 and assigning nonzero values to at least 
one of the remaining r-k elements in each column of Wid 
(c£. eqn. (3.54)). By contrast, when k=r some closed-loop 
eigenvalues cannot be arbitrarily assigned. In this case 
Theorem 3.5 requires that W, = W = 0 for disturbance 


2 12 


localisation but when W = 0, it follows from eqn. (3.56) 


Ae 


that the n-r eigenvalues of A, are the eigenvalues of P 


and thus they cannot be arbitrarily assigned. Consequently, 


when k=r the stability of the closed-loop system depends 
on the eigenvalues of P. Necessary and sufficient con- 


ditions for simultaneous disturbance localisation and 
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closed-loop stability have recently been reported by 


Fabian and Wonham [17] and Chang and Rhodes [18]. 


geb2eebblustrative Examples 


To illustrate the concepts and design procedures 


of Section 3.11 two examples will be considered. 


Example Bia ics 


A hypothetical third order system with 2 inputs, 


2 outputs and 2 disturbances is shown in eqns. (3.58) and 


G3 6590: 
S095 Gays) 0.75) ie aa 
A =| 1.0 eS =i B =il.0e- -=1 (3.58) 
T.0 =i te25 ii 1 
and 
: ree OF OF Pl 
‘chaps pee tA (3.59) 
Or Soames: 0 2 


The eigenvalues of A are -0.5, -1.25 and -2.25. Theorem 
3.3 and inspection of the system in eqns. (3.58) and 
(3.59) reveal that all of the state and output variables 
are disturbable with respect to each disturbance, since 


A contains no zero elements. 
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Suppose that the primary control objective is to 
make xX undisturbable with respect to disturbance, fic 


For convenience, let hy and Wig denote the eeaies 


y 
elements of H and W respectively. Then k=l and from 


Theorem 3.5 1t. follows that W, = [wy 5 


= 0. This choice means that hi> = hi3 = 


Wy 3! must be chosen 


So that Wi 2 = Wi 3 


0 and thus H and d the first column of D, are in the 


dle 
required form of eqn. (3.30). The eigenvalue/eigenvector 
assignment method of Section 3.11 can be used to calculate 


the remaining elements of W. Since in this example, r = 2, 


all six elements im the first two rows of W (22e. W and 


imi 
Wi> in eqn. (3.54) can be arbitrarily assigned. For 
example, Wi and Wy > can be specified as: 
il 0 0 
Wii = ‘ Wid = (3.60) 
1 i il 
Since Wid # 0 all three closed-loop eigenvalues can be 


assigned arbitrarily. If the closed-loop eigenvalues are 


selected as -3, -4 and -5, then 


il ’ A, = [-3] beso) 
0 =) 


With Wii “2° A and No specified, Woy and W,, can be 


evaluated from eqns. (3.55) and (3.56) and K calculated 


from eqn. (3.54) to give: 
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The closed-loop response to a unit step change in 
bi is™’shown in “Fig: 3342 and illustrates. that xy has been 
made undisturbable with respect to Ei: It is interesting 
to note that X3 is only slightly affected by & - This is 
due to the presence of a small (but nonzero) element in 
hy5- 


Example 3.2: 


Again consider the system in eqns. (3.58) and 
(3.59) but now assume that the control objective is to 
make x1 and X3 undisturbable with respect to both Ey and 
Go- Clearly, this is a much more ambitious design objec- 
tive than the one considered in Example 3.1. 

In view of the structure of the matrices in 
theorems 3.3) and 3.55 at 2S Convenient to reorder the 
state vector so that state variables, xy and Xo which 


are to be made undisturbable appear as the first two state 
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variables. Thus the new state vector, x, is defined by: 


er (3.63) 


i 
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hk 
x 
NO 
x 
1S) 
T 
* 
ke 
x 
Go 
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Similarly, let A and D denote the appropriate permutations 
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of matrices A and D corresponding to the state vector x, 
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and let H denote the corresponding permutation of H. 


~ ~ 


Theorem 3.3 indicates that in order for xy and Xo to be 


undisturbable with respect to aT and Eo: it is necessary 


~ 


ea ho = 0 and the first two rows in L must also 


be zero. Since ia = do4 


K so that x4 and Xo are undisturbable with respect to 1° 


that, h 


= 0, it is possible to design 


~ ~ 


are nonzero, Xy and Xo cannot 
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However, Since d and d 


L2 22 
be made undisturbable with respect to E> by using only 
state feedback control. In this case feedforward control 
of Eo must be used in conjunction with feedback control. 


Consider the feedback-feedforward control law of 


eqn. (3.65)!: 


BE 


u(t) =(Kx (Ce) 4K Eo (t) (36) 
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Combining eqns. (3.65) and (3.1) gives the closed-loop system 


x = H*(t) + LE(t) (3.66) 


The corresponding closed-loop system expressed in terms 
(BS er eae eSB 


~~~ 


x = Hx(t) + LE(t) C367) 
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~ nw ~~ ~ 


where H = A + BK, L = Be + D. 


In order for *, and x, to be undisturbable with 


nf 2 
respect to Sor the second column of itp Los must have the 
form eg = [0,0,a] where a is to be determined. Then ge 


and a can be calculated from 


FF ae 
B K + d., =: L, (3.68) 


to give 
FF ii 
K || | , and ® = 3.0 (3.69) 


To complete the controller design, the state 


feedback matrix K must be calculated. The undisturbabil- 
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ity objectives and Theorem 3.5 require that W. must be 


~ ~ 


chosen as Wo = Wid = 0 while Wid can be selected arbit- 


~ u 0 
Ww = (3% 70) 
J i 1 1 


Since Wi> = 0 and k=r, it is not possible to arbitrarily 


rarily as 


assign all of the closed-loop eigenvalues, as discussed 
ine oectvon 3.lls in fact can. (355/) indicates * that P 

is a scalar quantity with a numerical value of -4.25 and 
consequently, one of the closed-loop eigenvalues, Aor has 
a value of -4.25. Then from eqn. (3.56) if follows that 
Ww = 1. After specifying the remaining two closed-loop 


values to be the desired values of hy = -4 and do = -3, 
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W and K can be calculated using the eigenvector/eigen- 


value assignment procedure of Section 3.11. This gives 
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Figures 3.5 and 3.6 compare the open-loop and 
closed-loop responses to unit step changes in ey and bor 
respectively, and show that xy and X3 have been made 


undisturbable. 


Example 3.3: 


Another 4th order hypothetical state-space system 
will now be considered to illustrate the design procedure 
for achieving output undisturbability. The form of the 
Open-loop system is as defined in eqn. (3.1) and the 


coefficient matrices are: 


er one 0.850 ules yl site} 0.854 


=) 060 Oe 2.060 =1. 030 


ji = 
eo G 0.958 -1.810 0. 7-20 
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The eigenvalues of A are -3.0, +0.1, -1.0 and -0.51 and 
hence the system is open-loop unstable. Also, since A 
does not contain any zeros the pair (A,D) does not satisfy 
the conditions for undisturbability defined by expression 
(3.30) or alternatively by Remark 3.4, and hence all four 
of its state variables are disturbed by either of the two 


disturbance variables. Since Y5=x the second output 


3 
variable is disturbable but since Vy is a linear combina- 
tion of xy and Xo it is not possible to draw any immediate 
conclusions as to whether Yi is disturbable or not. 

The control objectives for this example are 
defined as: (i) to make the closed-loop system stable; 
(ii) to make both Yy and Y> undisturbable by aT and E54 
(this will illustrate the design of a system for output 
undasturbabality) ici. 9section 3:9)? (iil) sto specizy 
complex eigenvalues and eigenvectors (just to illustrate 
that they can be handled directly); and (iv) to use a 
control law of the form u=Kx, plus feedforward control 
where necessary. 

Because the specification is for output undis- 
turbability and y is not a subset of x, it is necessary 


to transform the system. The transformation (cf. eqn. 


37,50) salS 
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Z=Q ° x = J[---|x = ae (Bn 2) 
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and the coefficient matrices of the transformed system 


defined by eqn. (3.51) become 


Seon LA = (3.73) 

ti cacy uae 

ENG Ge ae Tine 

1 

1.3 

asl aes 
1.3 
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Comparison of the system characterized by {a', p'} 
with the conditions in eqn. (3.30), reveals that state 


variable Z4=Yy is "almost undisturbable" with respect to 


disturbance ey since al =0 and al and a 


14 2re very 


at 
aba ia eles 

small relative to the other elements of A. 
However, feedback control will be required to gen- 


erate the conditions necessary for complete output undis- 
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Since eV) and al are both zero, the condition 


dg ae 


implied by eqn. (3.40) is already met by at 


and Y5 can both be made undisturbable to 1 without the 


and hence Yi 


use of feedforward control. However, feedforward control 
will be required to give ta the proper structure to make 
the system undisturbable to Eo: 

In the transformed system the number of inputs, 
r=3, and the specified number of undisturbable states, 
k=2. Therefore it is possible to arbitrarily specify 
all the eigenvalues and three elements of each eigenvector 


as follows: 


A = diagonal {-2+j, -2-j, -3.0, -4.0} (3.74) 


1 ae, eo 0 
| 
leet | ak ee x8) 0 
W = SSeS eae meee G3e37) 
O-5=400,5453) 1. 
{ 
| 
Wipe AO as) aA 


Note that the unspecified elements wig are to be calcul- 
ated as part of the design procedure and that the complex 
eigenvalues and elements of the eigenvectors must occur 
as complex conjugate pairs. 

It is obvious that W defined by eqn. (3.75) meets 
the condition for undisturbability specified by eqn. (3.39) 


and that the system is closed-loop stable. The feedback 
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controller matrix is then calculated from eqn. (3.53) 


and is shown in eqn. (3.77) below. The variables 2),Yj 


of the closed-loop system are undisturbable by aT when 
this feedback controller is used, but a feedforward con- 
troller must be designed to make them undisturbable by Eo: 
In order to give is the structure required by eqn. (3.40) 


it is noted that in this example: 


(o's 7.0) 


jie is sec equal "to {0,0,1.3,8) where the Iva is arbitrary 


then KEE and 8 can be calculated from eqn. (3.76). Note 


FP 


that the solution for K is exact. 


oF 


The feedback controller K is then transformed 


to correspond to the original system by using the relation- 


ship redler and the final controller matrices are: 


0.1342 0.4032 -1.3094 0.9886 
K =| -0.4469 025135 PaoisG =. 8509) (32/7) 


aes SIS )S: 2 O40 4.2935 saat 3) OH 3) 


er Onesie s0e | 0.69231 


Simulated Results: Figure 3.7 (left) shows the simulated 
response of Yy and Y> to a step change in aT under both 
open-loop and closed-loop conditions. As expected the 


open-loop response of Yi is essentially zero. The 
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unstable open-loop response of Y> is made undisturbable 
by the feedback controller and undisturbability of Yi 
maintained. Note that the closed-loop system has zero 
offset even though it uses proportional feedback control 
Only; is subjected to a sustained step disturbance; and 
has no natural integrating modes. 

Figure 3.7 (right) shows the simulated system 
responses to a step change in G5: The open-loop responses 
of Yi and Yo are improved considerably by the feedback 
controller but as predicted by theory both feedforward 
plus feedback control is needed to make Yi and Y5 undis- 


turbable by Go: 


Selo eConclLusions 


The system property of undisturbability has been 
formally defined and related to the concepts of uncontrol- 
lability, non-accessibility, and structural uncontrolla- 
bility. Necessary and sufficient conditions for undisturb- 
ability have been expressed in terms of both the structure 
of the coefficients matrices of the state-space model 
and also in terms of the structure of the system eigen- 
vector matrix. These results provide new insight into 
the structure of linear multivariable systems and have 
interesting design implications for the synthesis of 


regulatory controllers. Specifically the results allow 
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(i) determination of the property of undisturbability 
of specific input-output pairs by a simple inspection 

of the zero entries in the appropriate matrices, and 
(ii) they lead to a simple and constructive design pro- 
cedure for the design of regulatory controllers to achieve 
undisturbability. An important characteristic of this 
design procedure is that it makes use of the information 
contained in the model about the effect of disturbances, 
whereas other design schemes for feedback controllers do 
not. In comparison to existing design schemes for dis- 
turbance localisation the main advantage of the proposed 
design procedure is its computational simplicity. The 
main alternative design scheme besides requiring complex 
algorithms suffers from the hazards of numerical ill- 
conditioning. 

Undisturbability is closely related to the system 
property of uncontrollability, and also structural un- 
controllability but is not identical. It is also re- 
lated to observability. The controllers designed to 
produce output undisturbability can be thought of as 
modifying the closed-loop system so that disturbance (s) 
affect only the unobservable modes, i.e. those that do not 
contribute to the output variables of interest. 

Three numerical examples illustrate the applica- 
bility of the new theoretical results and the proposed 


synthesis procedure. 
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Transient responses of the closed- 
loop system in Example 3.1 toa 
unit step change in E41: 
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CHAPTER 4 


PRACTICAL APPLICATIONS OF MULTIVARIABLE CONTROLLERS 


DESIGNED TO PRODUCE UNDISTURBABILITY 


Abstract 


The design procedure for undisturbability devel- 


oped in Chapter 3 is applied to the design of feedback 


d and a order 


models of a double effect evaporator, an Se order model 


and feedforward controllers for the ee 


Obranelatray binary distillation column, andiva ote order 
model of another 8 tray, pilot scale, binary distillation 
column. Experimental application of feedback and feed- 
forward controllers designed to produce undisturbability 

of the main variables of interest in the computer con- 
trolled, pilot-plant evaporator gave results superior to 
conventional controllers. The problem of simultaneous dis- 
turbance localisation and asymptotic set point tracking 

is also discussed. Some results on disturbance local- 
isation using proportional plus integral feedback are 


also included. 
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4.1 Introduction 


In this chapter the design procedure for undis- 
turbability, developed in Chapter 3, is applied to the 
design of multivariable regulators for a double-effect 
evaporator. Feedback and feedforward multivariable reg- 
ulators designed to produce undisturbability were eval- 
uated experimentally by application to a computer-controlled, 
pilot-plant evaporator at the University of Alberta. A 
detailed description of the evaporator, its models, and 
a schematic diagram is included in Appendix C. The 
performance of the evaporator with these multivariable 
regulators implemented is compared with optimal-quadratic 
multivariable controllers, and a multiloop scheme that 
have been applied previously to the same evaporator. 

Each control system was evaluated by introducing 
a step change in one of the disturbance variables of the 
System. Three different disturbances can be used to upset 
the pilot-plant evaporator: feed flow, feed concentra- 
tion and feed temperature. Only the first two have been 
used in this study. Feed flow has been found to be the 
most severe disturbance in the system and this was the main 
reason for using it. The feed concentration disturbance 
has been used to evaluate the degree of interaction in 


the control system. The controllers were implemented on 
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an IBM 1800 digital data acquisition and control computer, 
which is interfaced to the pilot-plant evaporator, using 
a computer control package developed by Newell [1]. 

The simplicity and practicality of the design pro- 
cedure of Chapter 3 was also demonstrated by its applic- 
ation to the design of multivariable regulators for an 
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order models of two different distillation 
columns. The performance of these regulators was eval- 
uated by a digital computer simulation. 

Two areas that become important when actual 
applications were investigated are: (i) provision for 
making setpoint changes and (ii) possible addition of 
integral control to the disturbance localisation control- 
lers. While regulatory control in the presence of external 
disturbances is the prime control objective in most pro- 
cess control systems, asymptotic setpoint tracking or 
following is also an additional important requirement in 
an increasing number of supervisory computer control 
applications. In this chapter the problem of simultaneous 
disturbance localisation and asymptotic setpoint following 
is also considered. Because of the quasi-triangular 
structure of the closed-loop system matrix for undisturb- 
ability (cf. Theorem 3.3) and the assumption of non- 


singularity ‘or Bo) (cha esection 3.11); 1t is shown here 


a 
that the additional requirement of asymptotic setpoint 
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following is almost always satisfied. Using these results 
setpoint controllers were also designed for the double 


2 h order models of 


effect evaporator, and the reps and Bo4 
the two distillation columns. These controllers were evalu- 
ated experimentally by application to the evaporator, and 

by digital simulation of the distillation column models. 

Model inaccuracies, parameter fluctuations occur 
in most practical applications and thus the original de- 
Sign objectives may not be achieved, e.g. in the case of 
disturbance localisation the state variables of interest 
may no longer be undisturbable. In such situations the 
possibility exists that the introduction of integral feed- 
back of the undisturbable state-variables would compen- 
sate for the effect of model inaccuracies or unexpected 
disturbances, especially sustained disturbances. In 
this chapter such a conjecture is shown to be false. 

The present chapter is organised as follows: The 
results on simultaneous disturbance localisation and 
asymptotic setpoint tracking are first discussed in 
Section 4.2. In Section 4.3, the possibility of intro- 
ducing integral feedback of undisturbable state variables 
to compensate for the effect of parameter variations, model 
inaccuracies, and/or sustained unexpected disturbances is 
shown to be non-existent. The experimental evaluation 


of multivariable regulators designed to produce undisturb- 


ability is discussed in Section 4.4. The application of 
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the design procedure of Chapter 3, to an oe and 20" 
models of two different distillation columns is considered 
in Section 4.5 followed by concluding remarks in Section 


a. 0 


4.2 Simultaneous Disturbance Localisation and Asymptotic 
Setpoint Tracking 


The problem of regulation when a system is subject 
to external disturbances with simultaneous decoupling 
was first considered by Mclane and Davison [2], and more 
recently by Fabian and Wonham [3]. In [3] geometric nec- 
essary and sufficient conditions for simultaneous disturb- 
ance localisation, decoupling and stabilizability are 
given. In addition to [3], this problem has also received 
attention in [4]. The requirement for complete decoupling 
as well as disturbance localisation results in a set of 
restrictive necessary and sufficient conditions; steady- 
state decoupling is a suitably weaker requirement and 
yields a setpoint controller of constants which is easy 
to implement. Davison [5] has recently introduced the 
notion of steady-state invertibility and has derived nec- 
essary and sufficient conditions for a system to be steady- 
state invertible. This section is concerned with the 


investigation of the problem of simultaneous disturbance 
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localisation and steady-state decoupling. The main prop- 
erty of interest here is the steady-state decoupling of a 
quasi-triangular (or undisturbable) closed-loop system 


characterised by matrices C, B, H and L, where 


| 
1 
as | _ | 
{ — 2 
(Ara) 
Hy m0 0 
H = |-=~---| and L = }-- 
LE) Lo 


The output vector, y, for the above system is 
assumed to be of dimension m S r (r = number of inputs) 


and y is also assumed to be a subset of the state vector. 


mxr mxm 
e R e R r 


Further in) eqn. (4.1), 8 , Rank (B,) =m, Hy 


By application of 


uF 


p (n-m) x (n=m) and L, € 


H, € pect es 
Theorem 3.3 it is clear that the outputs of the closed- 
loop system characterised by matrices C, H and L in eqn. 
(4.1) are undisturbable with respect to the q disturbances 
bi bor auerety ago Let the system of eqn. (4.1) be governed 


by the following control law, 


u(t) = Kx(t) + yt (4.2) 


such that A + BK = H. The transfer function between the 
outputs, y(s), and desired setpoint values, Yq(S)e is 


given by: 
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C(sI - H)~*BK°Py. (s) 


y(s) 


(sI - H,)~*B,K*Py,(s) Coe) 


From eqn. (4.3) and the final value theorem, it 
is clear, for the class of inputs where the 


ae tey One (sy, (s)) is a vector of constant terms, that 
= 


y(t) = Yq (t) as t > ~ if: 


-1 SP (4.4) 


For expression (4.4) to be satisfied with the given 


SP 


assumptions, K has to be calculated from the equation: 


B.K =- H (45.5) 


With r =m and By assumed to be nonsingular (cf. Section 


oy De eae is uniquely defined by 


pees Ta (4.6) 


For r>m, eqn. (4.5) results in a set of (mxm) equations 


in (rxm) unknowns. Hence, there is no unique solution 


fox ieee One way to solve for Kee in such a case is to 


assign arbitrary values to (rxm - mxm) parameters in 


oe and then solve for the remaining (mxm) unknowns. 
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Alternatively out of r available control inputs one could 
choose the most suitable m inputs and thus solve for a 
unique cons For the case where r<m, i.e. the number of 
Outputs is greater than the inputs, eqn. (4.5) cannot 
be satisfied exactly. However, such a case is not very 
common in control applications. 

It is clear from eqn. (4.5) that for the more 
general case where r2m, and where undisturbability of 
the m output variables has been achieved, i.e. H, > 
exists by design, and Rank (By) =m, then simultaneous 
disturbance localisation and asymptotic tracking of a 
certain class of inputs is always possible. Design ex- 
amples to demonstrate degree of decoupling achieved 


using this procedure will be considered in Sections 4.4 


and 4.5. 


4.3 Disturbance Localisation Using Proportional Plus 
Integral Feedback 


Model inaccuracies and parameter variations 
occur in most practical applications and thus the orig- 
inal design objectives for which a state feedback con- 
troller may have been designed e.g. for undisturbability, 
may not be realised. In such situations the possibility 
exists that the introduction of integral feedback of the 


undisturbable state variables would compensate for the 
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effect of model inaccuracies, and/or parameter variations, 
and/or sustained unexpected disturbances. In the follow- 
ing analysis it is shown that theoretical difficulties arise. 

The design procedure for undisturbability as 
developed in Chapter 3 is based on eigenvector/eigenvalue 
assignment. Moore [6] has recently characterised the 
class of all closed-loop eigenvector sets which can be 
obtained with a given set of distinct closed-loop 
eigenvalues, by providing conditions which the closed- 
loop system eigenvectors must satisfy (cf. Section 2.2). 
In the ensuing analysis use is made of these conditions 
to show why integral feedback of the undisturbable state 
variables is not possible when H, = Ol. 

Integral feedback compensation will be considered 
by augmenting the original state vector x(t) with a @xl 


vector z(t) where z(t) is defined by: 


z(t) = T x(t) (4.7) 


fxn and it consists of 2 appropriate 


In eqn. (4.7) Te 
rows to include a subset of state variables or a linear 
combination of state variables requiring tntecrat feed- 
back. Note that because of controllability conditions 
on the augmented system, CN Let A and B now denote 
the augmented system state matrix and input matrix, 


respectively. Let the pair (A,B) be first partitioned 


in the £orm: 
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1 
Le 12 
aa 22 eS eees 
A51 1 Ayo (4.8) 
1 
i 
Aadndlan 353 
Where Adi E Rexk A55 Ss Ty E pixk and 
B gkxr ; 
1 E - If the state variables Kas Xor eeay x are 


to be made undisturbable then the structure of the 
closed-loop system matrix H and the corresponding 


eigenvector matrix W must be of the form (cf. Theorem 


By 3 eandic . 5): 

| Oo 

| 0 Wy | 0 0 

a Sa. eS ao ee 
Beri | Ape os ey ie .” anexl od 1 Woo. “oa 

1 | 
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Bee eels, BW a bi et ef 

(4.9) 


In (4.9) matrices H and W have been partitioned so as 
to be compatible with the partitions of A in (4.8). 

The design procedure for eigenvalue/eigenvector assign- 
ment requires that the closed-loop system eigenvalues 
be distinct which in turn implies that the closed-loop 
system eigenvectors must be linearly independent. This 


requirement will play a key role in the following 
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investigation. The main result of this study with 
respect to integral control is this: The elosed-loop 
etgenvectors of an undisturbable system wtth state plus 
integral feedback will not be linearly independent if 


Rank (T,) ZL, “where Be hae been defined on Ge. 67: 


2 
Remark 4.1: From an inspection of expression (4.8) 
the main consequence of the above result is this: If 
integral feedback of an undisturbable state variable, 
x, (ik), is required then the (i,i)* (isZ) element of 
T) should be equal to unity; however, for the design 
procedure to be possible the closed-loop eigenvectors 
must be linearly independent, i.e. Rank (TS) = £, which 


implies that the let 


row of T. should at least have one 
nonzero term in it. This means that integral feedback 
of x. with a linear combination of other state variables 
x, (k<jSn) is possible, and integral feedback of x, alone 
is not possible. 

The proof for the above result proceeds in the 
following way. For the eigenvectors in (4.9) to be lin- 
early independent, van must exist. Since W is quasi- 
triangular this implies that w, and w, must exist, 


where using the notation of Theorem 3.5, Wi and Wy are 


defined by: 
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Ww. = (4.10) 


Now for kSr, W., can be chosen arbitrarily and thus 


1 


existence of w,* Can be assured. For w, to exist it 


is necessary that the £ rows formed by [W55 W33] be 
linearly independent. To find conditions under which Rank 


[W = 2, consider the class of eigenvectors that 


opin icky 
characterise Wy by using the result of Moore [6], i.e. 


the a4 eigenvector, Wi, must satisfy the equation: 


il 
4 \ Ww. 
AT Adi Aj > 0 By 1 
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(4.11) 


In eqns. (4.11) w. corresponds to hs and M, | is defined by eqn. 
ay 


ub 
(eS ee Fork < 2. ne, Was 0 (cf. eqn. (4.9)). 
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This means that w,* and Ws for all k <i<Zn+ 2 must 


satisfy the equation: 
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THUS Lf Rank T, ZtL pethen Rank [W,. 
-1 


implies that Wy does not exist. This completes the 


W335] < £ which 


DEOGE* 

The important point to note is that T, can have 
no zero rows. This means that it is not possible to 
have independent integral feedback control of an undis- 
turbable variable by using eigenvector assignment tech- 
niques. For example, if x. (isk) is to be made undis- 
turbable and it is also desired to have an integral feed- 


back 2#of£ Xe alone, then this is not possible. However, 


n 
with xX undisturbable integral feedback of (x. + 2 x.a.), 


j=kt+1 J J 
with at least one oe #0 is possible. If H., is made non- 
zero (eg. to indicate parameter variations in Ay, By or K,) 


then it is possible to design an integral controller for 


x. alone. 
i 


4.4 Evaporator Application 


The design procedure for undisturbability devel- 
oped in Chapter 3 was evaluated by application to the 
computer-controlled, pilot-plant evaporator at the 
University of Alberta. This evaporator has been used in 
a number of other studies to evaluate different modern 
multivariable control techniques (Fisher and Seborg, [7]) 
and hence it is possible to make direct comparisons be- 


tween alternative control techniques. 
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A description of the pilot-plant evaporator and 
its models is given in Appendix C. A schematic diagram 
of the evaporator and the conventional multiloop control 
scheme used as the 'base case' for comparison is also 
shown in Figure C.1, Appendix C. The primary control 
objective is to keep the product concentration, C2, 
constant despite disturbances in the feed flowrate, F, 
the feed concentration, CF, and/or the feed enthalpy, 
HF. It is also necessary to keep the two liquid holdups 
Wl, and W2, within operating limits, but small varia- 
tions in these variables are acceptable. The control 
(manipulated) variables are the steam flowrate, S, and 
the bottom product flowrate from each of the two effects, 
Bl and B2. In summary, the evaporator has 3 output var- 
iables, 3 control variables and 3 disturbance variables. 
These variables and their normal operating steady state 
values are defined in Appendix C. 

A number of different models of the evaporator 
have been derived in previous studies (e.g. Newell and 
Fisher, [8]). The models used in this study are the 


rd h 


3 and se order, discrete, state-space models defined 


by the coefficient matrices listed in Tables C.2 and 


Coo. “The 3rd order model is not as accurate as the 
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5 and es order models used in other studies, but is 


in reasonable agreement with experimental data and makes 
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it much easier to follow the design procedure. Hence, 


the design of a controller to produce undisturbability 


d 


will be) first considered for this vie order model. 


4.4.1 Controller designe torsethe order model 


The ate order model has as its state variables, 
Wl, W2 and C2; as its input variables S, Bl and B2; and 
as its disturbance variables, F, CF, and HF. Previous 
Operating experience has shown that the most frequent 
and severe disturbances in product concentration, C2, 
are produced by variations in feed flowrate, F. There- 
fore for purposes of this application the design ob- 


jectives for the controller are: 


i) to make C2 undisturbable by F. 
ii) to assign all closed-loop eigenvalues closer 
£o the origin. 
iii) to preserve the open-loop undisturbability of 


Wl and W2 by CF. 


The linearised model in the discrete form as 


used in this study can be represented by: 
x(n+1) = x(n) + Au(n) + @d(n) 


with the coefficient matrices, ®,A and 8 and elements of 


x, u and d given in Appendix C. 
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According to the gee 


order model, the evaporator, 
under open-loop conditions, is only conditionally stable 
Since the eigenvalues of @ are 0.96, 1 and 1. The ate 
unit eigenvalues in 6 follow directly from the "integrat- 
ing nature" of the two liquid holdups. Inspection of 
the coefficient matrices for the gate order model indic- 
ates that under open-loop conditions, all three state 
variables are disturbable with respect to F and HF, 
while Wl and W2 are undisturbable with respect to CF. 
(AS noted in Remark 3.7 of Chapter 3, the results in 
Section 3.7 can easily be extended to discrete-time 
systems.) 

Since eet #0 it is not possible to make C2 


undisturbable with respect to F by feedback control alone. 


However, a feedback feedforward control law of the form: 


F 


u(n) = K x(n) + K Pe(n) 


is required. In the discrete form of the closed-loop 
system equation, L = 6 + eee and hence it is possible 
by feedforward control to produce the required structure 


Treas Setting ee = la, 6, Ol with Cand erarbitraraly 
FFL 


chosen as 0.12 and 0.0032 gives, K as: 
0.32 
eae 
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Note that if a and 8 were chosen as zero, than all three 
State variables would be undisturbable with respect to F. 
This is possible because in this particular example A 

is non-singular and consequently the equation, Ly = ory “le 
i has an exact solution for aes To complete the 
controller design, the feedback matrix K is calculated 
using the eigenvector/eigenvalue assignment procedure 
described in Chapters 2 and 3. Since in this example, 
r=n, it is possible to assign all elements of the eigen- 
vector matrix W. If W is chosen to be the same as the 
Open-loop eigenvector matrix (a diagonal matrix since ®? 


is diagonal), then the following design objectives are 


realized: 


i) C2 is undisturbable with respect to F since W 


and L, are inethe form of eqns. (3.539) “and (320).. 


ii) Wl and W2 are undisturbable with respect to CF 


(as was the case for the open-loop system). 


Thus choosing W to be the open-loop eigenvector matrix, i.e. 
W = I, and assigning the closed-loop eigenvalues to be 
Qo2c, Ona7 and 0265 leads to the following state feedback 
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2-70 0 = 97.09 
K =" | oes 0 3.86 (4.13) 


4.29 13.05 Ley Bike, 


Table 4.1 contains a summary of this and other 


controllers used in the experimental evaluation. 


4.4.2 Signal flow graph analyses 


The signal flow graph in Figure 4.1 (left) shows 
that in the simplified, epee! order evaporator model, all 
three states are accessible to a disturbance in feedflow, 
F, but only the product concentration is accessible by 
a disturbance in feed concentration, CF. In other words 
Wl, W2 and C2 are disturbable by F and Wl and W2 are 
undisturbable by CF. 

Figure 4.1 (centre) is the signal flow-graph of 
the closed-loop evaporator system using the disturbance 
localisation controller discussed above. Comparison 
with the open-loop system flow graph shows that the 
feedforward controller eliminates the path from F to C2 
and hence makes C2 inaccessible and undisturbable by F. 

Figure 4.1 (right) is the signal flow graph of 
the evaporator system using the basic multiloop control 


scheme shown in Fig. C.1l. The multiloop controller, 
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Ke! used in this comparison is shown in Table 4.1. 
Note that in using this controller all states become 
accessible to both disturbances, and that a number of 
interactions are introduced when previously there were 
none. This emphasizes the importance of examining 


the structure of open and closed-loop system matrices. 


4.4.3 Experimental results 


Figure 4.2 shows the response of the evaporator 
to , 30% step changes in feed concentration, introduced 
at the times indicated by the arrows on the time axis. 


Controller, was designed to retain the open-loop 


Kae’ 
undisturbability of the liquid levels to concentration 
changes and the results in Figure 4.2 confirm that there 
is)no -significant ettect of CF on Wl or W2. ihe con— 
troller was not designed to make C2 undisturbable by CF 
and a small offset in C2 is obvious in Figure 4.2. Note 
that this controller is an ideal modal controller since 
the open-loop and closed-loop eigenvectors are ident- 
aca. 

The evaporator response to + 30% changes in 
feed concentration using the multiloop controller, Ke! 


is shown in Figure 4.3. This experimental data confirms 


the conclusion derived from the signal flow graph analyses, 
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that the multiloop control scheme introduces significant 
interactions. 

Figure 4.4 shows the evaporator response to +20% 
step disturbances in feed flowrate using the feed- 


FFL 


back plus feedforward controller, plus K . Figure 


Kael 
4.5 shows the evaporator response to +20% step disturb- 


ance in F using feedback control, K alone (equivalent 


a gay 
to the case where disturbance measurements are unavail- 
able). Another controller, Kapoe which retains the same 
Open-loop eigenvector structure but assigns eigenvalues 
at 0.59, 0.43 and 0.34 was also designed. The response 
of the evaporator with a +20% disturbance in feed flow- 
rate using only the feedback controller Kap is shown in 
Figure 4.6. Comparison of Figures 4.5 and 4.6, and 

and K 


controllers K serves to illustrate the design 


Vail Qie2 
freedom available, i.e. arbitrary eigenvalue/eigenvector 
assignment. 

Figures 4.7 and 4.8 show the evaporator response 
to the same disturbance using (i) the multiloop con- 


Erolver, K and (ii) the feedback plus feedforward 


me 
FF2 : 
controller, Kil and K as shown in Table 4.1. Also 


listed in Table 4.1 is a third order optimal proportional 


controller, For the purposes of comparison with 


Koptl’ 


Figures 4.3 to 4.8, two experimental runs were performed, 


using this optimal controller. The results for a +20% 
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disturbance in feed flowrate are plotted in Figure 4.9 
and Figure 4.10 shows the evaporator response to a +30% 
disturbance in feed composition. 

The response of the evaporator to a +20% disturb- 
ance in feed flowrate or a {30% disturbance in feed com- 
was superior 


position using controller, K or K 


all ak2 
than the performance of the multiloop controller, Ke! 
and better than,or at least comparable to, the perform- 


ance of optimal controller, K 


Opt 
With the feed forward control mode implemented 
the overall performance of the controller K sel plus es 
is better than the performance of Kel and roe It should 
be noted that when implementing feedforward control with 
the multiloop control configuration, Ly = 8 or ene =—FOs, 


that is the signal flow paths from F to Wl, W2 and C2 
have all been '‘'cancelled out" by the effect of control. 
In contrast to this the feedforward control action in 


conjunction with K Ors requires that only the sig- 


agli ax2 
nal flow path from Ff to (C2 be ‘cancelled out'. The 


net result of this is that feedforward controller, oie. 


With the multiloop control configuration must generate 


more control action, i.e. use higher gains than the 


feedforward controller, eae This is also evident 


from a comparison of the steam, S, to feed flowrate, F, 


Pie F2 


gains in K and KF ‘ 
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In general, the design procedure for ‘undisturb- 
ability' was found to be practical, was easy to use, 


and performed well experimentally. 


4.4.4 Controller design for the 5th order model 


h order model of 


The main reason for using the ee 
the evaporator was to demonstrate the use of the design 
procedure of Chapter 3, for the general case where the 
number of inputs, r, is less than the number of state 
variables, n. Application of the design procedure to 


the 5th 


order model is also useful in demonstrating the 
use of the stability criterion for disturbance locali- 
sation, i.e. to check if the closed-loop system is 
Stable (cf. Section 3.11.2). A discrete hae order 
model of the evaporator is shown in Table C.2. ‘The 


a order model are: Wl, Cl 


state variables for the ue 
(first effect composition), Hl (first effect enthalpy), 
We-end C2. ‘The control inputs, disturbances and the 


d order model. 


Outputs are identical to those for the oh 
The main control objective as before is to make C2 and 
W2 undisturbable with respect ot disturbance F and also 
to retain the open-loop undisturbability of Wl and W2 to 
CF. | 


As a preliminary first step the state vector x 
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was reordered as x! = [eo wow CL, HIT” and the state, 
control and disturbance matrices were permuted accord- 
ingly and designated by oe nN and Q. Since the fiest 
three rows of A are linearly independent it is possible 
to assign arbitrary entries in the first three rows of 
the closed-loop system eigenvector matrix using the 
eigenvalue/eigenvector assignment technique. In order 
to realise the above mentioned control objectives the 


and W were selected: 


following values of Wy 5 


ul 0 0 0 0 
Wa = |0 ih 0 and Wo = 10 0 (4.14) 
0 0 dl 0 0 


Observe that since W, = 0 it is not possible to assign 
eigenvalues corresponding to these two eigenvectors (cf. 
Section 3.11.2). This necessitates use of the stability 
criterion to check if the resulting closed-loop system 


is stable. Equation (3.57)shows that for this case the 


matrix P is: 


0.90 0 
— (4.15) 
-0.31 -0.76 


and thus has eigenvalues: 0.9 and -0.76 which are inside 


the unit circle. This means that two of the closed-loop 
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system eigenvalues are -0.76 and 0.9 while the remaining 
three can be chosen arbitrarily as 0.28, 0.47 and 0.65. 


The final closed-loop eigenvector matrix W was evaluated 


£nom tens. (3..55).and o(3356)..and is 


1 0 ie 0 
| 
0 1 0.4 0 
: | 
w=! 0 0 a =O (4.16) 
{ 
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W as shown above is non-singular and in the re- 
quired form for undisturbability (cf. Theorem 3.5). The 


feedback controller K 5’ was then calculated using the 


ae 
pseudoinverse of A and is shown in Table 4.1. The re- 


sulting closed-loop system matrix is: 


(4.17) 


Since eo, does not have zero entries in appropriate 


places a feedforward controller is also required to 


complete the design requirements for undisturbability. 
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For the purposes of illustration, the design of two dif- 


EES FF4 


ferent feedforward controllers K and K was consid- 


ered in this application. ae was designed to make C2 


and W2 undisturbable with respect to F (i.e. a = 


F4 


LapOeasd #6 yi )omond ie was designed to make all 


three output variables C2, W2 and Wl undisturbable with 


~ 


respect. to F (i.e. a = (0 0 0 v yj). These two 


feedforward controllers are defined in Table 4.1. 


4.4.5 Experimental results 


Figures 4.11, 4.12 and 4.13 show the evaporator 
responses to a +20% step disturbance in feed flowrate 


with feedback controller, Kaps’ implemented, with Raps 


PES 


plus K implemented, and with Kaps plus lee implemented. 


For comparison purposes the evaporator response to a 


h 


+20% disturbance in feed flowrate using a we order pro- 


portional optimal controller is also shown in Figure 4.14. 


The evaporator performance using any one of these se 


order state feedback controllers implemented in each of 
the four experimental runs is better than that obtained 


using a 379 order controller. This would be expected 


rd 


because the 3 order model is not as accurate as the 


eee order model. From an overall performance point of 


view the evaporator response to feed disturbances using 


a sath ah RS ds 
Lis stem Ox PONEsma, see BO pay “tty 
dec! &faetqoseLbre tot Bee ie er eancat 


as yaar a tw ee O° G1 + 2 oe ty Toe 


‘ties bamiktob aye 0 oo Buea 
-_ 


IOIEACWES Dis YOR: jPeivheas LI k fied aennplt ee 
bars? thw sedaysuce sapere: ODER, Apert 


x) aoke DagnEmns qe a M9 FEO ood Nae 


ASHRAM, t 


256 
Bass: bits ‘ive pee tite’ .oe2 


ipa orienta Aaa 
. gd Senid?st cit RaQueVE oe speog wy. a 

-924 24075 or, » pRseH Ti | BeSh ad, 
Sty supe of awowe wel inn telblosgnen | 
8 ue. Ast tO Ste a sie sonentot seat 


SG “ifce (il he QOS Lane i 


hemisede cn baie 


feedback controller, Kaes! plus feedforward controller, 
Rees or that obtained using only feedback controller, 

Kapse is at least comparable to or even better than the 
proportional optimal control. Another Bae order con- 


troller, as designed by using a simultaneous 


Rie’ 
eigenvalue assignment procedure (Park [9]) was also 
implemented on the evaporator. Some difficulty was 
experienced in implementing this controller. Its per- 
formance as shown in Figure 4.15, was found to be 
escillatory: Difficulty was experienced in implement- 
ing other controllers designed by the same method. A 
3rd 


order modal controller, based on Rosenbrock's 


ley 
approximate modal control method [10] was also implement- 
ed on the evaporator. The performance of this controller 
to step disturbances in the feed flowrate, as shown in 


Figure 4.16, was satisfactory but not as good as the 


response of controller, Kei! as shown in Figure 4.5. 


4.4.6 Design of setpoint controllers for the double- 
efleéct evapordcer 


The 3xd 


order evaporator model has 3 state- 
variables, 3 inputs and 3 outputs. Because of the dis- 
crete nature of the evaporator models, with only minor 


modifications, the results of Section 4.2 can be applied 
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directly to calculate the setpoint controller, K ; £OL 


d 


the discrete 3*° order evaporator model. The condition 


for asymptotic setpoint tracking or following in dis- 


d 


crete systems for the 3°° order evaporator model 


requires that: 


(ft one eee a he ST (4.18) 
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The resulting, K , for the 3°° order model is listed 
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in Table 4.1. For the discrete eo order evaporator model, 


I : : 
K° 2 1S given by: 


al 
1 


DPS 


K = A (I - Hy) (4.19) 


~ 


where Hy Ls. the. upper detti3ex 3) matrixes partition.of 


~ 


aes) as shown in eqn. (4.17); A, is the 


Upper Syxieumatrixesparts tlonsof Ar with) its finst three 


H(= 6 DK 


rows corresponding to state variables C2, W2 and Wl. 


The resulting setpoint controller, ee for the pth 


order model is listed in Table 4.1. The evaporator 


response to a +10% step change in the setpoint of C2 


SP1 SPs 


with controllers. Kk plus K , and Kaes plus K 


acti 
implemented is shown in Figures 4.17 and 4.18. As 


expected the evaporator response to +10% step change in 


the setpoint of C2 with controller Kips plus Sas was 


better than that of the oe order controller, Kael plus 


eu It should be noted from Figures 4.17 and 4.18, 
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that in each case the asymptotic setpoint tracking 


property of controllers ae and ee 


is approximately 
satisfied. For the purpose of comparison, the evapora- 
tor response to a +10% step change in the setpoint of 
C2 using an optimal controller, Boel plus aga was 
also evaluated experimentally and is shown in Figure 
4.19. The evaporator response in Figures 4.17 and 
4.18 is better than or at least as good as that in 
Figure 4.19. 

In this section results on experimental eval- 
uation of the design procedure for undisturbability 
have been presented. In general the design procedure 


for undisturbability was found to be practical and con- 


venient. It was easy to use, gave considerable insight 


into system performance, placed the closed-loop eigen- 


values and eigenvectors in the desired locations, pro- 
duced practical controllers with reasonable gains, and 
performed very well experimentally. To the best of the 
author's knowledge, this is the first reported applic- 


ation of this design approach. 


Experimental evaluation of controllers, Kael plus 


Sel 


K , and K plus Boe also confirmed the results of 


ats 


Section 4.2 that simultaneous disturbance localisation 


and asymptotic setpoint following is possible. The 
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design procedure for calculating such setpoint controllers 
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is easy to use and gives easily implementable constant 


Setpoint controllers. 


4.5 Distillation Column Application 


1 loge bl yeh order distillation column model of Davison [11] 


A state space model of an 8 tray binary distilla- 
tion column with variable column pressure has been de- 
rived by Davison [11]. The linearized equations of the 

eh 


acd § order system with 3 inputs, 3 outputs and 1 main 


disturbance are as follows: 


x(t) =A x(t) + B u(t) + Dw (t), y(t) = C x(t) 


(4.20) 


The state vector, x(t), and the output vector, y(t), are: 


44 
x 


It 


(XOr Xyr Xor seer Xgr Xe p) and 


y = (Xor Xp Pp) 


with oes composition of more volatile component in con- 
denser (mole fraction); 2S ae composition of more vola- 
tile component on plates, i = 1,2, ..., 8 (mole fraction); 
x, = composition of more volatile component in reboiler 
(mole fraction); u, = reboiler temperature; u, = conden- 


ul 


ser temperature; u, = reflux flow (lb-mole/sec); p = 
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pressure in the system (atm) and Mes feed composition. 
Numerical values of matrices A, B and D are given in 
Bate: 

According to Davison [12], the main control ob- 
jective in the design of a controller for the column is 
to regulate the top (vy) and bottom product (y5) com- 
positions, and the pressure (y,) in the presence of com- 
position fluctuations in the input feed stream. The 
open-loop distillation column is stable. Inspection of 
Matrices A, B and D reveals the following: 

i) Elements of D corresponding to state variables 


X_, X, and p are zero, hence no feedforward con- 


= b 
trol will be required in the design for undis- 


LULDabD EL ty. 


ii) The row of B corresponding to state variable Xo 
has all zero elements. Consequently the corres- 


ponding partition B, is singular. This means 


i 
that the design technique for undisturbability 
Cannot be applied directly,and hence the control 


objectives have to be modified. 


iii) In modifying the control objectives use is made 
of the following important characteristic of 
the system: that is since x, cannot be in- 


fluenced dtrectly by any one of the three inputs, 
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the corresponding row of A will also be the row of the 
closed-loop system matrix, H. Inspection of the first 
row of A corresponding to xX, also reveals that, Xe is 
accessible to the control cwpnte via a coupling 


between Xo and Xy° This is important because it means 


that ie xX is made undisturbable then Xo will also be- 


come undisturbable, especially since the first row of A 
is invariant under any feedback. Thus the modified 


control objective will be to make x x. ands p.undis-— 


Le 0 
turbable with respect to Wa Then once xy becomes non- 


accessible to we undisturbability or nonaccessibility 


of x, to Wr can also be assured. Observe that the first 


~ 


three rows of B, as denoted by B, corresponding to the 


dl 


modified system with x = [x,, Xr pr Xor Kor Xagreees Xo] 
are now linearly independent. The rest of the design 


procedure is now simple. With a choice Wy = T,, and 


Wo = Ouellet. Theoremi3.5)', and Ay 


a check is made to see if the eigenvalues of matrix P 


=) 0180 .-he .o7ene oy see, 


(c£. eqn. (3.57)) are in the left half of the complex 
plane (Section 3.11.2). The eigenvalues of P do indeed 


turn, out, to, be, in -thesteft) half of the: complexs plane. 


With this check complete, K is very simply evaluated from 


the equation: 


tae o 4 e 
K =B [H, | O]), with H, = A 
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Notice that computation of K only requires inversion of 
a 3X3 Matrix, and matrix multiplication. The computed K 


was found to be: 


-4271.6 2975s (4 soe 1 /6.5. =00.7 °°") Omeo Oet0 S01. 5 


ze 
ll 


4271.6 -3031.1 -147.2 Zed 68.7 700 RO 080s 301 
=37/ 56 -54.5 cag =e =0.6 (090 20F70,0 -5.4 
a2) 


The resulting closed-loop system matrix, H = A + BK and 


the disturbance matrix D are shown in Table 4.2. Both 


~ ~ 


H and D are in the required form for x xX P and x. 


ees 
to be undisturbable. Figure 4.20, shows the response 


of Xr x and p to a unit step disturbance in wae Also 


b 
superimposed on each of plots in Figure 4.20, is the cor- 
responding open-loop system response to a step disturb- 
ance in war For the purpose of comparison an attempt was 
also made to reproduce the response of the distillation 
column with a conventional controller and another con- 
troller designed by Davison [12]. However, diffi- 
culties were encountered in reproducing these two 
results, and hence each plot in Figure 4.20 has an insert 
from [12] showing the nature of these responses... For, 
convenience the controllers depicted as 'conventional 


control' and 'new control system' in the figure caption 


in the inserts are also reproduced here: 


* Use of the numerical values in [12] led to unstable 
responses. 
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New control system [12]: 


c 
II 


1 = 000 xX» + 600 X57 1000 .P >=" 4552 Whe 1 2010 we 


Uy 1000 X5 600 Xo LOOO RE sez Wh arink 49900 Wes 


G 
II 


= 1000 x. =) 1000 x) — 286 w. + 64900 w 
3 2 7 m m 


Conventional control system [12]: 


uy 100,000 x4 


u, = -100,000 p 


Notice that Davison's controller requires feedforward 
COntErOL of We and also computation of Win 

Using the results of Section 4.2, a setpoint 
controller, xSP was also designed for the above distil- 


lation column. This was calculated to be: 


Miss 2 = aystea 9 272.6 
pf eide7e. 3 Bats) 162.2 (4.23) 
128.3 67.2 vo) al 
~sp 


Implementation of K as shown above results in an 
almost completely non-interacting system as evident by 
the response of the distillation column to unit step 
setpoint changes in each of the variables Xun Xp and p 


(Figure 4.21). 
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4.5.2 204 order distillation column model of McGinnis 


[3] 


McGinnis [13] has derived a 20th order linear 
state space model of a pilot scale distillation column 
at the University of Alberta. This model has been de- 
rived by considering stage enthalpies and liquid flows 


nd order model 


at each of the ten stages. A reduced 2 
has also been obtained from this higher order model, 

and has been used in previous control applications, 

Wood and McGinnis [14]. The 20th order model has 2 
inputs, 2 outputs and 3 disturbances. The outputs of 
interest are the liquid phase enthalpy in the reboiler, 
h, (=yj), and the liquid phase enthalpy in the condenser, 
hy (=y19)- The 3 disturbances of interest are feed 
flowrate, feed enthalpy and condenser heat duty. Of 
these, feed flowrate has been found to be the most severe 
disturbance and hence the main control objective in this 
Study is to make the two state variables or the outputs, 
hy and Aig: 


ance. The 2 control variables are reflux flow (R) and 


undisturbable with respect to this disturb- 


Steam flow or reboiler duty (QR). Numerical values of 
matrices A, B and D are given in [13]. 
As would be expected feed flowrate disturbance 


directly affects only the feed stage, i.e. state variables 
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hn, and Le (liquid phase enthalpy and liquid flowrate 

at stage 5), and then propagates to all of the remaining 
state variables by dynamic coupling between each of the 
stages. It is then clear that as far as feed flowrate 
and enthalpy are concerned the rows of D corresponding 

EO hy and hyo have zero entries. In other words feed- 
forward control action is not required for undisturbabil- 
and h 


ey Oe. with respect to feed flowrate and feed 


h 


1 10 


enthalpy. Just like in the Wins order model, the 


row of B corresponding to state variable hig has two 


zero elements. Consequently, the partition By is sing- 


ular. However, the effect of control directly reaches 


hy and since there exists a signal flow path from Ng to 


Ayor the original control objective is modified to make 
hy and Ay undisturbable. It was known that in carrying 
Out such a strategy, the original control objective would 


be realised because the row of A corresponding to hig is 


invariant under any feedback. The design procedure for 


undisturbability was then applied to the augmented system 


~ 


Wiehe = [eee Opa ees. oe ah L 


g? Migr Agr Agr Se eee aig ae 


With a choice W, = Io, Ww. = 0, and is diag(-2.5, -16.0), 
eigenvalues of the matrix P were evaluated (cf. eqn. 
(3.57)). The eigenvalues of P were found to be in the 
left half of the complex plane. With the check for 


~ 


Stability complete, K was evaluated from the equation: 
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K=B, ({A, Aj] - [H, 01), with H, = A, 
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(4.24) 


Notice that the computation for a state feedback con- 


troller for a 20th 


order system only requires inversion 
Of a 2x2 diagonal matrix and its multiplication by a 
Sparse 2x20 matrix. The point to note is that the de- 
Sign procedure for undisturbability is very simple and 
easy to use, and as illustrated by the above example it 
can often be carried out without the aid of a computer, 
even for such large order systems. A setpoint con- 
troller based on results of Section 4.2 was also de- 


signed for this system. The final control law for the 


above system was found to be 


Riars 9Sbh este adio Lnwereel 6.29 hes=nayess h,”? 
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eae paach ie Cech oe eee my keh ieee 
QR = : ioe’ 2 : 2 : 10 

(4.25) 


The response of hy and hig to unit step disturbances in 


feed flowrate with the control law of eqn. (4.25) im- 


plemented is shown in Figure 4.22. The performance of 


the setpoint controller was also evaluated by application 


of unit setpoint changes in hy and hio- The remarkably 
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large degree of non-interaction achieved with a simple 
asymptotic setpoint controller is also evident in 
Figure 4.23. 

A remark on the practicality of the control 
scheme of eqn. (4.25) is in order here. State variables 
h h 


h h, and h are directly proportional to the 


LS os aac: Baas: 10 
temperatures at each of corresponding stage and hence 
these are eaSily available for state feedback. State 


Variable, L represents the liquid flowrate from the 


2" 
bottom tray to the reboiler. This measurement is not 
available at the present time on the pilot scale dis- 
tillation column. However, a reasonable estimate of 

Los can be obtained by writing a material balance equa- 
tion around the reboiler, and assuming that the vapour 
flowrate is directly proportional to steam flowrate into 
the reboiler. It is also important to emphasize here, 
that by using the structural aspects of the results on 
undisturbability, the form and structure of the control 
scheme (such as in eqn. (4.24)) for this or a similar 
system could have been predicted with a knowledge of 
only the zeroelements of matrices A, B, and D. Such 

a preliminary analysis prior to actual design and 


modelling can be useful in planning of sensor locations. 


It should also be noted that the control scheme 
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shown in eqn. (4.25) requires only measurement of 6 
State variables. The reason for this is that because 
of the sparse open-loop system matrix ne and the zero 
elements of Hy as required for undisturbability, a 
large number of columns of the difference matrix 


(TA, Ay] [H, 0]) are zero (cf. Section 2.4). Thus 
from eqn. (4.24) it is easy to see that feedback of 
state variables, corresponding to zero columns of the 


difference matrix, is not required. 


4.6 Conclusions 


The design procedure for undisturbability 
developed in Chapter 3 was applied to the design of 
feedback md feedforward controllers for the computer 
controlled, pilot-plant evaporator, an ye order model 
of an 8 tray binary distillation column, and a 20 order 
model of an 8 tray pilot scale binary distillation 
column. The problem of simultaneous disturbance local- 
isation and asymptotic setpoint tracking has also been 
discussed here. The proposed solution to this problem 
is both practical and simple. The possibility of having 
an independent integral feedback control of an undis- 
turbable variable was also investigated and such a 


possibility was shown to be non-existent. 
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Experimental results from the computer-controlled 
pilot-plant evaporator at the University of Alberta 
demonstrated that the use of design procedure: produced 
the expected undisturbability; provided better control 
than conventional single variable controllers; and the 
performance of the resulting controllers was at least 
comparable if not better than controllers designed using 
other modern multivariable control techniques such as 
optimal-quadratic control. Evaluation of the controllers 
for the two distillation column models by digital sim- 
ulation also demonstrated the effectiveness of the de- 
sign method in eliminating the effect of specific dis- 
turbances on outputs of interest. The computational 
requirements for the design procedure were shown to be 
Only marginal even in its application to large order 
systems. Experimental evaluation and digital simulation 
of setpoint controllers designed to retain the disturb- 
ance localisation property and at the same time asymp- 
totically track a certain class of inputs gave good 
results and confirmed the practicality of such an 
approach. 

In general, the design procedure for undisturb- 
ability was found to be practical and convenient. It 
was easy to use, gave considerable insight into system 


performance, placed the closed-loop eigenvalues and 
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eigenvectors in the desired locations, produced practi- 
cal controllers with reasonable gains, and performed 
well experimentally. To the best of the author's 
knowledge, this is the first reported experimental 


application of this design approach. 
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in we (feed composition). 


155 


Ve (ncnt-oangla) St ¢ 
| Como t <mvenh) 78” 1 


~*~ 
~~ 
Ww 
. 
~, 
A 
‘| - 

"1 od eS EM 9 
¢ 7 ' ie 18 4 ‘TM —_ —- = - “aw 
{ | io) cA S : —_": 

in L004 / 
a 
1” 
% 
po Ae Se > ep" ri al 
Pdi ee 
aii 


% a naenety  e ha : ~ 
\ a 


Lael whan) 
(facl~sahor NY ve 


Figure 4.21: 
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Figure 4.22: Simulated responses of the ayers order dis- 


tillation column model comparing closed- 
loop and open-loop response of hy (top plot) 


and hyo (bottom plot), to unit step dis- 


turbance in feed flowrate. 
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Figure 4.23: 
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Simulated response of the a order dis- 
tillation column model (with feedback plus 
setpoint controller), to unit step setpoint 
changes in h, (top plot), and h (bottom 

1 10 
plot). 


CHAPTER 5 


DISTURBANCE MINIMISATION AND POLE ASSIGNMENT IN LINEAR 


SYSTEMS WITH APPLICATIONS TO OBSERVER DESIGN 


Abstract 


A method is presented for designing a state feed- 
back controller to assign eigenvalues in a closed-loop 
system and also minimise a quadratic performance index 
which includes a measure of coupling between disturb- 
ances and the outputs. Two illustrative examples are 
provided to demonstrate the use of the method. An 
experimental application of the proposed procedure was 
carried out on a pilot plant double effect evaporator 
and this evaluation is also included here. An applic- 
ation of the method to the design of full order observ- 
ers for systems with unmeasurable disturbances, is also 


considered. 
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5S. Introduction 


Practical control systems must operate in environ- 
ments where unknown external disturbances are invariably 
present. Despite the inevitability of external disturb- 
ances, few multivariable synthesis techniques make use of 
the information usually available in process models that 
describes the effect of external disturbances. For example, 
optimal control systems in which a quadratic performance 
index is minimised, or conventional pole-placement tech- 
niques typically ignore disturbances or only consider 
impulse-type disturbances. Consequently such systems 
give satisfactory control when the disturbances occurring 
in the system are of an impulsive-nature, but may not 
perform well when the disturbances are sustained, which 
unfortunately is frequently the case. 

The present chapter poses and deals with the fol- 
lowing two questions: 1) Is it possible to formulate 
the regulator problem from an optimisation point of view 
such that the effect of arbitrary and unknown external 
disturbances on system outputs of interest is minimised? 
2) How does one construct a suitable feedback strategy 
which provides maximum attenuation of disturbances in the 
resulting closed-loop system and which ensures that the 
system eigenvalues are at specified locations in the left 


half of the complex plane? 
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The main difference between the approach adopted 
here and the related results o disturbance localisation 
(or rejection) by Wonham and Morse [1], Bhattacharyya [2] 
and in Chapter 3 (see also Shah et al. [3]) is that the 
results derived here can be applied to any linear, con- 
trollable multivariable system, i.e. even systems which 
do not satisfy the necessary and sufficient conditions 
for disturbance localisation. 

This chapter is organized as follows: Section 
5.2 provides an outline of related previous work. Section 
5.3 considers the formulation of the performance index 
or objective function to be minimised in a general form, 
namely, a functional which is a linear combination of 
(1) a measure of coupling or transmittance between the 
disturbances and the outputs, and (ii) a quadratic per- 
formance index. Section 5.4 shows how stability and/or 
arbitrary pole assignment constraints can be incorporated 
into the above formulation. The design algorithm is also 
outlined in this section. Section 5.5 contains two 
numerical examples to demonstrate the use of the proposed 
design procedure. The application of the design procedure 
to the design of observers, so that the effect of unmeas- 
urable disturbances is minimised, is the subject of dis- 


cussion in Section 5.6. 
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5.2 Previous Work 


Wonham and Morse [1], and Bhattacharyya [2] have 
outlined geometric necessary and sufficient conditions 
for the existence of state feedback control to eliminate 
the effect of external disturbances on system outputs. 

In Chapter 3,the term undisturbable refers to a system 
in which selected state variables are invariant to arb- 
itrary changes in a specified disturbance, and necessary 
and sufficient conditions for undisturbability are ex- 
pressed in terms of the structure of the coefficient 
Matrices in the state space model. 

By a geometrical characterization of the regulator 
problem, Bhattacharyya et al. [4] have derived conditions 
under which an initially zero output of a system is main- 
tained at zero, and any non-zero output due to non-zero 
initial conditions or impulse-type disturbances is made 
to approach zero. An extensive treatment of the regula- 
tor problem based on an algebraic approach can also be 
found in Silverman and Payne [5]. Johnson [6,7] has 
also studied, from an optimiSation point of view, the 
regulation problem in the presence of a class of unknown 
disturbances that can be characterized by certain waveform 
types or modes. 


Davison and Smith [8] have derived necessary and 
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sufficient conditions for the existence of a minimal 
order, realizable, state feedback controller such that 
the eigenvalues of the closed-loop system take on pre- 
assigned values in the left complex plane, and such 
that the outputs tend to zero as t+ ~, in the presence 
of constant unknown disturbances. Ina later paper, 
Davison [9] has extended the solution of the above problem 
to also take into account a class of unmeasurable arbit- 
rary disturbances that satisfy a differential equation. 
In a more recent paper Willems [10] has derived a set of 
criteria such that the effect on the output of unmeasur- 
able arbitrary disturbances can be reduced as much as 


desired by means of suitable stabilizing state feedback. 


5.3 Statement of the Problem 


Consider the following linear time-invariant mult- 


ivariable controllable system 


Ko = Axte he Bue DE 
(eb; 
y = Cx 
a oe q 
whetce) s.e Rk - iis: the state, u'e.R is the input, € € RK 
is the disturbance and y € Re 7S Che CULDUL aay eo cue 
D are constant matrices of appropriate dimensions. For 


thessystem in veqn, (5.1) let the control objective be to 
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design or construct a state feedback matrix K, to generate 


a tcent rol: 
lioe=" koe (5.2) 


so that (i) a specified performance index or an objective 
functional is minimised and (ii) the resulting closed-loop 
system is stable. It is well known that of the r x n 
degrees of freedom available in the design of a state 
feedback matrix for a controllable system, only n elements, 
or degrees of freedom are required to assign all n eigen- 
values of the system to the desired locations. This leaves 
(cr x n - n) free elements or degrees of freedom in K that 
can be utilized towards minimising a specified performance 
index. For the regulatory control problem a widely 
applicable and useful form of a performance index is a 
linear combination of a quadratic performance index and 

an index or measure of coupling between the disturbances 
and the outputs. An obvious index or measure of coupling 
between a disturbance, Bae and output, am is the mag- 
nitude of the transmittance between bs and ee There- 


fore in its general form the objective functional to be 


minimised can be represented by: 
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where J, is a measure of the coupling between the dis- 


turbances and outputs, Jo is a quadratic performance 


criterion and a oe are weighting factors. J. is defined 


L’ L 
as: 
J. = trace ( 5 6. 1Ge(jo.)- |) (5.4) 
i f=]. °2i' GL a G, (J4;) . 
where || denotes absolute value of elements of GG. the 
index i=l, ..., &£ denotes the range of frequencies Wr Wor 


-++/W, OVE which the objective functional is to be min- 
imised; Bs is the weighting factor to take into account 
the contribution of each measure of disturbance/output 
transmittance or coupling at a particular frequency; and 
G,(j4;) is the weighted transmittance at frequency ws 
between the disturbances € and outputs y, and in its 


general form is defined as: 
, = : Be =] 
G, (ju,) PC(jws1 (A + BR)) DM. (S7e5) 


where P. and M; are m x m and q x q weighting matrices 
for the outputs and the disturbances, respectively. By 
assuming a known dynamic model of the disturbances, 
Muller and Luckel [22] have also considered the design 

of state feedback controllers that minimise the effect of 
external disturbances on the state variables of the 
system. They have quantified the influence of disturb- 


ances by minimising a measure of modal disturbability - 
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that is the inner product of the reciprocal eigenvectors 
of the closed-loop system and the disturbance vector, d; of 


interest.* The well known quadratic performance index, 


a = | (xtox + u Ru) dt (5.6) 
O 


where Q is a symmetric, positive, semidefinite matrix, 
and R is a symmetric, positive, definite matrix. A remark 
regarding the form of the system defined by eqns. (5.1) 


and (5.2) and the form of J. and ae defined by eqns. (5.4) 


L 


and H5e6) iS. in order there. 


Remarks 5.1 + The control law as defined by eqn. (5.2) 
and the corresponding performance index, J, as defined by 
eqn. (5.3) are in no way limited to a closed-loop system 
with only proportional state feedback. Dynamic (integral) 
feedback compensation is easily considered by simply aug- 
menting the original system of eqn. (5.1) to include 
additional state variables representing, for example, 
integrals of the outputs. 

Note also that in general an objective functional 
such as the one defined by eqn. (5.4) suffers from the 
following disadvantage: For finite values of Jy, to exist, 


stability conditions need to be incorporated into the 


* Note: a zero row in Vd; implies an undisturbable mode. 
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EneeTunctlonalsdefined by eqn. 9(5 53). 9 Thais fis: the main 


subject’ of discussion in the following section. 


5d Stability Considerations 


To incorporate stability constraints into the 
system of eqn. (5.1) and (5.2) use is made of an important 
property of the transformation of a system to a phase- 
variance canonical form, that is the invariance of the 
transformation under state feedback. The ensuing analysis 
briefly follows the approach adopted by Gourishankar 
and Ramar [11,12] in their minimisation procedure for 
eigenvalue sensitivity, and also, indirectly, the eigen- 
value assignment procedure of Topaloglu and Seborg [13]. 

By considering state feedback of the form in 


edn. (Sta), butesonePthatvonly utilizes? (r°-" i), control 


inputs, i.e. a control law of the form u' = K'x where u' 
2 ib ees 
[u,, Uy se Be ar era Ae uJ and K [ky k, Ae 
de 
ee iL Ket ne kL] , the system of eqn. (5.1) can be put 


into the following single input form: 
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where b, is the oe column of By» and Ky is the 
oe row Of K. If now the single input system as 
characterized by the pair ye) is controllable and 
feedback via ve is admitted, then by a choice of suitable 
Ree the resulting system can be assigned desired eigen- 
values. It is known that almost any feedback matrix K' 
will make the system of eqn. (5.7) controllable with 


“respect to input us (Davison and Wang [14]). The result- 


ing closed-loop system matrix, H, is then given by: 


t T al Wi 
He=4 (Awe ebs ke tbe k. = (E + b.k. ie 
i=] 11 a aN) pec: 
iF au 
Now denote by matrix T(K') (i.e. a matrix T which is a 


£unction of K') the nonsingular transformation that trans-— 
forms the matrix pair ay au) to a phase-variable canon- 
ical form. Wilkie and Perkins [15] have shown that such 
a transformation matrix is invariant under state feedback 
In the present investigation this fact was found to be 
useful for eigenvalue assignment in the single input 
system Vie) Note that there are many other ways of 
assigning eigenvalues for the pair Ea ce (for example 
Via veqns =(259) jp ,and thensevaluating H. However; in this 
case the procedure via the phase-variable canonical form 
was found to be convenient, with H given by the following 
equation: 


H = T(K')HG(K') 


Where for a desired set of closed-loop eigenvalues H is 
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the phase-variable form of H, and G(K') = mae With 
H defined by eqn. (5.9), J, can be evaluated from: 
J. = trace ( , Be {Pec CG WI=Hj2DM.) tp Uae )t-H) -DM,) |) 
L Lei al i 1 a a ~ 
(510) 
Likewise ue can also be evaluated as a function of K', 
Ramar and Gourishankar [12]. With the control law of the 
form: 
u' = K'x (S21) 
the quadratic performance index, J {K*) is given by: 
Pee = (xtox + u' Ru') dt = xT (Q + K' RK") xdt 
fe) O 
(S22) 
Hor initial conditions x(0) = Xor Ue can be expressed 
as 
5 (K') = xt WCC (3:43) 
q fe) fo) 


where V(K') is the symmetric positive-definite solution 


of the matrix Lyapunov equation 


ak 


Vik MH ee oe eK) RK" So (5.14) 


A criterion which is not dependent on the initial condi- 


tions 48 and is 'optimal' in the average sense can be 


a — — 2 7! nan ; oi 


by + PAD 
dorw i ead 
| otnsliva oes mip on 
arent rie ae ral a 
tt a* * tne eee a |, ar Ja 
are 


ot, ©) | : 7 1 
_ \w. te i kis ed cst ao 5 wort | | 


uy 40 Wh ream nei diet kl) 38 


a 


a, 
_ < A 
toh : x'A = hd 


o 


wears od sippy (178) gill apse, (0) ecole RD 
rete ‘ * 


| , | fe 
(Er .@) a, RV “x = ; 
”, 


destulor satntt 


considered by modifying (5.13) to [25]: 
Liga = trace[V(K') ] (57.1.5) 


Thus eS can be evaluated for different values of K' 


by wsclVving for Vik" )etineean,. (5.14). 


o.4,° Thedesign algorithm 


Thus, very simply, the main step in the Been 
procedure for minimising J requires J, and oe to be’ com- 
puted for a given set of desired eigenvalues and a 
choice of gain parameters in K' as determined by a numer- 
ical search procedure. In detail the algorithm consists 
of the following steps: 

i) An initial choice of K' is made. 


ii) T(K') and G(K') are evaluated for the matrix pair 


(E, a 


iii) With the desired closed-loop eigenvalues specified, 


the closed-loop system characteristic polynomial 


and hence H are evaluated. H is then computed 


£Eromveqn. §(5.9))% 


iv) J. is computed from eqn. (5.10) and ths is computed 


L 
from eqn. (5.04) sand (5.15). 


v) A numerical search procedure is used to estimate 
a new set of values for K' and steps (ii) to 


(iv) repeated until a Ko is found that gives a 


pt 


170 


fe2.0) — 


') Xo sovfay 308 | 
‘abe ooo nt (a 
ei v4 e ne 


i) 


spiaeh odt mt quis akan acts »yaqule ae veut : ~ 
mi oth Gd ot Bitte, oply set itipsy t ie kn bene 8 a <i 03g 
- pre aevlewiepie) BertesS to Jee aevie s sot oa 

15min & Yt  benimeodet mae "8 ai ee et bath 
nse anao maetropis git) yoniet 7, ore ee 
| pres ae 

sdtpuled “A So so sera ets Lent am) an 

£434 ml ga eid’ tom ‘oye oan (sD baw eer 
fartisveda sor Lovgnegas' Aaddimicis botieob vee 
iq-emony Lod ene glia oo At 
basting creel sain exe H # 


- 


> 


pierre Gt 
a : 


hy fal & 


minimum value of J = J, cu ee 


vi) k, is computed from the final T (Kept) and the 
desired set of closed-loop eigenvalues to yield 


a final K 
OvuE. 


Note that unless J = 0 the final aes that is sel- 


ected is not truly optimal, since J is also a function of 
oe and the measure of controllability (eg. cf. Simon and 
Mitter [26]). Thus for the true optimal solution pee eros 
cedure is repeated for j} = 1, ... r (assuming that (E, o-) 
is always controllable) and the final Seay is selected to 
be that corresponding to the minimum of J(bs)y qe Le ce 
It is easy to extend the minimisation procedure 
proposed above to the design of a unity rank feedback 
matrix of the type: K = gti, where g is an r x 1 vector 
and f is an n x 1 vector. In this case the minimisation 
procedure selects parameters in the vector g to minimise 
the performance index defined in eqn. (5.3), while 
the vector f is used to assign desired eigenvalues. 
Such a procedure has clear advantages over conventional 
pole-placement schemes where the choice of g depends 
only on the controllability of the pair (A, Bg) and is 
otherwise arbitrary. In addition this design procedure 
can be used as an alternative when the minimisation 


procedure for the design of nonunity rank feedback matrix 
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cannot be easily applied to high order systems because of 
computational difficulties usually encountered in high 


dimensional search problems. 


5.4.2 Disturbance localisation vs. disturbance 


minimisation 


The time response of the pos system output, y,(t), 
in the presence of external disturbance, ae with the 


initial state vector, x(0) = 0, is given by: 


t n T T 
Vac ey. = : ey cw vd. See ~ t) FE, (t) dt 


(Sil) 


th row of C} Wp and vp are the ee eigen- 


where cy as the a 
vector and reciprocal eigenvector of H, > is the jae col- 
(Shiibo eae hae By dp is the ae eigenvalue and ae LS ‘the — com- 
ponent of disturbance, €. The necessary and sufficient 
condition for disturbance localisation or invariance of 


y, (t) in the presence of ane) is, chah et. al. (Lol (cr. 


Chapter 3)" 


Ay fi 


CpWev pd. = 0 POR aLL par Be eee ey 


(Sil) 


In contrast to this, disturbance minimisation over a wide 


range of frequencies requires that the expression of eqn. 
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(5.4) be a minimum. For example, UB eel) che sien ga = 90 


implies asymptotic tracking, in the presence of step 


disturbances this is 


. ai T 
ey see ee 


ee 


Ths =—nt) (5 213) 


fe 


Note that frequently for low order systems with r =n it 
is possible to carry out the minimisation scheme analytic- 
ally for s = 0 by using expression (5.18) and thus achieve 


zero offset without integral control. 


5.5 Illustrative Examples 


To demonstrate the use of the proposed design 
procedure, a numerical search procedure of Fletcher [17] 
as based on a quasi- Newton method (and available as 
subroutine ZXMIN in the International Mathematical and 
Seacistical Mibranies, IMol bibraty 1, ehdle1ene sto 7), 
was applied to a second order system with coefficient 


Matrices A, B, D and C given by: 


3 -4 1 0 2 ih 
A= B p= 
2) 2 0 Z 2a. 


(Bas) 
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simplicity, consider only the minimisation of J. that 
is, in eqn. (5.3) choose a = 1 and or = 0. Furthermore, 
assume that ey is the main disturbance variable of 
interest in the above system so that the main objective 
is to regulate Y= X} well in the presence of disturbance 
aT and to a lesser extent Eo: For the purpose of com- 
parison of the present design scheme with other control 


schemes it was decided to limit the absolute value of 


the gain elements in the controller to 15. 


9.9.2 Proportional, control 


The values of parameters (, Woe Bis Pi, M, as 
defined in eqn. (5.4) and the resulting disturbance min- 
imisation controller, Ky are shown in Table 5.1. Also 
shown in Table 5.1 is an optimal controller, Koy designed 
for the system of eqn. (5.19) by minimising a quadratic 
performance index. The choice of state and control 
weighting matrices Q and R respectively is shown in 
Table 5.1. Since ae exists for the second order system, 
it was possible to also design a controller, K, by a 
combination of pole-placement and analytical minimisation 


technigues. Details of controller, K3, are also shown in 


Table. 5 sic. 
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Figures 5.1 and 5.2 show the response of x, (t) 
to a unit step disturbance in Eo (t), respectively, when 
1! Ky and K are 


implemented. In each case, as would be expected, the 


each of the three feedback matrices K 


response of the system with the disturbance minimisa- 
tion controller Ky is better or at least comparable to 
that of the optimal system. Figure 5.3 and 5.4 show 
the response of x5 (t) to a unit step disturbance in 

Ey (t) and BE, (t) respectively. In Figures 5.3 the res- 
ponse of x, (t) to a unit step disturbance in E, ft) is 
worse than that of the system with optimal or pole- 
placement controller. This is intuitively expected 
Since in trying to reduce the effect of Ey (t) on 1, (4) 
as required for the design objective with limited 
gains and therefore limited control energy, the effect 
of the disturbance is necessarily diverted to the 
remaining state variable of the system, which in this 


case is x(t). 


52.2 Proportional pluseintegral control 


By augmenting the original system of eqn. (5.19) 
to a third order system where X, =X) integral feed- 
back Of Y, Or x, can also be considered. A disturbance 
minimisation controller, Kar was designed for the above 


system and is shown in Table 5.2, together with the 
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values of required parameters as defined in eqn. (5.4). 
For the purpose of comparison an optimal PI controller, 
Key designed by minimising a quadratic performance index 
was also calculated and is shown in Table 5.2. Finally 
two different unity rank state feedback matrices Ke and 
Ko were also designed. These are also shown in Table 

a2a0, Controller Ke was designed by using a conventional 


pole placement scheme with an arbitrary choice of 


g*="[0 rf eae where Ke = gtl. K. was designed by applic- 


V 
ation of the minimisation procedure in the optimal 
choice of the g. 

Figure 5.5 and 5.6 show the response of y, (t) 
to unit step disturbances in E, (t) and Eo (t) respectively. 
In each case, aS would be expected, the disturbance 
MinimiSation controller, Kae outperforms the conven- 
tional pole placement controller Ke and is better than 


or at least comparable to the performance of the optimal 


System. 


5.6 Applications to Observer Design 


An immediate use of the above method is in its 
application to the dual or adjoint of system of eqn. 
(5.1) and thus to the design of observers in which the 
effect of ummeasurable disturbances on the state estim- 


ates of interests can be minimised. The method is of 
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course only directly applicable to the design of ident- 
ity, or full-order observers, and not to the design of 
reduced order observers which have necessarily limited 
degrees of freedom. The flexibility offered by the 
method in allowing arbitrary eigenvalue placement is 
especially useful here. 

Johnson has considered the synthesis problem of 
unknown input observers by modelling the unknown inputs 


by a dynamical system or approximating them by a poly- 


nomial function; all of his results have been succintly 


summarized in [18]. 

Basile and Marro [19] have discussed the design 
of full order observers for systems with unknown inputs 
and disturbances. Their method requires a recursive 
algorithm for design, but the procedure for determina- 
tion of stability conditions requires a complex 
algorithm. Mita [20] has recently proposed a method 
for the design of a reduced order observers for multi- 
variable systems such that the required state estimates 
are made insensitive to unknown inputs or disturbances. 

More recently Gourishankar et al. [21] have 
also considered the problem of designing observers for 
systems subjected to unmeasurable inputs that can be 
approximated by a polynomial type function. They 


have also derived a set of necessary and sufficient 
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conditions that can be readily applied to determine the 


existence of an observer for such a case. 


the design of ann 


eqn. 


i) 


a1) 


eis) 


For completeness the main steps involved in 


order observer for the system of 


(5.1) are now summarized: 


If the output matrix C in the system of eqn. (5.1) 
does not satisfy C = [tO] then the system of 

eqn. (5.1) is transformed such that the final 
transformed system has an output matrix C', where 
c' = [IO]. For convenience of notation drop the 
Superscript ' and always assume from here on that 
¢ is) in the mequired ‘ror. 

Let the dynamics of the observer be governed by 


the following equation: 
2 = Fz + Gy + Bu (5.20) 


with F = A - GC, and the required estimate, yee de- 


; A “w~ nw ~ T A 
fined by Y, a [xX aq! Xmt2! coos x] aS 
T T 


Thus the observer output equation is 


y. = Nz, where Nn? is a basis matrix (5.21) 
iL er ee 


The error dynamics of actual states Xt! Xm! 


seer KX and the observer states Zt! 2me2! oct! 2y 
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are given by the transfer function matrix: 


a(e) eect GA = caylee. ete snes! 
where e\(t)iacs ree) ~ Z(t) (xy (t) = Zag lt))e 
rete, (x(t) _ z(t))1. Thus for a system with 


state matrix ibe input matrix cy disturbance mat- 
ae Bt Noe and output matrix oo the method of 
sections: 5.3 and 5.4. can be applied directly £0 
construct ee such that the following measure of 


coupling or transmittance between &(s) and e(s) 


is minimised: 


ee 
ee iy 4B = 
Tro = trace (2 81 (Po (sI - (A = ohne Bs = néM.) 


i=l 


=2 


B,D (sI - (AT-C'G"))~~ N™M,) | (5.23) 


9.6.1 Example 


Consider the following state space matrix A, 
the output matrix C and the disturbance matrix D of a 
third order system for which an observer to estimate 


x, (t) is to be designed. 
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£ 
and D = 0 (5.24) 
2 


If the disturbance 4 is unmeasurable, then it is clear 
from the error dynamics of the full order observer (cf. 
eqn. (5.22)) that the error in the estimate of x, (t) 
has a non-zero value for constant disturbances. Thus, 
let the main objective in the design of an observer for 
the system of eqn. (5.24) be to obtain a good estimate 
of x, (t) even in the presence of unknown disturbances 
in &,(t). 

By considering the dual of the above system 
where the state matrix is Ae input matrix is Coe 
output matrix is De and disturbance matrix is [1 0 ghee 
the proposed procedure for disturbance minimisation can 
be applied to determine an optimum observer gain matrix, 
G. Three different observer gain matrices, Gy G. and 
G3 were designed incorporating different design para- 


meters and specified closed-loop eigenvalues. Values 


G. and G, and the corresponding details of the 


of G), G 3 
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parameter values used in the design scheme are shown 

in Table 5.3. Figure 5.7 shows the error e, (t) for the 
estimate of x, (t), in the presence of a unit step dis- 
turbance in €, (t). The time response of e, (t) on 

Figure 5.7 graphically illustrates the design freedom 
available in the disturbance minimisation scheme, 
namely, arbitrary eigenvalue assignment plus minimisation 
of disturbances at one or more frequencies. To compare 
the design of the observer gain matrix by a disturbance 
Minimisation scheme to other conventional design schemes 
for a full order observer, the optimal control scheme 
was applied to the dual of the system of egn. (5.24). 


The resulting gain matrix G, as well as the weighting 


4 
matrices Q and R are shown in Table 5.3. A reduced 
order observer was also designed and the resulting 


gain matrix G. with other design details is shown in 


5 
Table 5.3. Figure 5.8 shows the error e, (t) in the 
presence of a unit step disturbance in Ey (t) when three 
different gain matrices Gar G, and G. are implemented 
for the system of eqn. (5.24). As would be expected, 
the gain matrix Gor designed by using the disturbance 


minimisation scheme performs better than the two 


other schemes. 
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Brad. Experimental Evaluation 


The design approach developed here to minimise 
the effect of external disturbances on outputs of 
interest was evaluated experimentally by designing a 
controller for the reduced third-order model of a com- 
puter-controlled, pilot-plant evaporator at the Univer- 
sity of Alberta. A schematic diagram of the evaporator, 
with details on its mode of operation, and a table of 
important variables together with its discrete third- 
order state space model, are given in Appendix C. This 
evaporator has been used previously to evaluate a 
number of modern multivariable control methods (Fisher 
and Seborg,[23]) and hence it is possible to make direct 
comparisons between alternative control techniques. 

The third order evaporator model has 3 state variables 
(which are also the 3 output variables), 3 inputs and 
3 disturbances. 

The main output variable of interest is the 
second-effect or product concentration, C2. Previous 
operating experience has shown that the most frequent 
and severe disturbances in product concentration, C2, 
are produced by variations in feed flowrate, F.. There- 
fore for the purpose of this application the design 


objectives) for the controller are: 


snags te babe snes 4 peouhes ets 08 el. 
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_ 


a) to minimise the effect of feed disturbance, F, 
Onw che =producteconcentration, €2, and 
b) to assign all closed-loop eigenvalues closer to 


the origin. 


A simple inspection of matrices $6,A (for the 
third order model as tabulated in Appendix C) reveals 
that the input variable, B2, can only control or man- 
ipulate the state variable, W2, i.e. C2 and Wl are 
structurally uncontrollable by B2: For this neason B2 
was first used to assign the eigenvalue of the mode 
corresponding to W2. The resulting ae order system, 
with 2 inputs, 3 disturbances and B2 = 13.05*wW2 is 


shown below: 


0 0 Egeoszse 00ers. 
or Fe Ondyeeton te im 32 Ego s77 0.0854 
0 Mo f96 00507 §20.0441 

Ovr20 0 20 toss 

Puen | Wai ase 0 -0.0156 


UO Gomer Oyo S00 0.0218 


The above step serves two main purposes. First, it 
changes the analytically simple control problem (be- 
cause A can be inverted) into a not so straightforward 


problem. This is useful in demonstrating the practicality 
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of the method. Secondly, it reduces the dimension- 
ality of the numerical search procedure from 6 to 3, 
which is computationally easier to handle. With the 
twowcontrol) objectives: way “andi, b)) in mma, ne 
proposed design procedure was used to minimise the 
effect of disturbance , F, on C2, at steady state, i.e. 
at w = 0. The desired eigenvalues of the closed-loop 
system were chosen as 0.65, 0.47 and 0.28. The result- 


ing Gontroller, “K is as shown in Table 5.4. Also 


dmin’ 


shown in Table 5.4 is controller, K designed by 


ope! 
using optimal control method, Wilson [24]. The values 
chosen for matrices, Q and R in computing Rope are 
also shown in Table 5.4. 

The response of the closed-loop system to a 
+ 20% disturbance in feed flow rate with controller 5 
Ko © and K implemented is shown in Figures 5.9 and 
dmin opt 
5.10 respectively. As expected, the controller Se ae 
reduces the steady state offset in C2 to almost zero. 
The performance of this controller is better than or 
at least as good as the performance of the optimal 


Gontrovler LS. Operating experience on the evaporator 


pt’ 
has shown that disturbances in feed concentation, CF, 
does not significantly affect the product concentration, 


C2 ra To rcheck 1£ minimising, the efiecteor F on C2ihad 


not made C2 more sensitive to disturbances in CF, a 
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+ 30% disturbance was introduced in the feed concentra- 
tion, CF, with controller Bae implemented. Figure 
5.11 shows the results of this run. Even when most of 
the control effort (via design) was concentrated in 


minimising the effect of F on C2, the product concentra- 


tion sensitivity to CF had not increased. 


5.8 Conclusions 


i) The suggested procedure allows simultaneous 
eigenvalue assignment and disturbance minimisation 
and has immediate applications to the design of 
observers for systems with unmeasurable disturb- 


ances. 


ii) The objective functional defined by eqn. (5.3) 
is more easily evaluated than in Willems, method 
[10]. The basic computational requirements for 
the design procedure are modest except for the 
numerical search procedure which does not always 
perform well, especially for higher order 


systems. 


iii) The illustrative example shows the advantages 
of the proposed method in comparison to methods 
such as pole assignment by unity-rank state 
feedback which do not take into account the in- 
formation available in the disturbance matrix, D, 


of the process models. 
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iv) An experimental evaluation of the proposed design 
scheme on the pilot-plant double effect evaporator 
gave good results and also demonstrated the 


practicality of the method. 


Some aspects of the proposed procedure are worthy 
of further considerations. For example, it would be 
useful to investigate necessary and sufficient conditions 
for a minimum Jy, = 0 to exist (at certain w # 0) for 
a class of unknown disturbances that satisfy a differ- 
ential equation. Obviously the results in [1-3] are 


sufficient but not necessary for J. to be identically 
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CHAPTER 6 


EIGENVALUE INVARIANCE TO SYSTEM PARAMETER VARIATIONS 


BY EIGENVECTOR ASSIGNMENT 


Abstract: 


A method is presented for designing constant 
multivariable feedback controllers that make selected 
closed-loop eigenvalues tnvartant to unknown perturbations 
Of arbitrary magnitude in system parameters. The suffic- 
ient conditions for eigenvalue invariance are expressed 
in terms of the structure of the closed-loop system 
matrices, and a numerical example is included to show how 
eigenvalue/eigenvector assignment techniques can be used 
to produce the desired eigenvector structure and thereby, 


eigenvalue invariance. 
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6.1 Introduction 


In most industrial applications, perturbations or 
errors in the system parameters and/or feedback gains are 
quite common. These perturbations or errors might be due 
to changes in the plant model, faulty instrumentation, 
etc. If the control system design is based on a specific 
set of assumed model parameters, then perturbations in 
these parameter values can cause deterioration in the 
plant performance, or may even result in an unstable 
system. 

Eigenvalue or pole-placement algorithms have been 
widely used in both single variable and multivariable 
design applications because they give the user the ability 
to specify system stability and general characteristics 
of the time domain response. Hence it would be desirable 
when using pole placement (or other modern design techni- 
ques that rely on system models) to ensure, as much as 
possible, that perturbations in system parameters do not 
negate the original design objectives. 

The purpose of this chapter is to show how, and 
under what conditions, selected closed-loop eigenvalues 
can be made invariant to perturbations in the system 
parameters. The question of how many, and which eigenval- 
ues can be made invariant to perturbations is also con- 


Sidered. 
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This chapter is organised as follows. Section 

6.2 discusses relevant previous work carried out in the 
area of eigenvalue insensitivity, Section 6.3 formulates 
the control problem and Section 6.4 is concerned with 
derivation of results to achieve eigenvalue insensitivity 
to small perturbations. In Sections 6.5 and 6.6 we con- 
sider the more important result, namely achieving eigen- 
value invariance to arbitrary perturbations by eigen- 
vector assignment, Section 6.7 deals with a numerical 
example to illustrate the design paper and is followed by 


Section 6.8 on Conclusions. 


6.2 Previous Work 


The problem of sensitivity reduction of closed- 
loop eigenvalues has been previously considered by 
Tzafestas and Paraskevopoulos [1], who designed a state 
feedback controller to assign closed-loop eigenvalues 
and to reduce eigenvalue sensitivity to parameter varia- 
tions. More recently, Mita and Ngamkajornvivat [2] have 
outlined a method for single-input systems that achieve 
eigenvalue invariance to a class of unity-rank perturba- 
tions in the system matrix. Seraji [3] has proposed a 
method in which unity rank state feedback control is 
employed to assign closed-loop eigenvalues and also make 


the dominant eigenvalues insensitive to given variations 
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in the system parameters. This procedure is not suitable 
for programming on a computer, and also leads to various 
difficulties when applied to high order systems. Gouris- 
hankar and Ramar [4] have outlined methods for the design 
of unity and non-unity rank state feedback controllers 
which assign closed-loop eigenvalues and also minimise a 
measure Of eigenvalue sensitivity to given parameter 
variations. Karlin et al. [5] have considered the 
synthesis of feedback control laws which make linear 
system state or output trajectories insensitive to small 
parameter variations. The existence of such control laws 


necessarily requires restrictive conditions. 


6.3 Statement of the Problem 


Consider the following nie order, completely 
controllable, linear multivariable system with r inputs, 


ult), GQ disturbances et). and m Outputs, y(t): 


x(t) = Ax(t) + Bu(t) + D&(t) 
(62.1) 


CC) 


y (t) 


and the state feedback control law: 


u(t) = Kx(t) (6.2) 
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The resulting closed-loop system is described by: 

x(t) = Hx(t) + DE (t) (6.3) 
where H = A + BK. 
In terms of individual matrix elements: 


= te ame: AG (6.4) 


where MMe and Sg represent the Cee elements of H and 


A respectively; b; represents the yee 


h 


COW, Ol eb and 


represents the 4 column of K. Let the perturbations 


Or variations in H, due to variations in A or K, be denoted 


by SH. In this investigation the main questions of 
interest are: 


i) Under what conditions can the closed-loop 


eigenvalues be made invariant to perturbations 


6H? 
ii) How many and which eigenvalues can be made 


invariant? 


iii) What design procedure will produce eigenvalue 


invariance using only proportional feedback 


CONnELOL: 


6.3.1 Practical Significance 


The case where the perturbations 6H must be known, 


constant values is of little practical significance 
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because it would then usually be possible to carry out 
the design of the control system based on the corrected 
system matrix (H + 6H). However, the case where the 
perturbations are unknown is of considerable importance. 
For example, modelling errors and/or real-time variations 


th 


can produce errors in the i column of A, which directly 


ane ; : ‘ 
affect the. i h column of H. Variations in the transducer 


gains in the jee loop, and any other perturbations that 


affect only the ade column of K, can also be handled dir- 
ectly since they contribute only to the yee column, of Hh. 


thi 


Thus variations in the i column of H can be represented 


as: 


“il 
G. Se Arab (Se Symes 0 
Nice recs SA 0 
1.663 + b9(6ky) --. 
ee 
Oe. Cae pats ann 


where by represents the een row of B. 


Note that for the general case, variations in 
elements of B cannot be handled by the procedures pre- 
sented in this paper. However, Ramar and Gourishankar 
[6] have recently proposed a method for the design of 
unrestricted rank state feedback that minimises the sen- 
sitivity of closed-loop poles to variations in matrices 


Bor K. For special cases where the number of inputs is 
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equal to n then complete eigenvalue invariance to perturb- 
ations in B can be achieved by considering the dual or 
the adjoint of the system in eqn. (6.1) (Chapter 2, see 


auso, onan et -al. = iv ii. 


6.4 Eigenvalue Insensitivity to Small Perturbations 


Faddeeva and Faddeev [8] have shown that the fol- 
lowing first-order relationship can be used to approximate 
the change in a distinct eigenvalue, hes due to small 


perturbations in the elements of the system matrix H: 


meet 
5A, = v, (SH) Wa (6)=6)) 


th 


where Wi and v, are the i eigenvector and reciprocal 


eigenvector of H respectively. Therefore, for a given 


eh eigenvalue will be invariant to small 


system the i 
perturbations, if the right hand side of eqn. (6.6) is 
identically zero. The four conditions under which this 


Can occur are: 


a) (6H) w,; = 0 
b) v, (SH) = 0 


(AT) 


CG} (SH) Wi ow. where i #jandaisa 


non-zero scalar 


d) oe ashe) bv; where i # j and 8 is a 


non-zero scalar . 


; vO 8 
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6.4/1. Eigenvalue Insensitivity by Eigenvector Assignment 


The approach taken first in the following analysis 
is to examine the system eigenvectors and determine under 
what conditions eqn. (6.7a) is satisfied. 

Obviously, if 


{ We } C Ker{éH} 1s 1 ae ey Oe Oe ey) 


where p = dimension Ker{éH} then eqn. (6.7a) is satisfied 
and at least the p closed-loop eigenvalues corresponding 
to the eigenvectors {was i=l, ...,p} are insensitive to 
any small perturbations, 6H. 

One approach to the design problem would be to 
specify the closed-loop eigenvectors such that they satisfy 
condition (6.8) directly. Unfortunately when constant 
State feedback control is used to assign closed-loop 
eigenvalues and eigenvectors, then only r elements (where 
r is the number of inputs) of each of the n eigenvectors 
can be assigned arbitrarily, and this is possible if and 


Omari is 


i) the resulting n eigenvectors are linearly 
independent, and 
ii) the corresponding r rows of the input 
matrix B are linearly independent. 
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Because (n-r) elements of each eigenvector cannot be 
assigned arbitrarily, and because it is important to 
distinguish between arbitrary and known system perturba- 
tions, it is convenient to reorder the columns of 6H 


(if necessary) and partition 6H such that: 
} CaS) 


where OH, is an nxk matrix containing all the unknown 
perturbations (ki < 4x) 
6H, is an nx(r-k) matrix of known constant per- 
turbations 
SH is an nx(n-r) matrix of known constant per- 
turbations. 

Then, if the Moat eigenvector is partitioned into an rxl 


vector Wl, and a (n-r)xl vector, W2., eqn. (6.7a) can be 


rewritten as: 


! t Wl, 
eye od ay Wo. = 0 (Ga0} 
1 


Since the r elements of Wl, can be assigned arbitrarily, 


it is usually possible (subject to the two constraints 


noted above) to satisfy the condition: 
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However, the second part of the sufficiency condition in 


egn. (6-10) requires that: 


6H, W2, uO mr GO tetas) pn ane 1 (Ge 2) 


and W2, cannot be specified arbitrarily. Sufficient 


conditions for eqn. i(6.12)) to hold are: 


6H. A, = 0 

SH, Ae 0 (65.3) 
and 

6H Bo = 0 


where A A, and B. are partitions of A and B as shown 


oer | 2 
in the following partitioned form of H = A + BK = WJV: 


! ' : ) t a | 
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(6.14) 

1 ga (Glare) ea Qelsie iy 

where Ay: Bie Ky and Hi€& R : A, and H,€ R ; 
Bae Re a matrix; Kie€ prx(n-r) , ancdewW, a ences Wee 
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Broo; FOr Susy pietencey 


Ther sulfliciene, conditions for eqn. (6.12) wo 
hold, as given by eqn. (6.13), can be derived starting 
with the following two equations which are a consequence 


Otvgeqn. (6. 14): 


Wade os 3 ky 
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Assuming eqn. (6.13) to be true, then pre- 


mMuLtiplying both of eqns. (6.15) by oH post-multiplying 


3/ 


the first equation in (6.15) by W, and the second by W 


i) Be 
and then adding the two equations yields: 
SH,W,J, (VW, + VW) + SH4W,J, (VW, - VW) = 0 
(6.216) 
From the identity VW = I it is known that 
ViWi + VW. = I. and VW, + VAW3 = 0 and hence, know- 
ing that Jy Z# 0, egn. (6.16) implies 
6HW, = 0 (6,17) 


Likewise pre-multiplying both equations in (6.15) by 


6H post-multiplying the first equation by W. and the 


Bu 
second by Wai and then adding the two equations yields: 
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429 + VAWy) = 0 


SHW.J, (VW, + VW) + SHW J (VW. 


(Go) 


By the use of the identity VW = I, eqn. (6.18) implies: 


SH4W, = 0 (6 39) 


Since, 


[W2, W2, Sanpete W2, A Oren W2 J = [W. WwW] 


eqns. (6.17)*and (6.19) snow that eqn. (6.12) holds wt 
the conditions expressed by eqn. (6.13) are correct. 
This proves sufficiency. 

The sufficient conditions for eigenvalue invar- 
iance expressed in eqn. (6.11) and (6.13) are obviously 
very restrictive but it should be noted that they are 
satisfied for the special case when 6H = 0, i.e. when 
there are no variations in at least (n-r) columns of H. 
In this case at least p eigenvalues can be made invariant 


where; 
p = n-rank { SH, } - rank {6H} (6.20) 


Normally, since all k columns of 6H, contain arbitrary 
elements, it can be assumed that rank (oH, } =a eeand 
that its columns are linearly independent with those of 


6H. 
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The above development of conditions governing 
eigenvalue invariance to small perturbations in H was 
based on eqn. (6.7a). However, it can easily be shown 
thatt iteeqn.. (6. /a)eas satisfied then the use of eqn. 
(6.7b) will not lead to any additional invariant eigen- 
values. This follows since eqn. (6.8) implies that rank 
(SH) = n-p and that the eigenvectors {was Dawid eegecy. (TiS 
span 6H. Therefore v, (6H) 40 fOr 26= (Oey ely eee 
LE COnNCLeE1OnN. (6./a) Usesatistied, condition (6.70) cCan— 
not be met for i>p. Also it was not possible to develop 
an improved design procedure based on eqns. (6.7c) and 
(6.7d) and hence they are not considered further in this 
paper. 

The results of the above discussion can be sum- 
marized as follows. Sufficient conditions for (n-k) 
eigenvalues of the closed-loop system to remain insen- 
Sitive to unknown perturbations in elements of H are: 

i) that the perturbations be small, 
ii) that the perturbations occur in only the 
first k columns where k is less than or 
equal to the number of inputs, r, 
iii) that the elements in the remaining (n-k) 
columns of H be known constant values, 
iv) the first k elements of each of the (n-k) 
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In the following sections, sufficient conditions 
under which eigenvalues can be made invariant to unknown 
perturbations, 6H, of arbitrarily large magnitude are 
developed and a design procedure is presented that will 


produce a closed-loop system that meets these conditions. 


6.5 Eigenvalue Invariance to Arbitrary Perturbations 


The following sufficient conditions for eigenvalue 
invariance are derived based on the properties of quasi- 
triangular matrices rather than starting with eqn. (6.6) 


which is only a first order approximation. 


THEOREM 6. 1 : 


A sufftietent conditton for (n-k) of the etgen- 
values of the closed-loop system matrix, H, to rematn 
tnvartant to arbttrary and unknown perturbattons tn the 
frrerarn columne Of H, ts enart Hf pe Of; tne Fovlouing j7orm, 
OF Thar te can be Prougne tnvo the fobLLowing form py 


reordering of the state variables: 
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ee (aut) 
i ee 
where H, € peak He pin-kiatn-k) the (n-k) tnvartant 


eigenvalues {Ay yze hee? eel cas ai are the etgenvalues 


7 grok kine: 
natin th OF oe 

ote shun. i: 

Liiw dada Beoneas yf al ' soord dp 


eo ares — mate ya gost 


sro bincdtyusi 19 


giitevnapha 1ot Zn* a nbnIo9 peranere eubwotios set, ils mo 
‘go 2eQene oy ite hoasd bovineb ain wart secs 
daly cored sets” yots ot asolexsam aelupn noi 


=heGup I¢ 


i 
1 wi - Ser ° 


; ine aamesconnen p «sb x0 reek s ying ak 7% 
sane ie ae 


a 


(Ro a} (hie 


a . a 


Fil 


es ) , ne phen 
ea whi ee wt heh). ene merasbAgo “enheera 


mit ~ 1 


Meanie Of oF een natant Ge" yi -peaaia wht Yo 
ihe mer 7 nie: ven enn ako 


siny * . ve 

; a 

err: "v ary Vy Gila ay Nhe 2%," ‘ye fata a “4 on vernon 8 aod : ae 
ail whey aati otal aie ee iiewone eA De ay . 

. 2 geben eo eres” att Arig at 


= 


re - n j 9 
4 sa a yy, i= “8 
pinesie ft 


-_ 


- a ie 
; 


i 
’ a) ay 


214 


OF Hy and all perturbattons are assumed to occur in H, 


and/or H 2. 


Proof or Theorem 6.1: 


Ef H is ‘quasi-triangular as defined by eqn. (6.21), 


then the n eigenvalues of H are the k eigenvalues of Hy 


and the (n-k) eigenvalues of H By assumption the per- 


4° 
1 and/or H and hence Hye and 


its eigenvalues are invariant. Alternatively, from eqn. 


turbations occur only in H 


(O8L4)) it. follows ‘that Hac] Ac+BoK Stand: His 2 +BOK, and 


2 Zea mle 4 4 


L’ A and K,- 


hence are invariant to perturbations in A 


Remark 6.1: ite suriicrenticondi tons moldege1Omrk sn 


butr-design using »eigenvector assignment requires k<r: 


Remark 6.2: ise H. = 0 then the eigenvalues of H are 
the eigenvalues of Hy plus Hyr which are both independent 


of A Therefore all of the n closed-loop eigenvalues 


3° 
are invariant to perturbations in A. 

Remark (G0 © The time domain trajectories of the state 
variables {x,, L = kel. ent due to initial, condition 
or disturbance excitation are invariant to perturbations 
in the first k columns of H if the theorem is satisfied 
and the initial conditions and the disturbance matrix, D, 


satisfy the conditions: 
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[H,: H,]x(0) = 0) 
[H,H,]D = 0 (6-22) 
where Hz =-0 


Cruz and Perkins [9] have also derived a similar result. 
However, their result is not constructive, it is essent- 


Vally sa structural. result. 


Remark 6.4: TE Ene VOuUEpUt Matrix C as, of therronm 

C= bees 0] then the system modes associated with the 
invariant eigenvalues are unobservable. Furthermore, 
for such a C, the invariant zeros defined and discussed 
by Bengtsson [10] are the same as some (or all) of the 
invariant eigenvalues as defined by Theorem 6.1. The 
relationship between invariant zeros and eigenvalue 


invariance is explored in detail in Appendix D. 


Remark 6.5: An H matrix which satisfies the conditions 
of Theorem 6.1, also satisfies the sufficient conditions 
for insensitivity to small perturbations discussed in 
Section 6.2. However, Theorem 6.1 provides sufficient 
conditions for invariance to perturbations of arbitrary 
magnitude. 

As shown in the next section and by the numerical 
example, it is possible to design state feedback control- 


lers that produce a closed-loop system matrix, H, with the 
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structure defined in eqn. (6.21), and hence invariant 


eigenvalues. 


6.6 Design Procedure for Eigenvalue Invariance 


It has been shown in Chapter 3 that if H, = 0 
then Wo = 0 and vice-versa (see also Shah et al. [1l1]). 
Therefore an alternative condition to that expressed 


in Theorem 6.1 is that W have the structure: 


W, 10 
W = [e=-1--- (6.23) 
Wala 
Thus any eigenvector assignment procedure that sets Wo = 0 


will result in a system with (n-k) eigenvalues that are 
invariant. tojperturbations insthe first Kwcolumns of Hh. 
The algorithm developed by Srinathkumar and Rhoten [12], 
discussed in Chapters 2 and 3, permits arbitrary assign- 
ment of up to r elements of each eigenvector plus all of 


the system eigenvalues and thus can be used to achieve 


eigenvalue placement and invariance. The main restrictions 


are: 


i) The eigenvalues must be distinct; V = wt 


must exist; and By (cf eqn. 6.14) must be 


non-singular as noted earlier. 
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ii)’ If k,;, the number of columns of H “containing 
arbitrary and unknown perturbations, is less 
than the number of inputs, r, then Wo can be 
set to zero and all eigenvalues can be 
assigned arbitrarily. However, if k =r 
then (n-r) eigenvalues of the closed-loop 
system are the eigenvalues of [A, ~ Byer eas! 
and cannot be specified arbitrarily (cf. Chap- 


ter 3). 


6.7 A Numerical Example 


A linear, third-order system with two inputs 
(Gourishankar and Ramar, [4]) is used to illustrate the 
design procedure for eigenvalue invariance. The system 
is open-loop, unstable with eigenvalues at 0, 1 and 1 and 


COctticilent matrices (cf. can. 6-1): 


Oe ene) i 0 
Ace | 07 a Bes Oe (6-24) 
BP MOL ali ay peal 


For the purpose of this paper it is assumed that there 
are errors in the transducer gains associated with X5 
plus perturbations in’ the second column of A. (The design 


Objective 1S to place the cigenvalues at*-1, -2 and —3 
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respectively and also to make the dominant eigenvalues 
invariant to the errors associated with X5 and to the 
parameter variations in A. 

The first step is to reorder the state vector 
such that xt = [x5, X10 x3] and hence that all of the 
arbitrary and unknown variations are in the first column 
of H. (The Superscript ~ denotes the reordered system). 
Then eqn. (6-20) indicates that at least p = 3-1-0 = 2 
eigenvalues can be made invariant. The design procedure 
for eigenvalue/eigenvector assignment using constant 


state feedback (Srinathkumar and Rhoten [12]) can then be 


used to assign the eigenvalues to -3, -2, -1l and the 


first two elements of each eigenvector as (cf. eqn. (6.14)). 


W, = W, = (6-25) 


The resulting closed-loop eigenvector and system matrices 
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The corresponding feedback control matrix calculated 
from the relationship H = A + BK, using the pseudoinverse 
Of B Ys denoted by Wl and 1s. shown in Table 6.1. “Rt is 
obvious that W, and hence H, in eqn. (6.26) have the 
structure specified by Theorem 6.1 with k=l. Hence, 
Theorem 6.1 shows that the (n-k) eigenvalues -2 and -1l 
will be invariant to arbitrary perturbations in the first 
k=1 columns of H. (Note the following discussion is in 
terms of the original rather than the reordered system.) 
Table 6.1 includes the numerical values for 


BAG one i.e. the sensitivity of the 7 


eigenvalue to 
Vabtations in the ae element of the second column of H. 
These values confirm that eigenvalue/eigenvector assignment 
design procedure did make the two dominant eigenvalues, 
=eaind= —2, 1nVat lant seo saboltbaryevablattonus ein snes second 
column of H. For comparison purposes three other feed- 
back matrices that position the closed-loop eigenvalues 

at -l, -2 and -3 were also considered. The results of 
Gourishankar and Ramar [4], who employed a numerical min- 
imisation approach to reduce eigenvalue sensitivity to 
variations in hy> and ho» for the same example, are in- 
cluded in Table 6.1 where K2 denotes their state feedback 
Matrix. Table 1 also includes results for two unity- 

rank feedback matrices, K3 and K4, designed using conven- 


tional pole assignment techniques. 


The sensitivity of the eigenvalues of the closed- 
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loop system using K2 as the state feedback controller 
(Gourishankar and Ramar [4]) is small as would be expected 
from using a minimisation procedure. Note, however, that 
the elements of K2 are generally larger than those of Kl. 
The proposed eigenvector design method is relatively 
straightforward and requires much less computational ef- 
fort than does the minimisation procedure. Furthermore, 
the proposed method avoids the convergence and ite he 
ality problems encountered in minimisation procedures 
such as the one used in [4]. 

It is well known that when a state feedback matrix, 
K, is constrained to have the unity rank structure, 
K = gf!, the choice of the rxl vector, g, is arbitrary as 
long as the single-input system (A, Bg) is controllable. 
In, computing K3..the choice, g = [1 Tha was made. For this 
controller the sensitivity of the eigenvalues to variations 
in the second column of H was found to be high (see row 3 
of Table 6-1). Different values of g were then selected 
by ‘trial and error and it.was found that, for ig = i Cine 
the resulting controller, K4, had reduced eigenvalue sen- 
sitivity to variations in the second column of H. The 
point to note here is that most conventional pole place- 
ment design techniques give little guidance as to how g 
should be selected so as to reduce eigenvalue sensitivity. 

The significance of the numerical values of the 


eigenvalue sensitivity coefficients is illustrated by the 
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Simulated responses in Fig. 6.1. Here, the closed-loop 
response of x, to the itil CONG. tLon, x5 (0) = 1.0 is 
shown for each of the four controllers (solid line). 

Also included in the figure is the closed-loop response 


when a 20% error in x. exists due to varratrons in. the 


2 
transducer gain (dashed curves). The figure clearly 
illustrates the advantages of using an eigenvalue invar- 
iance, or a minimum eigenvalue sensitivity procedure, 
when designing for pole placement. 

The results are not reproduced here, but the 
time-domain trajectories of xy and X3 (using K,) CO anpelal 


conditions excitation of x, (0) = 1.0 and x, (0) 1.0 were 


completely invariant to the errors in Xo Since the system 
meets the conditions of Theorem 6.1 and Remark 6.3. 

Note that Theorem 6.1 means that even if the con- 
troller design is based on a model that contains errors 
in the first k columns of A (or of H), then n-k eigenvalues 
will still retain their assigned values and be invariant. 
In this example, as implied by Remark 6.2, all of the 
system eigenvalues are completely invariant to arbitrary 


variations in Az, a5. Ci ani A5y and asi: Variations in 


Se 


1 would affect only the eigenvalues of we 
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6.8 Conclusions 


Theorem 6-1 provides sufficient conditions for 
invariance of (n-k) eigenvalues of a closed-loop system 
to arbitrary and unknown perturbations in k columns of 
the closed-loop system matrix, H. By contrast, most 
previous work is based on first-order approximations to 
the variations in closed-loop eigenvalues and provides 
only tnsensittvity to small and/or known perturbations. 
However, the important result from a design point of 
view is that in many practical applications, selected 
eigenvalues can be made invariant to factors such as 
poor parameter estimates, modelling errors, variations 
in transducer gains, etc. The illustrative example 
shows the advantages of the eigenvalue/eigenvector assign- 
ment technique in comparison to other methods for reduc- 
ing the effect of parameter variations on selected closed- 


loop eigenvalues. 
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FIGURE 6.1: 


no measurement error 


20% error in transducer 
gain for x, 


TIME 


Simulated closed-loop responses of the 
state space model defined by equation 
(6.24) and the four proportional feed- 
back, conerollers in Tables 6.ile 
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CHAPTER 7 


CONCLUSIONS AND RECOMMENDATIONS FOR FUTURE WORK 


Vole CONCEILDUELONL OL Lhe Thesis 


The main contributions of this study are: 

i) Derivation of necessary and sufficient conditions 
for undisturbability, and the use of these and 
related results is the basis for the development 


of a simple systematic, design procedure; 


ii) A basis for direct analysis of the structure of 
systems which is easy to apply, and is useful from 


a control system design point of view; 


iii) A detailed evaluation of the new design tech- 
nique, by experimental application to the double 
effect evaporator, and by digital simulation of 
He Bes and 20 order models of two different binary 


aistiliationm columns. 
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iv) Development of a new procedure for designing 
State feedback controllers to assign closed-loop 
system eigenvalues and also minimise the effect 
of external disturbances on system outputs of 


interest; 


v) Development of a simple and constructive result 
for designing constant multivariable feedback 
controllers to make selected closed-loop eigen- 
values invariant to arbitrary variations in 


system parameters; 


vi) An algorithm to compute the invariant zeros of 


a system. 


7.2 Conclusions 


In this thesis the role of eigenvalues and 
eigenvectors in multivariable control systems design has 
been illustrated by relating the properties of undis- 
turbability and eigenvalue invariance to the structure of 
the open or closed-loop system eigenvector matrix. 

The concept of undisturbability has been formally 
defined and related to such concepts as uncontrollability, 
non-accessibility, and structural controllability. Nec- 


essary and sufficient conditions for undisturbability 
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have been expressed in terms of both the structure of 
the coefficient matrices of the state-space model and 
also in terms of the structure of the system eigenvector 
matrix. These results provide new insight into the 
structure of linear multivariable systems and 
provice a basis for the analysis and synthesis of 
regulatory controllers. Specifically, the results 
allow determination of the property of undisturbability 
of specific input-output pairs by a simple inspection 
of the zero entries in the appropriate matrices. Exist- 
ing design procedures for eigenvalue/eigenvector 
assignment have been reviewed and conditions for eigen- 
value/eigenvector assignment have been discussed here 
in terms of range space restrictions. The results on 
undisturbability, together with an existing algorithm 
for eigenvalue/eigenvector assignment lead to a simple 
and constructive design procedure for synthesizing 
multivariable controllers which provide disturbance 
localisation. An important characteristic of this de- 
sign procedure is that it makes use of the information 
contained in the model about the effect of disturbances, 
whereas many other design schemes for feedback control- 
lers do not. 

Experimental evaluation of feedback plus feed- 
forward controllers designed to produce undisturbability 


was carried out by application to a computer controlled, 
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pilot-plant evaporator. The design procedure was 
further evaluated by digital simulation of a 
and aon order models of two different distillation 
columns. These applications demonstrate the simplicity 
of the design procedure in its application to large 
order systems. The problems of simultaneous disturb- 
ance localisation with asymptotic setpoint tracking, 

and disturbance localisation using proportional plus 
integral feedback have also been investigated. 

A design procedure for minimising the effect of 
external disturbances on specific system outputs of 
interest and assigning eigenvalues in a closed-loop 
system, has been developed. This method can also be 
used in the design of full order observers for systems 
with unmeasurable disturbances. Experimental evaluation 
of the proposed minimisation scheme on the double effect 
evaporator gave good results and also demonstrated the 
practicality of the method. 

Simple constructive results for the design of 
constant multivariable feedback controllers, to make 
selected closed-loop eigenvalues invariant to unknown 
perturbations of arbitrary magnitude in system para- 
meters, are given. In general, these results show the 


advantages of having a quasi-triangular (or decomposed) 
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system structure either for undisturbability, or for 
making dominant eigenvalues invariant to arbitrary 
perturbations in specified system parameters. A discus- 
sion on the relationship between such concepts as invar- 
iant zeros, (A,B)-invariant subspaces, and parametric 
Sensitivity is also included. One outcome of this dis- 
cussion has been a simple algorithm for computing the 


invariant zeros of a system. 


7.3 Recommendations for Future Work 


As a result of work done in this study, several 
new research areas worthy of further consideration have 


been identified. Some of these are: 


eal) structural observability 


Lin [1] first introduced the idea of structural 
controllability in 1974. Since then Shields and Pearson 
[2], and Glover and Silverman [3] have extended this con- 
cept to also apply to multi-input systems. The duality 
of these results for investigating 'structural* observ- 
ability has not been formally established. Because of 
the algebraic nature of the results on undisturbability, 
duality between undisturbability and unobservability is 


not difficult to establish. When a one to one 
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correspondence between undisturbable modes ce state 
variables exists, the undisturbable modes in a dual or 
adjoint system are the unobservable modes or state 
Variables. In fact, the discussion on undisturbability 
and structural controllability can be immediately trans- 
lated by using duality to relate undisturbability to 
"structural' unobservability. However, the limitation 
here is that the dual of undisturbable subspace is 
smaller than or equal to the ‘structurally’ unobserv- 
able subspace, and therefore the complete class of 
structurally unobservable subspaces cannot be identified. 
Tt 258 not difiticult tojsee that: the concept of 
structural observability can be of great use in sensor 
location and measurement studies in the plant design 
stage (see remarks on the control scheme for the 20en 
order distillation, column. = Sections 4.5.2) -eeotiuctural 
observability considerations can also influence sensor 
location and measurement in systems that require state 
estimation or filtering, i.e. the object here would be 
to locate sensors such that the ‘(structural) observa- 


bility index' is maximised. 
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controllability (and 'observability') in a system 
requires very little information - i.e. only knowledge 
of the zero elements in the system. As a result such 
an analysis prior to the actual plant design and con- 
struction stage can be of help in avoiding potential 
pitfalis (e.g. in avoiding ineffective icontrol inputs) , 
or it can help identify and suggest modifications 
that may help the final control of the plant (e.g. 
in modifying the feed system to make a critical plant 
process undisturbable). Such an analytical technique 
needs to be systematised for use as an effective design 
tool. | 

TEeMsenotr difficult to see that the structural 
results on undisturbability etc. also apply to graph 
or node related structures, for example computer flow 
diagrams or critical path schedules. Consequently, 
it would be worthwhile to investigate the application and 
relation of these concepts to such areas as the design 
and development of large fault-tolerant (= undisturbable) 
computer programs, or in the planning of critical path 


scheduling in face of uncertain tasks. 


7.3.3 Disturbance minimisation as a design tool in 
the frequency domain context 


The absolute magnitude of the transmittance 
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between a disturbance and an output has been used as 

a Measure of disturbance/output coupling in the design 
procedure for disturbance minimisation in Chapter 5. 
Extension of this idea in the frequency domain context 
would consist of plotting a Nyquist array of the load 
transfer function matrix, G(s), and then computing a 
dynamic or constant precompensator, K(s), to minimise 
the effect of disturbances on specific outputs at a 
specific frequency range of interest. There is a close 
parallel between the approach suggested here and the 
procedure due to Hawkins [4] which is concerned with 
pseudo-diagonalisation for the Inverse Nyquist array 
method of Rosenbrock [5]. 

Obviously the final design procedure would. 
involve working simultaneously with Nyquist arrays of 
the process and the load transfer function matrices, 
G(s) and G(s), respectively. Even though the comput- 
ational requirements for such a procedure are significant 
(i.e. it would require use of an interactive computer 
graphics facility) there seems to be considerable scope 
for a design procedure that simultaneously considers 
the two important design requirements: regulation and 


setpoint following through a realisable controller, K(s). 
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GENS Eigenvalue/eigenvector assignment using output 


feedback 


In this thesis some of the difficulties assoc- 
lated with eigenvalue/eigenvector assignment using 
Output feedback control have been identified. It 
appears from the work carried out here that such a 
direct procedure is prohibitively complex. An alterna- 
tive procedure worthy of further consideration is the 
control law reduction techniques of Wilson [6] and 
Bengtsson [7] or the modification of these methods such that 
the reduced order output feedback control law guarantees 
system stability and also approximately retains the 
favourable eigenvector structure of the closed-loop 


system matrix (with state feedback). Development of 
new or modified control law reduction techniques, that 
would preserve desired structural characteristics is a 


worthwhile area for future work. 
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NOTATION 


General: 


The field of real numbers and the field of 
complex numbers are denoted by R and ¢ respectively. 
The n-dimensional vector space over the real numbers 
is written R”, and over the complex numbers fee Upper 
case letters, e.g. A,B,C,D are used for matrices 
(linear maps). The range space of a matrix D is denoted 
by the corresponding script letter D (or D = Image D). 
Script letters, e.g., V, W are also used to denote linear 
subspaces. The distinction between range spaces and 
subspaces is clear from the context. The null space of 
a matrix C is written as Kernel (C) or sometimes as 
Ker(C). Lower case letters with or without subscripts, 
e.g. b., oy Les denote vectors. Sometimes lower case 
letters also denote scalars. The distinction between 
vectors and scalars is clear from the context. The 
symbol, +, will be used for summation of linear subspaces, 
and the symbol @, for the direct sum. o(:) is used to 
denote the spectrum or eigenvalues of the matrix (argu- 
ment). (A/W) denotes the matrix representation of the 


restriction of A to W and X/W denotes the quotient space 


OF X modulo W. 
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state coefficient matrix 

Control ior input coefficient matrix 
Dabti tions Of ec 

pseudoinverse of B 

output coetficient matrix 

disturbance or load coefficient matrix 
partition of a transformation matrix 
transfer function matrix 

closed-loop system matrix 

identity matrix 

canonical state matrix 

state feedback control matrix 
feedforward control matrix 

number of undisturbable state variables 
disturbance (or load) coefficient matrix 
for the closed-loop system with feedforward 
control 

Vee column of L 

dimension of output vector 

dimension of state vector 

dimension of disturbance vector 
dimension of control vector 


Laplace operator 
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woeeit Moun(continued) 


u 


vi 


control*vector, oxl 

closed-loop system reciprocal (or left) 
eigenvector matrix esa 

closed-loop system eigenvector matrix 

state vector, nxl 

state response of the closed-loop system to 
disturbance 

output vector mxl 


output response of the closed-loop system to 


disturbance 


Greek Letters 


p 


state coefficient matrix for the discrete 
system 

control coefficient matrix for the discrete 
system 

disturbance coefficient matrix for the 
discrete system 

nxn diagonal matrix whose diagonal elements 
are the eigenvalues of H 


disturbance vector, qxl 
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Superscripts 


FF feedforward 

—1 matrix inverse 

Ak transpose of a vector or matrix 
ar transform of a matrix 

ts pseudoinverse of a matrix 
ie oy 4. OF, 


} refer to partitions of a matrix 
eee cree 


Subscripts 
. component of a vector or row of a matrix 
5 component of a vector or column of a matrix 


Notation for Table 4,1 and Computer Plots of Evaporator 


Response 


Each figure caption includes a string of char- 
acters delimited by parenthesis and/or semicolons which 
completely define the design and application conditions. 

The general form and order of this information is 
as follows: 

(run number) 
(controller type (Table identification); model order; 


CONtCrOl type; nun, cOndtt1ons):. 


The codes used for each of these specifications 


are defined in the following sections: 
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Run number 


DL - disturbance localisation run 

ML = mudltulooperun 

OPT — (optimale contra run 

EA - eigenvalue assignment run 

MC “Bamodal COnLerel nun 

P - proportional feedback control 
P+FF - proportional feedback plus feedforward control 
P+SP - proportional feedback plus setpoint control 


Controller type 


SuDSCcLript 
ae - disturbance localisation controller 
mL = multiloop controller 
opt - optimal controller 
Superscript 
FF - feedforward controller 
SP — setpoint controller 


Model order 


MO = 3 Signifies controller design based on the oa 


order model 
MO = 5 Signifies controller design based on the oa 


order model 
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Control type 


FB - unless specified this denotes state feedback 
FF - feedforward control of feed flowrate 
SP - setpoint control 


Run conditions 


B20 7c Cr — denotes —2 037 step in Cr followed by a 130% 
step in CF. 
E20 cor - denotes +20% step in F followed by a -203 


step in F. 
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APPENDIX A 
PRELIMINARIES 


Ae > Matrix Pseudoinverse (ci. Greville [A.1]): 


The matrix inverse is defined for square matrices 
Only; it is used in the solution of sets of linear 
equations of the form Mx=y, so that x=M y. In such 
linear equations the number of unknowns is equal to 
the number of equations. For nonsquare matrices used 
to describe systems of equations where the number of 
equations is not equal to the number of unknowns, the 
equivalent operator is the pseudoinverse. If a matrix 
B has more rows than columns, the pseudoinverse is 


defined as 


B* = (B’p) tpt 
provided (BB) is nonsingular. In the solution of 
linear equations, the set of equations, Bx=y are the 
so-called overdetermined case - where there are more 
equations than unknowns. The resulting solution 
x28 y is best in a least-squares sense. If a matrix, 
say C has more columns than rows, its pseudoinverse is 


then defined as: 
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CGC a= CGE 


This corresponds to the undetermined case - there are 


fewer equations than unknowns. 


A.2 A-Invariant Subspace  (Wonham and Morse [A.2]) 


A subspace VC X is said to be A-tnvartant if 


vel ->Av ec V, i.e. AV GC V.. The restricted map A:>V 


is defined so that Av = Av for all ve V. The restric- 


tion of A to V is denoted by A/V. Some examples of 


A-invariant subspaces are: 


a) 


ie 


se Oi 


Lf Wi is the ee eigenvector of A, i.e. 
ALi = AyWar then clearly the subspace 
Y= {Wy 1Woe eet w,} is A-invariant. Further- 


ar ceer7W; are linearly independent 


then, dimension (UV)! = 1, 


more if Wi9W 


Let denote the subspace of controllable 


states. Ro is given by: 
Ro = B+ AB + AB toh eweuenome ct aa 


By Caley-Hamilton, ARS SG Roe L.e@. Ris 


A-invariant. 


If 86 denotes the subspace of unobservable state, 


where 9 is given by: 
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n sis 
6 = fy) (Kernel (ca* 1)) 
= 


then by Caley-Hamilton A6 C 6, i.e. 6 is also 


A-invariant. 


Such subspaces play an important role in deter- 
mining the internal structure of linear dynamical systems. 
For example if Ro @ So See TiC ARC he then it is 
possible to choose a basis (see Lemma 3.1, in this Appen- 
dix) such that o(A) = o(A/R) correspond to the eigen- 
values of the controllable modes and o(A/R) correspond 
to the eigenvalues of the uncontrollable modes (cf. 


Kalman Canonical form, [A.3]). 


A.3 (A,B)-Invariant Subspace 


The idea of an A-invariant subspace discussed 
earlier can be generalized to take into account the 
effect of state feedback as follows. For a state 
feedback controller, K, the closed-loop system dynamics 


are described by the equation: 
(tw eeee BR) x(t) 


A subspace VC X is called (A,B)-tnvartvant if it is 
(A + BK)-invariant, that is (A + BK)V C V. One import- 


ant point about V being (A + BK)-invariant is the 
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following: There exists a K such that (A + BK)VC V 
if and only if AV C.B + V. One obvious example of 
(A,B)-invariant subspace is the space spanned by the 
closed-loop system eigenvectors. 

The main usefulness of (A,B)-invariant sub- 
spaces arises out of the fact that many control problems 
can be cast eee of (A,B)-invariant subspaces. 
Once this is done, the solution for K can be obtained 
Since (A,B)-invariant subspaces can be made (A + BK)- 
invariant by a suitable choice of K. For example if V 
has been constructed to satisfy some conditions then a 


Suitable K can be computed as follows: 
AVC V+B > AV= [VjBIR, R = [VIB] “Av 


where V is now a basis matrix of V, i.e. the columns 

of V are linearly independent and they span V, and the 
* 

superscript denotes the pseudoinverse. If now 


* : A 
[ViB] is partitioned as 


ate “i 
vere | = 22] 


so that the dimensions of Ly and L, are comparable with 


those of V and B, then 


L 
AY =e VeIBa = AV = ViL.,AV + BLAAV 
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The above equation can be written as: 


(A - BLA) V = VL,AV 


comparing this to the definition for V as: (A + BK)Vc JV 


yields, K as shown below: 


Further discussion on this and related topics can be 
found in Wonham [A.4], Basile and Marro [A.5] and 


Morse [A.6]. 


Proof of Lemma 3.1 
Lemma 3.1 


Let M be a nonempty set of nxn matrices in 
ee if W is an invariant subspace (WC R") under M ¢ M 
and is of dimension (n-k) over R, then there exists a 
basis of R” over R such that every M ce M in this basis 


is of the form 


M 0 
; Cant) 
Roe ae 
where M, € RKXK ong My € p(n-k)x(n-k). Furthermore, in 


this basis M, is the matrix representation of the 


4 


restriction of M to Weand M- is the restriction of M to 


1 
R/W. 
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Proof: If dimension (X)=n then let X = W@Z. Let 


(Zy, Zon seer ) be a basis’ for Z and (Weaye Wea! ae 


ee be a basis for W. Then (ZyrZor seer Zr Weiye 


Wear seer w,) is sa basis: tor XX. Dt as clear that 
k 
p Aas um Atak tly lt I oak) Pa Pe Al (A.2) 
: aerate 
i=1 
and that 
n 
weWrwe= Ff OW (AGS) 
i=k+1 


Thus when Zs and Ww, are operated on by M, it is easy 


to see that: 


k 
MZeee=) Sy aetiee ee » TM. Ws (Ac4) 


k 
and Mie = oa ee 2 oe DS m..W. (A.5) 
n Ww ce ata feng? 1d 4 
In eqn. (A.5) the elements Mea are the elements of the 
matrix representation of M to W. Because W is 
M-invariant, in) eqn.) (A2>)), ie =\0 for 1 = 1,k and 
j = k+1,n(>M, = 0). This means that M in the basis 


(Z)1Zoe oer Zp rWe agrees W) is of the form: 
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and that M, is the matrix representation of the restric- 
£10n Of M-toml.. PAdse [M, M,]7 is the matrix representa- 
tion of the restriction of M to Z and since W represents 
the zero element in the quotient space ROW, My is the 


matrix representation of the restriction to RO /Ws or M, 


is the matrix of the linear map on RD /W induced by M (in 


the basis { (Zz, te We Saeeeats W) oe & vats 
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APPENDIX B 


EIGENVALUE/EIGENVECTOR ASSIGNMENT USING 
OUTPUT FEEDBACK CONTROL 


This appendix is in support of the material 
covered in Chapter 2. The following discussion outlines 
the algebraic manipulations required to arrive at an 
algorithm for eigenvalue/eigenvector assignment using 
Output feedback. It also points out the difficulties 
associated with the proposed algorithm. 

Consider a linear nth order multivariable 


system with r inputs and m outputs as given by 


x(t), =) Ax(t) + Bu(t) 


(B49) 
Vit man Cx t) 
and subject to the output feedback control law 
UNC SRV tema Ke CXtt) (BZ) 
The resulting closed-loop system matrix is given by: 
(A + BKC) = WIV (Bec) 
Partitioning eqn. (B.3) gives 
ae Sad ra, | Kic.ic een a Vocea ail vee ey 
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Assuming r = 


m, then in the above equation A 


po 2 Way? 


J, and Vay are mxm matrices; Anos Woo J and Vo9 are (n-m)x 


(n-m) matrices; K is mxm and B, and C 


] 


are assumed to be non- 


Singular, obtained if necessary by reordering of the state 


Variables. 
From-eqn. (B8.2:)< 
ee ee Te ee 221 (B.4) 
Roy + BoKC) = WoqJqVa71 + WoadoVoy (B.5) 
and hence 
7 4 | 
By = (Wy pJyVaq + Wy odoVoq - Ay )Cy (B.6) 
, nS, 
and Bok = (Woydq¥1q + WoadaVoq ~ Aq Ie (B.7) 
Also 
Ago + Boke, = Wo1I Vio + WoodoVon (B.8) 
and Ary + ByKCy = Ward Va5 + Wy odoVon (B.9) 
Substituting (B.7) into (8-3) gives: 
a1 ©. 
Rog + (Wo154Vq_ + WopJoVoq - Ao )Cy Co = WoydVyo + WoodoVoo 
ele . ...(B.10) 
and substituting (B.6) into (B.9) gives 
=a : 
Aro + (WyyJqVyq + WyadoVng - Ayq Oy Co = WyqdqVy2 * Wy adoVon 
cp ete eles | Tecaniinive Gan Dhaene -+-(B.11) 
Aaa Ea ice: Pay ee Yowell 
Va ayia + Va aeta¥er® - Vapaaa® = Yarra “ie * 
Vay atovon - VaqAr2 (B.12) 


using the identity Vay ay “5 Va Woy = I gives 


JV 579 - My yA 70 + VaqAqe = JyYao = My oWoglyMi12 - Vay otovor® - 


V U54 


pa oy 12 


+ VN one 


(Be 13) 
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Pre-multiplying (B.10) by Vin gives 


ViaAge ~ Yoho 2 = VyoWoqJ Vy + VyaWoodaVon - VyaWoq91V419 - Vy aWoodaVo,Q 


(B. 14) 
Using the identity Viq"15 = “V1 Woo gives 
Vyahoe ~ Vyaho1® = VyaWoyty Yio - Vay MyataVoo - YyoWoqeyy 118 + Vay otovey9 
(Be15) 
substituting (B.15) into (B.13)° gives 
JyVy2 ~ VyalAge - ApqQ) = 94¥q4 + ¥yqCAy = Ay 79) ee a 
Let L = Ayo - Ay, and R= Ags - ApqQ 
Then (B.16) becomes 
P42 pe SesNig cakopins tans aed 


Premultiplying eqn. (B.10) by Vo9 gives 
Voahoo ~ Vooha19 = VooWo 5 Vyo + VooWoodoVan - VooWoq91V 119 - VooWoodaVo79 


(Biers) 
Using the identity Voq hay = Volo) gives 


VooAoo - VooAo19 = -Vo Wy q9qVyo + VogtaodoVan + VoqWy IV 418 - VooWoodoVny@ 


(B.19) 
Premultiplying eqn. (B.11) by Vol gives 
a tee 2) 1 eet elie ene elie ie 2 ese ee 
(B. 20) 
Using the identity VoqMa0 = [- Vo ooo gives 
VorAy2 ~ Yoyhy19 * Jo¥o7® = Yovan = Vopty yyy Mia ~ VoaWoota¥oo - 


Vo Wy 191712 * VoaWon4aVo79 eet 


Jae 


(20.8) a | a { af ; 


asviy (f 1.9) odnt neat 


ol ci) (Dy ph = gehhigy + Pratt * (Dru = gghbg - OT is ; 
Gegh ~ sch * Abie Dep = oA = 188 
eamooed (af.6) 
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substituting (B.19) into (B.21) gives 


UeVie5, =v 


phee ee Paaee 71 Memeo eh aioe Te aber 


Using expressions for L and R gives: 


oR = JoV51Q + VoL CeIeASD, 


From eqn. (B.4) 


ae 
Ke Baa ot ov ay = Aaa? aa 
and from eqn. (B.9) 
Sees 
Ky. = By! (Wy yJqVpq + Wy oIVon = Aqy) (B.25) 
Substitutinge (6.24) initOm ea.o) 
u i 
Let S = BB, 
Roy + SWayJy¥ yy + SWyadoVoy ~ SAqq = WoqdqVq7 + WoodoVoy a 
Substituting (B25) ainto WiB.i3)) 
Ae sie J.V (B.27) 


22 TW’ 12 * SWyatovos - SAqo = WoqI Vio + WoodoVos 


Post-multiplying (B.26) by Way ands (Bea7y) by Woy gives: 


Roy Wyz > SAq aq + SW CE - VyoWog) = Woy d (CE - Vy alloy) + WoodoVoq Way 
Pst die (B.28) 
which on simplification gives: 
PoqWay ~ SAaqWay + SWyqdy > Woqdy = SW JyVy ably - Woqdq Vy oWoy 
* WooDVor¥yy ~ SMa 2IoVoq (B.29) 
and 
AooWo, ~ SAqoWo, = WoqdyVyabloy + WoadaVooWoy - SWy Iq Vy obo, 
Fa (B.30) 


cee cee 
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using the identify V = -\ 


aa 22Hoy 


AoW, ~ SAyoWoy = ~SWyy9qVyoWoy + WoqIqVyoWoy - WoodVoqhay 


ss SW) oJ oVoqWa, (B 
substituting (B.31) into (B.29) gives: 
Ge) eadial ls nna le ieeete hen enela reser al 
that-is 
Way, ~ (Aga - SAyo)Woq = Sih yd, + (Ap) - SAY) Way ‘Be 
Letting T = Ao, SA, 4 
and P= Ayo - SAy 5 
Eqn. (B.32) can be rewritten as: 
Woy - PWoy = SW, 49; + ™, (eke 


In the same manner post multiplying (B.26) by Wao and 


(Bered ie by Woo gives: 


W Jo - PW TW + SWaou 


epee) i272 


Equations (B.17), B.23),(B.33) and (B.34) are relevant 
for the design procedure. Inspection of eqns. (B.7) and 


(B.33) shows that these two equations have to be solved 


G35 
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oo) 


32) 


33) 


34) 


iteratively with the identity: Va qWay +: Va Woy = alae Pal Nast; (1.5: 


having made an arbitrary choice of Way and Vaqe Woy and Vio 
are evaluated from eqns. (B.17) and (B.33) respectively, and 
a check is made to see if the identity Vaqhay at Vy Woy = | 


is Satisfied; if not, new choices are made for initial 
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values of Way and V and the procedure repeated. Clearly 


1] 
convergence cannot be guaranteed for such a procedure 

and it becomes at once evident that the procedure can be 
quite lengthy. Once an appropriate choice of Way and V4y 
is made and Woy and Vio have been computed such that the 
required identity is satisfied, then a choice of Wi and 
Voy is made and eqns. (B.23) and (B.34) are solved for 
Vio and Woos respectively. Now a test is made to deter- 
mine whether the following identities are satisfied: 


V,4W TV eel = 0, V,,W aan | = 0 and 


Wee 12% 90 1811 22Woy 


VoqWi - Vo oW oo = I. If not, the procedure is repeated 
for different choices of Wi and Voq° It is tnus clear 
that the algorithm can be extremely lengthy and that 
there are considerable difficulties involved with the 
design procedure for eigenvalue/eigenvector assignment 
using output feedback. 

Note that the procedure outlined briefly at the 
end of Reference Blmay at first appear to be slightly 
different than the one above. However, both are equiv- 
alent. BY Initially choosing only Woy and Vay arbitrarily 
then computing Woy and Vio (from eqns. (B.17) and (6 .33),) 
and finally checking to see if the identity Vaya + 
Va oWoy = I holds, the procedure as proposed above should 


be slightly easier to follow. Note that convergence, 
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i.e. getting the correct values of V and W such that 
Vie = Weholdsmas: well aceeqns, (8. 1/)..(823)., (Base) and 
(B.34), is not guaranteed. 

This appendix essentially points out the lengthy 
and difficult design procedure that follows as a result 
of difficulties imposed by the rather restrictive set of 
necessary and sufficient conditions for eigenvalue/ 
eigenvector assignment using output feedback as outlined 


in Chapter 2. 


REFERENCE FOR APPENDIASE 


Penola saeSe ines ok NSe aL). Gis an Gp SA UOT O56). Lie, 
'Eigenvalue/eigenvector assignment for multivar- 
iable systems and further results for output 
Feedback control , Electronics’ tett.; 1975, 11, 
pp. 388-389. 
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APPENDIX C 


PILOT PLANT EVAPORATOR AND ITS MODELS 


Cll she Evaporator 


The pilot plant evaporator used in this work is a 
double effect unit with the two effects operating in ser- 
ies. This pilot plant has been described in detail by 
Fisher and Seborg [C.1]. The major pieces of process 
equipment are shown in the schematic diagram in Figure C.l. 
The control loops shown in Figure C.1 represent the multi- 
loop control scheme applied to the evaporator in previous 
Sscudies [C.1].. 

The first effect has natural circulation through 
Leo LS anche lond,. 3/74) inch O.Ds tubes. Ib as heated by 
process steam. The second effect is a long tube vertical 
Unit which Ys run an, ies torced circulation mode. Lt 
has three, six foot long, one inch O.D. tubes. It is 
operated at a lower pressure than the first effect and is 
heated by the vapour produced in the first effect. 

The evaporator is fully instrumented and can be 
controlled by either Foxboro electronic controllers or 
under Direct Digital Control (DDC) from an IBM 1800 Data 
Acquisition and Control Computer operating under MPX. 


Multiloop DDC can be applied directly using the computer 


267 


” P i 7 iz - 
7 
nynetOM ane awh fem ait oad |! ye 
poe - oa ~ 
| - 
ai Avow wide a2 sehr ros sHOge re sree sis pest 


tes ei paiteieen t0sTtte wt ais daiw shout 4 a ands ua 


) Eigjsab mei tagiegeab nasd nen same’ barat set 
baat er alt 


TI ol: meee vitamedoe eds al Labial ex dawg 2 


sua 1o S|eaoRk - van eat .tt.3) puodse 


mf ernt. ots SAatTes bd 9 etupl” ot crwoMer eqoot tox4n00 ® sl r 
Ap 


L 


Iveta. oo. to¢bY%omwsve sd3o- of beclqgs ne foaites JOO 


mid ‘Zhe 9 ) aetbuye 
ae aa 


i. 7 
dppotda np hts lvotio: Desuvem aan ssotte aexk? pit 


qi bectaot en at’ ohitith 0 atteat Wie: spat doat 8 


(soiaxesy edut onel oeel ° eF%s. Bro ae oan ms y ance octq | 
77 > .sbom soitsivouts Beret aoe fil i ‘ . os «sin yo 
escaut 10.0 feed Sno eet ‘hoe nim 4 a eel 
ae! foe tootte agaxk? ong nade ety cts woos oan bs 28 tea 
Jhobie seultesns wt tio fy 80 Har 9 ads. Yo hodeod 

at aso Soo hedascatenl| ethos adi Toss GSTS edt 
fo es2slioisnaas sun we oa dh chal x ee obs eit 743000 
fexsno 2 ae, Jipia go a tebs 
-X4M wober ade ay if eT sy oe | 


adnd OBL MAT os aoe, ‘(oa 


a ylise ott pbs Nags 


aii ‘ sence ous pitt 


268 


control package and advanced control schemes by user 
written programs which utilize a set of svstem programs 


to interface between the user and system control programs. 


C.2 The Evaporator Model 


The complete development of the double-effect 
evaporator model has been presented by Newell [C.2] who 
derived fifth and tenth order non-linear models. Based 
on these models Wilson [C.3] calculated discrete five- 
state and three-state, linear, time-invariant models 
using a linearization procedure and Marshall's model 
reduction method, 

The linearized models in the discrete form can 


be represented by: 


x(n+1) = Ox(n) + Au(n) + Od(n) (Cet) 
and 


y(n) = C x(n) (Cr 2) 


The elements of the vectors x, u, d, y are defined 


as normalized perturbation variables: 


W, - W 
mle _i_ iss (cis) 
W 
lss 
where Nites is the normal steady state value of Wi- 
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The vectors xX, u, Gd and y for the fifth order 
discrete model are defined in Table C.1. The coefficient 
matrices of this discrete time model, with a 64 second 
sampling interval, are shown in Table C.2. 

For the third order discrete model the state 


vector x is given by: 


PISeuienoeeD) © (Cc. 4) 


The elements of the state vector are defined in Table C.l. 
The vectors u, d and y are equal to the ones defined in 
Table C.1l. The coefficient matrices for this model, 


with a 64 second sampling interval, are shown in Table 


eos 
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DESCRIPTION OF THE EVAPORATOR VARIABLES 


State Vector, x 


x) = [W1, Cl, HI, W2, €21 


W1 
Cl 
H] 
W2 
C2 


First effect holdup 

First effect concentration 
First effect enthalpy 
Second effect holdup 


Second effect concentration 


Control Vector, u 


ive 
U = 


S 
Bl 
B2 


[Sal eee, 
Steam flow 
First effect bottoms flow 


Second effect bottoms flow 


Disturbance Vector, d 


d 
F 
CF 
HF 


de 


[F, CF, HF] 
Feed flow 
Feed concentration 


Feed enthalpy 


Qutput Vector, y 


rT 
y = 


[W1, W2, C2] 


Normal Steady State Value 


2525 iD. 
4.59% glycol 
18952-8107 \be 

AT 5rd. 
10.11% glycol 


2.0 1b./min. 
3.485 1b./min. 
1-587 Vbe/maine 


Se Oni be/min. 
3.2% glycol 
156.9 BTU/Ib. 
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PIFTH ORDER 


-0.0119 
0.0116 
0.0116 

-0.0138 
0.0137 


THIRD ORDER DISCRETE EVAPORATOR MODEL (T 


o — 
can) 


0 
0.0848 
-0.0432 


0 
0 
0.0400 


TABLE C.2 


DISCRETE EVAPORATOR MODEL 


-0.0008 -0.0912 
Qyeece 0.0871 
-0.0042 0.4376 
-0.0009 -0.1052 
0.039) 0.1048 
0 
0 
0 ) 
-0.0406 
0 
0 0 
0 ] 
0 0 
TABLE C.3 


-0.0135 
-0.0157 
OUZIg 


0. 1182 
-0.0351 
-0.0136 

0.0012 

0.0019 


-0.0326 
0.0378 
0.0529 


=) G4 SEC.) 


-0.0050 
0.0049 
0.0662 

-0.0058 
0.0058 


= 64 isec. ) 


-0.0811 
0.0854 
-0.0442 
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-0.0406 
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APPENDIX D 


ON INVARIANT ZEROS, INVARIANT SUBSPACES 


AND PARAMETRIC SENSITIVITY 
Dek «6Llntroduction 


This Appendix extends the discussion of Chapters 
3 and 6 and is particularly concerned with relating the 
following three relatively new concepts: (i) invariant 
zeros (as defined by Bengtsson [D.1]), (ii) (A,B) -invar- 
iant subspaces [D.2] (see also Appendix A), and (iii) 
eigenvalue insensitivity to variations in system para- 
meters [D.3](see also Chapter 6). 

In [D.1] the invariant zeros of a MIMO system 
have been defined as those zerosof a system that are 
invariant under state feedback in the sense that the 
open-loop system denoted by S(A,B,C) and the closed-loop 
system S(A+BK,B,C) have the same set of zeros for a class 
of linear maps K. The concept of (A,B)-invariant subspace 
has played a key role in a number of control problems [D.2]. 
The most familiar one is the disturbance localisation 
problem. In the following discussion these two concepts 
will be related. It will be shown that when certain set 
of conditions are imposed on (A,B)-invariant subspaces, 


then a simple algorithm to compute the invariant zeros of 


* for the general definition of zeros of multivariable 
systems see Macfarlane and Karcanias [D.4]. 
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a (left) invertible system can be derived. In order to 
show a use for such an algorithm we will relate invariant 
zeros to insensitive eigenvalues. [In particular, it will 
be shown that the invariant zeros are the same as those 
eigenvalues which are insensitive or invariant to 
arbttrary parameter variations in specified columns of 
the A or K matrices. In contrast to results in the first 
part of Chapter 6 which are strictly valid only for eigen- 
value invariance or insensitivity to small, arbitrary 
perturbations in specified columns of A or K, it will be 
shown here that the analysis carried out here allows a 
designer to make certain eigenvalues invariant to unknown 
variations of arbitrary magnitude in specified columns of 


Anor Ke (Ci. Theorem 6.L). 


D.2 Invariant Zeros and Invariant Subspaces 


In [D.2] it is shown that the maximal (A,8)=in- 
variant subspace, Vien contained in Kernel C, can be ob- 


tained by the following sequence (see also discussion of 


(A,B) - invariant subspaces in Appendix A). 
V = Kernel C 
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Then let p be the first integer such that aah = eae 
then U™ = os This sequence converges after at most 


Vv steps, where v = dim (Kernel C). 
Note that to simplify the calculation of Y" we 
impose the following condition on Ce + B), i.e. assume 


that dimension (ye + 8) =n. Then 


v= Un A eee) 42 Bye ay bay Sarnel! ¢. 


For the case where k = r (i.e. the number of undisturb- 

able state variables or outputs (k) is equal to the number 

GOL input Or control variables (7)), the condition: 
dimension (V° + B) =n, 


implies the following: 


i) The first k rows of B must be linearly independent. 
ii) In addition dimension VO = n-k since C can be 


written as C = [I O], and in order to realise the 


k 
property (A + BEI oS VS Seis required that (n-k) 
of the eigenvectors of (A + BK) lie in ise. i.e. (n-k) 
of the closed-loop eigenvectors lie in Kernel C. 
Using the eigenvector/eigenvalue assignment technique 
(of Chapter 2) it becomes clear that for (n-k) of 

the closed-loop eigenvectors to lie in Kernel C, 

Wid =O. The eigenvalues corresponding to these 


eigenvectors are then given by: 
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and the diagonal elements of Ao are the eigenvalues 


of P. 


D.2.1 A New algorithm to compute invariant zeros Of a 


System 
Note that for the above case o(A + BK] V™) =o (P). 
From [D.1], invariant zeros of a left invertible system 
are o(A + BK|U™) Hence we have the following 


new algorithm to compute the invariant zeros of a system: 
If dimension (v° + 6B) =n cand kK =r, dee. the 
first k or r rows of B are! linearly independent then the 
invariant zeros of the system are given by the eigenvalues 
=—] 


of Oo B.B, Ai) where: }A,58 15 the (n-k) x (n-k) 


lower right matrixipantition of A and Alo is (k x n-k) 
Upper right matrix partition of A. By and Bo are the top 
feet (nNOtLerKu= tN (Or ek x ik) “and: bottom (i =k) x xs 


(or (n-k) x k) partitions of B respectively. 


Example D.1: 


Consider a system with the triple (A,B,C) given 


by: 


220 


A 
; 7 
(iyo = “yiek.+.Alo eies ooh etd 208 Suda’ * on : Ca 
nedaye alrttotevnat tte, 8 Bo sere: stnetiavnt atte): late . 7 


i) 


tine 


u7e 6 Jo 2ecey Sodtveval ott TGS oo ad opie wort > 
ett J 1ebreyshnw 4 lca it iat »@ Fo, ao | a T 3 a ra ne 
revlia aft Yd. MeVvip os 1 moacave odd Xo eo75s snsiaaval 


Hn) « (A+a) tego res ee® Ss .ernw (gp gee - 


Si WIG ofl ene 8 A Re varepanlie eA ueqqu 


7 


i ae 


wollot ad’ gvet ow apmell (yx Soa A) o Pn 


2 


i, 


et (ted: Sa ae 3 + M) notentemit 3t) 


a 


> 


Mio ek pA Bae A Ae notustieg wht Bat tapts ue _ 


s * (4 -onh) edged bout (ot et at tod Mf A 4 etcon) . «x 
Vis ee ee @ to enorst tis (t & Oem) aay 7 


a = ial an 


RT ea 


To rewrite the output matrix of the above system as 


Co pear ol O], use will be made of the transformation 
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and with first two rows of ce linearly independent in- 


variant zeros of the above system are 
g a i feo ep 
=o (A) BoB, A,) 


II 


o (-3+4) = {1.0} 


Example D.2 


Consider the SISO system with 


given by: 
-1 0 0 iL 
A= 0 =] 0 Bye [1 C= 
0 0 =2 Ai 


using the transformation (Cf. Section 


ie rae ee i 
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1 Bi 
ml Siero alee ee 
Bt 
Z 
An 
0 0 
Now ee =p ig eee lea ii) , also aim. le + 8) = 3 
0 ie 
andy since the first row (* = L)Por Bt is linearly inde- 


pendent the invariant zeros of the above system are 


y SRM hs ak ae ay 
= 0 (A, BoB, A.) 
-1 0 ul 
i.e. eigenvalues of - (GF =| 
0 -2 1 
-1 1 
a6 = {-l1, -1} 
Oia 


Note that these are precisely the eigenvalues which can- 
not be assigned arbitrarily in the design of a state 


feedback controller to make z4(=y,) undisturbable. 


a, oe 
, 
2 « 
cE me ¥) smi sate", 
py Lt. Saas 
7 ; 2 
‘ 7 7 
yizsenit ai "R tm (fo *) wos sand? \edtt 
¥ ¥ - 
16 jdeye eves edt Jo SOhtear Mes a ay t 
: ; 7 
ie i 
| | 7 . 
E | {4 sells 
AA - 
| Ht 0 
== | —| 
ty io 
ee Vel 
| 
- i la 
ea a1 ® 
“A650 daidw aovlavnspis saz tai oT 


stata = to toteeb + 


280 


D.3 Invariant Zeros and Parametric Insensitivity 


For the class of systems considered here, namely 
where C = [I) O], the first k rows of B are linearly 
independent and k = r, and hence Ue = /"' = Kernel Cr. Le 
is clear that in order to satisfy (A + BW ae V°, (n-k) 
of the closed-loop system eigenvectors must lie in 
Kernel C. In other words the first k elements of the 


(n-k) eigenvectors must be zero, i.e. OS. Bor 


ee 
such a class of systems there are (n-k) invariant zeros. 
For such cases it can be shown that: 
i) it is possible to make as many eigenvalues invariant 
or insensitive to parameter variations (in sperseeuese 


columns of A or K), as the number of invariant zeros 


of the system, and 


ii) the invariant or insensitive eigenvalues remain un- 
moved despite arbitrary and unknown parameter 
changes in specified columns of K and hence the A 
Matrix. Note also that these parameter variations 
need not necessarily be small. This point will be 
elariived shortly. 

For a left invertible system {invariant or insensitive 

eigenvalues} = o((A + BK) |V_)- Denote by K" the class of 
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where Ke is @ Variable k x kK matrix (assume K = x) 
and Ky is a k x (n-k) matrix of constants. In other 
words for any arbitrary vartatton in K. (i.e. the first 

k columns of K) or equivalently in the first k columns of 


A, the eigenvalues o(A + BK | V°) remain invariant or 


insensitive. This point will be illustrated by the follow- 


ing example: 


Example De3 


Consider the system S(A,B,C) where 
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The open-loop eigenvalues are 0, 1, and -l. If the 


measurements of x4 and Xo 


meters in the first and second columns of A are likely 


are unreliable or the para- 


to vary then we can design a K_ such that (A + BK) VC Viner 
with the resulting invariant eigenvalue of (-1) being 
completely insensitive to any arbttrary vartatton in the 
first two columns of K or A. The remaining two 


eigenvalues of (A + BK) can be assigned arbitrarily. 


Suppose then Wii is chosen as 
1 0 0 
Wi = and Wio = 
0 a 0 


(for invariance to first and second columns of A and kK) 
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Where B is the pseudoinverse of B. 
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and the resulting eigenvalues are: -l, -2.2 and -1.2, 
where the invariant eigenvalue is -l. Notice that the 
variations in K and A are of stgntficant magnitude and 
usually cannot be handled by conventional sensitivity 
techniques. 

Another result that follows easily from the above 
analysis is that the initial condition response of the 
above system for all x(o) € ve is invariant regardless 
of any variations in the first two columns of A or kK. 
This corresponds to trajectory tinvartance of all x(t) 
for all initial conditions x(o) which lie in Vo° (Cr. 
Remark 6.3). For example for the above case the response 
(opis x, (t), X~ (t) and X(t) is unchanged for x(o)2 = 
[0 () a] 4 at0, i.e. x(o) é Oe Also for b, rnb 
(ba is the jth column of B) the class of K_ such that 
(A + BK) V2 C V®° results in trajectory invariance. (This 


result is really related to disturbance localisation.) 


D.4 Concluding Remarks 


In this Appendix the concepts of invariant zeros, 
invariant subspaces and parametric insensitivity have 
been related. The three illustrative examples show the 
use of the proposed algorithm for calculation of invar- 
iant zeros of a system. It is also shown how the invar- 


iant eigenvalues are the same as invariant zeros of the 
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system, and hence are insensitive to arbitrarily large 
and unknown parameter variations in specified columns of 
the A or K matrices. For single-input single-output 
systems the relationship between invariant zeros and 
parametric insensitivity can also be understood from a 
root-loci point of view - that is because of a pole-zero 
cancellation (invariant eigenvalue - invariant zero can- 
cellation) variation in certain elements of K does not 


move the invariant eigenvalue. 
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